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ABSTRACT 


This  report  provides  the  documentation  for  a  set  of  computer 
applications  for  the  evaluation  of  antenna  parameters.  The 
applications  are  written  for  the  Mathcad  personal  computer 
software  for  various  antenna  types  listed  in  the  thesis  index. 
Antenna  dimensions  and,  in  some  cases,  ground  parameters  are  the 
only  required  inputs  for  each  application.  No  new  antenna 
parameter  equations  were  developed  as  a  part  of  this  research. 

The  chapters  of  this  thesis  are  intended  to  provide  Mathcad 
antenna  application  users  with  the  background  information 
necessary  to  readily  use  and  interpret  the  software  for  each 
antenna  type.  Appendices  are  provided  with  examples  of  each 
antenna  application.  Each  application  has  an  introductory 
paragraph  and  a  table  of  required  inputs. 

The  Mathcad  software  provides  various  numerical  outputs  and 
performance  predictions,  as  well  as  a  graphical  representation  of 
radiation  patterns  in  the  far-field.  Mathcad  application  results 
are  consistent  with  the  predictions  of  applicable  publications, 
as  well  as  other  antenna  numerical  analysis  programs. 
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BACKGROUND  AND  P0RP08E 

This  report  and  associated  Mathcad  computer  software  are 
submitted  in  partial  fulfillment  of  the  thesis  requirements  for 
the  degree  of  Master  of  Science  in  Electrical  Engineering  from 
the  Naval  Postgraduate  School  in  Monterey,  CA. 

The  thesis  requirement  was  generated  by  a  statement  of  work 
from  the  Naval  Maritime  Intelligence  Center  (NAVMARINTCEN)  such 
that  any  IBM  compatible  personnel  computer  with  MS-DOS  version 
3.2  or  higher  and  a  math  coprocessor  could  run  Mathcad  software 
applications  to  analyze  the  parameters  of  different  antenna  types 
requested  by  NAVMARINTCEN. 

Required  user  inputs  to  the  applications  are  limited  to 
antenna  dimensions  and  ground  data,  although  in  some  cases  other 
data  may  be  estimated  to  provide  further  insight  into  the 
antenna's  performance.  The  Mathcad  applications  provide  various 
performance  predictions  as  well  as  a  graphical  representation  of 
the  antenna's  far-field  radiation  pattern.  The  corresponding 
thesis  chapter  furnishes  the  application  user  all  the  necessary 
background  information  needed  to  interpret  the  program's  formulas 
and  displays,  thereby  allowing  NAVMARINTCEN  to  interpret  the 
capabilities  and  limitations  of  antennas  of  interest. 

Dietrich  [Ref  1.]  completed  the  first  portion  of  this 
project.  This  thesis  will  be  the  second  in  a  series  of  three 
reports  intended  to  fulfill  the  NAVMARINTCEN  statement  of  work. 
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II.  INTROPnCTIQW 


When  a  foreign  country  develops  a  new  conununi  cat  ions  or 
radar  system  there  are  many  reasons  why  various  United  States 
agencies  may  want  to  be  appraised  of  the  new  equipment's 
capabilities  and  limitations.  Indeed,  if  the  country  is  hostile 
to  the  United  States,  the  need  for  rapid  threat  analysis  can  be 
urgent.  Unfortunately,  without  some  human  intelligence  or  other 
highly  classified  source  data,  input  to  any  threat  analysis  is 
constrained  to  dimensional  information  gained  from  photographs  of 
the  equipment's  antennas.  In  the  past,  intelligence  agencies 
analyzing  each  new  system  on  a  case  by  case  basis  found  this 
process  to  be  very  slow,  tedious,  and  man  power  intensive. 

With  the  advent  of  powerful  personal  computers  and  the 
availability  of  sophisticated  mathematics  software,  antenna 
analysis  using  data  obtained  from  photographic  intelligence  may 
now  be  achieved  in  a  rapid  manner.  The  goal  of  this  report  is  to 
document  the  software  developed  to  accomplish  this  type  of 
performance  appraisal.  With  this  report  and  its  programs, 
NAVMARINTCEN  is  supplied  with  a  user  friendly  tool  to  aid  in 
their  task  of  antenna  system  evaluation. 

It  should  be  recognized  that  computer  analysis  of  antenna 
parameters  has  its  limitations.  For  example,  it  is  impossible  to 
account  for  the  effects  of  adjacent  structures  on  far-field 
radiation  patterns  and  still  keep  the  Mathcad  applications 
moderately  simple.  In  addition,  without  knowledge  of  parameters 
such  as  feed  line  characteristic  impedance  and  antenna  materials. 
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it  is  impossible  to  precisely  assess  the  efficiency,  gain,  and 
radiated  power  of  any  antenna.  Nevertheless,  in  most  cases  this 
report  should  provide  the  tools  necessary  to  gain  an  excellent 
initial  insight  into  the  capabilities  of  systems  which  use  the 
antenna  types  considered  herein. 

Each  chapter  of  this  report  reviews  a  specific  type  of 
antenna  and  is  written  as  a  comprehensive  reference  for  the 
software.  Copies  of  each  application  are  included  as  appendices 
to  provide  the  user  with  a  printed  illustration  of  the  software. 
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III.  THE  HELICAL  AMTENNA 


The  helical  antenna  is  a  wideband,  highly  directional  device 
when  operated  in  the  axial  radiation  mode.  It  is  commonly  used 
in  satellite  and  communications  systems.  To  implement  the 
Mathcad  software  for  the  helical  antenna  the  following  dimensions 
are  needed: 

D  =  diameter  of  helix  (center  to  center  of  the  conductor) 

S  =  turn  spacing  (center  to  center  of  the  conductor) 

L  =  length  of  one  turn 
n  =  number  of  turns 
d  =  diameter  of  helix  conductor 
C  =  circumference  of  the  helix  =  ttD* 
a  =  pitch  angle  =  tan'^ (S/ttd) * 

The  first  four  dimensions  are  application  inputs,  the  *  indicates 
that  the  remaining  two  parameters  are  calculated  by  the  helical 
antenna  Mathcad  application.  The  helical  antenna  geometric 
relationships  are  illustrated  in  Figures  3.1  and  3.2. 
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D 


=  Distance  along 
conductor  betweei 
arrows 
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FIGURE  3.2  Pitch  Angle  of  a  Helix 

The  axial  mode  helical  antenna  has  a  highly  directive  main 
lobe,  negligible  mutual  impedance  with  adjacent  antennas,  a  low 
voltage  standing  wave  ratio  (VSWR) ,  and  a  resistive  input 
impedance  if  the  following  are  conditions  met  [Ref  2:  pp.  277- 
288  ]  : 

.8iCj^l.l5  {wavelengths)  (3.1) 

n>3  {turns)  (3.2) 


I2°sa^l4’  {degrees) 


(3.3) 


Subscripts  containing  (X)  indicate  the  dimension  in  wavelengths. 

Assuming  (3.1)  -  (3.3)  are  satisfied,  one  can  estimate 
directivity  (Dq)  as  follows: 

D.=12C^^nS^  ( dimensionless)  (3.4) 

For  a  long  helix  (nS^  >  1)  ,  the  relative  phase  velocity  of 
the  traveling  wave  (p)  is  the  key  variable  for  calculating  far- 
field  radiation  patterns  and  associated  parameters.  Although 
several  equations  can  be  used  for  determining  relative  phase 
velocity,  the  one  which  most  closely  matches  measured  results  is 
[Ref  2:  pp.  288-300]: 

L, 

p=— - f—. — -  {dimensionless)  (3.5) 

(  1  /  277) 

In  (3.5),  (m)  corresponds  to  the  transmission  mode  number  of  the 
antenna.  The  transmission  mode  is  a  term  used  to  describe  the 
manner  in  which  an  electromagnetic  wave  propagates  down  the 
helix.  The  number  assigned  to  a  given  transmission  mode  (T„)  is 
an  integer.  When  m  =  0  the  helix  radiates  in  what  is  termed  the 
normal  mode,  since  the  main  lobe  is  perpendicular  to  the  axis  of 
the  helix.  In  some  texts  a  helical  antenna  which  radiates  in  the 
normal  mode  is  called  an  electrically  small  antenna.  The  normal 
mode  is  not  commonly  used  and  will  not  be  covered  further  in  this 
report.  The  Mathcad  applications  analyze  only  the  non-zero 
transmission  modes  of  a  given  helical  antenna. 

The  mode  of  a  helix  is  determined  by  its  physical  size,  with 
higher  modes  corresponding  to  larger  antennas.  The  relationship 
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between  helix  circumference  and  spacing  for  m  =  1,2  is 
illustrated  in  Figure  3.3. 


FIGURE  3.3  Helix  Mode  Chart 


The  general  relationship  between  helical  radiation  mode, 
turn  length,  circumference,  and  spacing  is  provided  in  the 
following  equations  [Ref  2:  p.  289]: 

L^=C^+S^  (3.7) 
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—  =S,+/n  (3.6) 

P 

Once  the  relative  phase  velocity  has  been  determined  and  a 
transmission  mode  is  selected,  it  is  possible  to  resolve  the  far- 
field  radiation  pattern  of  the  helix.  As  long  as  the  helix  is 
long,  it  can  be  regarded  as  an  array  consisting  of  (n)  one  turn 
loops.  To  begin  radiation  pattern  computations  the  phase  shift 
(i|i)  of  each  equivalent  point  source  in  the  effective  array  factor 
of  the  helix  is  computed  as  follows: 


^=2n  (S^cosd- 


(radians) 


(3.8) 


In  (3.8),  (0)  corresponds  to  the  coaltitude,  or  deflection  angle 
from  the  axis  of  the  helix. 

As  a  result  of  the  symmetrical  nature  of  a  helical  antenna's 
main  lobe  the  following  relation  holds: 

iv/m)  (3.9) 

The  far-field  radiation  pattern  of  a  single  helical  turn  is 
reasonably  estimated  by  cos(6).  The  electric  field  pattern  (E) 
is  given  by  the  product  of  the  array  factor  and  the  individual 
turn's  pattern.  As  predicted  by  the  principle  of  pattern 
multiplication,  the  array  factor  corresponding  to  an  array  of 
isotropic  point  sources  dominates  the  field  pattern  generated  by 
a  single  turn  of  the  helix.  This  effect  can  be  seen  in  the 
following  formula  for  electric  field  [Ref  2:  pp.  294-295]: 

It  should  be  noted  that  unless  the  helix  is  very  short  (nS;> 
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(3.10) 


sin^^ 

£:=sin.:JL  ( - l-)cos0  (V/m) 

2^  sinl 
2 

<  .5),  ground  plane  reflections  and  their  effects  on  electric 
field  patterns  for  the  antenna  are  negligible.  Consequently, 
ground  parameters  are  not  required  for  this  application. 

The  radiation  intensity  (U)  at  any  far-field  observation 
point  is  a  function  of  and  per  the  following  equation  [Ref 
3 :  pp.  28-29] : 

J/=— ^  [  l£’el^+|£lp]  iw/ solid  ang)  (3.11) 

2t1o 

From  (3.9),  (3.11)  can  be  reduced  to: 

U=—\E\^  (w/solidang)  (3.12) 

Ho 

In  (3.11)  and  (3.12),  (rjo)  is  the  intrinsic  impedance  of  free 
space. 

The  average  radiated  power  (Prad)  antenna  is  given 

by: 

P^^^=jf^UdQ=j^''j^UsinBdBd^  <  v)  (3.13) 

In  (3.13),  (n)  is  a  sphere  in  the  far-field  surrounding  the 

antenna. 

It  is  impossible  to  determine  total  efficiency  (e^)  of  the 
helical  antenna  based  only  on  dimensional  information.  Therefore, 
an  antenna's  gain  (G)  cannot  be  precisely  determined  using  the 
following  general  gain  formula: 


10 


G=e^D^  (dimensionless) 


(3.14) 


However,  a  unique  feature  of  helical  antennas  is  that  input 
impedance  (ZJ  is  essentially  equal  to  input  resistance  (R)  when 
(3.1)  -  (3.3)  are  satisfied.  Fortuitously,  the  input  resistance 
of  the  helical  antenna  can  be  calculated  with  observed 
measurements  by  [Ref  2:  pp  277-278]: 

Axial  Feed:  R~140C^^  (Q)  (3.15) 

Peripheral  Feed:  {£i)  (3.16) 

If  antenna  feed  characteristic  impedance  (Z^)  is  known  or 
can  be  estimated,  then  reflection  efficiency  (e^)  can  be  computed 
from  the  voltage  reflection  coefficient  (F)  by  [Ref  4:  p.  460]: 

r=  ^  (dimensionless)  (3.17) 

ej.=l-|r|^  (dimensionless)  (3.18) 

Although  (3.1'/  and  (3.18)  provide  an  estimate  of  reflection 
efficiency,  no  other  helical  antenna  efficiency  terms  can  be 
determined  based  on  geometry  alone.  Thus,  all  other  components 
of  total  efficiency  are  assumed  to  be  unity  and  gain  is  expressed 
as: 

(dimensionless)  (3.19) 

EIRP  is  a  commonly  used  term  from  communications  that  is 
formally  defined  as  the  product  of  antenna  gai  i  and  total  power 


accepted  by  the  antenna  from  the  transmitter.  EIRP  is  determined 
as  follows  [Ref  5:  p.  62]: 


(W)  (3.20) 

A  functional  helical  antenna  will  exhibit  nearly  circular 
polarization  when  (3.1)  -  (3.3)  are  satisfied.  Axial  ratio  (AR) 
provides  a  figure  of  merit  for  circular  polarization  in  that  if 
it  is  equal  to  unity  the  polarization  of  the  antenna  is  exactly 
circular.  The  further  axial  ratio  is  from  one,  the  more 
elliptically  polarized  the  helical  wave  will  be  .  The  axial  ratio 
of  a  helix  is  [Ref  2:  pp.  301-307]: 

AR=  |L;^  (sin (a) -l/p)  I  (dimensionless)  (3.21) 

The  helical  antenna's  unit  polarization  vector  (aj  at  a 
given  point  in  the  far-field  is  computed  using  the  Cartesian 
components  of  electric  field.  The  Cartesian  components  of 
electric  field  and  the  antenna's  unit  polarization  vector  are 
determined  using  the  results  of  (3.9)  and  (3.10)  as  follows  [Ref 
5:  p.  555]: 


£^^=£■0003  (0)  CCS  V  i)'  -£^sin  (<J>) 

(V/m) 

(3.22) 

Ey^E^cos  (0)  sin(<J>)  +£^cos  (^) 

( V/m) 

(3.23) 

£^=-£0Sin(0)  (V/m) 

(3.24) 

a  ^(x,v,  z)  =  (dimensionless) 

(3.2C) 

y/\E(x,y,  z)  1^ 

In  (3.25),  (a^t  y  2)  are  the  Cartesian  unit  vectors. 


When  the  helical  antenna  is  used  for  reception  and  the 
incoming  wavers  electric  field  unit  vector  (p,)  at  a  given  point 
in  the  far-field  is  known  or  can  be  estimated,  the  polarization 
loss  factor  (PLF)  is  given  by  [Ref  3:  p.  51]: 

PLF=l5^d*p  (dimensionless)  (3.26) 

The  term  which  best  describes  an  antenna's  ability  to 
capture  incoming  electromagnetic  waves  and  extract  power  from 
them  is  maximum  effective  aperture  (A„)  .  Maximum  effective 
aperture  for  a  helical  antenna  is  [Ref  3:  p.  63]; 

A^^=[PLF]  (w^)  (3.27) 

In  (3.27),  (k)  is  wavelength  of  the  frequency  (f)  of  interest. 

Although  the  current  (I^,)  at  the  terminals  of  the  helical 
antenna  cannot  be  determined  from  dimensional  information  alone, 
if  the  current  is  assumed  to  be  unity  the  radiation  resistance 
(Rj)  and  maximum  effective  height  (hea)  are  estimated  by  [Ref  2: 
p.  42]: 

(Q)  (3.28) 


(m) 


(3.29) 


Helical  antenna  Mathcad  applications  are  valid  for  conductor 
diameters  given  by: 


.  005A.^di.05A.  (m) 


(3.30) 
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The  bandwidth  (BW)  of  an  operational  helix  is  determined  by 
the  high  and  low  frequencies  (fhigh*  flow)  corresponding  to  the 
dimensional  limits  of  (3.1).  Therefore,  bandwidth  can  be 
calculated  using  the  speed  of  light  (c)  by: 


f  _1.15c 

^  high  Q 

(Hz) 

(3.31) 

f  _  .8c 

low  ^ 

(Hz) 

(3.32) 

^high  ~  ^  low 

(Hz) 

(3.33) 

Most  of  the  parameters  calculated  by  the  Mathcad  helical 
antenna  applications  are  only  valid  if  the  observation  point  (r) 
is  in  the  far-field.  An  observation  point  is  considered  to  be  in 
the  far-field  if  all  of  the  following  are  satisfied  [Ref  3:  p. 

92]  : 

nl.eX  (/n)  (3.34) 

r^SnS  (m)  (3.35) 

(w)  (3.36) 

Table  3.1  and  Figure  3.4  compare  measured  data  to  that 
calculated  by  the  Mathcad  applications  for  a  10  turn  helical 
antenna  (D  =  .1074  meters  and  a  =  12.8°)  [Ref  6:  13-6  -  13-9]. 
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TABLE  3.1  Helical  Antenna  Data  Comparison 


ANTENNA 

MEASURED 

CALCULATED 

PARAMETER 

DATA 

DATA 

HALF-POWER 

BEAMWIDTH 

39° 

37  ° 

GAIN 

12.5  dB 

14.3  dB 

Electric  Field  Magnitude  in  V/rn 


Comparison  of  Hellool  Antenna  Electric  Fields 


IV.  THE  BEVERAGE  AMTENKA 


The  Beverage  antenna  is  a  single  wire  antenna  parallel  to 
the  ground  and  terminated  with  a  load  equal  to  the  characteristic 
impedance  (Z,,)  of  the  wire.  The  transmitter  or  receiver  of  a 
Beverage  antenna  has  one  end  connected  to  the  wire  and  the  other 
to  ground.  Because  of  its  matched  termination,  the  Beverage 
antenna  does  not  develop  a  significant  standing  voltage  wave 
along  its  length.  Therefore,  it  is  known  as  a  traveling  wave 
antenna.  The  relative  phase  velocity  (p)  of  the  wave  traveling 
down  the  antenna  is  typically  less  than  one.  Thus,  the  Beverage 
antenna  is  also  considered  a  slow  wave  antenna.  [Ref  3:  pp.372- 
374] 

Although  radiation  can  occur  at  any  non-uniformities  in  the 
device,  the  Beverage  antenna  primarily  generates  a  vertically 
polarized  cone  shaped  main  beam  that  points  in  the  direction  of 
the  traveling  wave.  The  geometry  of  a  Beverage  antenna  is 
illustrated  in  Figures  4.1  and  4.2. 
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FIGURE  4.1  The  Beverage  Antenna 
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Ground 


FIGURE  4.2  Elevation  Angle  (0)  of  a  Beverage 
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In  Figure  4.1,  (h)  is  the  antenna's  height  above  ground  and  (L) 

is  the  total  lencth  of  the  antenna.  In  Figure  4.2,  (6)  is  the 

angle  of  incidence  of  an  incoming  or  transmitted  wave  with 
respect  to  ground. 

Typically,  the  electrical  length  of  a  Beverage  antenna  (LJ 
will  be  on  the  order  of  0.5  to  2  wavelengths.  Maximum  length 

which  the  antenna  is  expected  to  operate  is  a  function 
of  both  arrival  angle  of  the  incoming  wave  and  relative  phase 
velocity.  A  precise  formula  for  maximum  length  is  [Ref  7:  p.  14]: 

= ; — -  [wavelengths) 

4(^-cos(e)) 

P 

Unfortunately,  use  of  (4.1)  is  normally  not  possible.  The 
wave's  angle  of  incidence  is  always  changing  and  can  only  be 
estimated  using  statistical  techniques.  In  addition,  relative 
phase  velocity  is  not  easily  determined  by  the  antenna's  geometry 
and,  consequently  is  not  generally  known.  Some  relative  phase 
velocity  measurements  have  been  conducted  over  the  following 
frequency  band  [Ref  7:  p.  19]: 

1.6  MHz  <,  f  ^  10.5  MHz  [Hz)  (4.2) 

If  the  frequency  of  interest  meets  the  criteria  of  (4.2),  then 
(p)  can  be  computed  by: 

p=. 65891  ( — ^_  )  038523821  ( diweDS  lonl  ess)  (4.3) 

1000 

If  inadequate  information  is  available  to  use  (4.1),  the 
Beverage  antenna  application  user  can  estimate  maximum  length 
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from  Table  4.1  [Ref  7;  p.  14]. 


TABLE  4.1  Maximum  Effective  Length  of  a  Beverage 

Antenna  in  Wavelengths 


6  in 

^Amax 

deg 

p=.89 

p=.91 

p=.93 

0 

2.02 

2.53 

3.32 

10 

1.80 

2.19 

20 

1.36 

1.57 

1.84 

30 

.97 

1.07 

1.19 

40 

.70 

.75 

.81 

50 

.  52 

.55 

.58 

Because  of  the  difficulty  determining  (6)  and  (p) ,  the 
Mathcad  application  assumes  that  the  difference  between 
frequencies  corresponding  to  0.5  to  2  wavelengths  is  the 
bandwidth  for  the  Beverage  antenna.  The  application  also 
computes  relative  phase  velocity  per  (4.3),  but  the  user  is 
cautioned  that  the  frequencies  of  interest  must  satisfy  (4.2). 

The  Beverage  antenna  transmits  or  receives  vertically 

polarized  waves.  In  the  case  of  reception,  the  question  might 

arise  as  to  how  a  wire  lying  parallel  to  the  ground  can  receive  a 

vertically  polarized  signal.  For  higher  frequency  operations 

that  utilize  sky  wave  propagation,  the  tilt  of  the  incoming  wave 

provides  a  horizontal  component  of  the  vertically  polarized 

electric  field  (E)  with  respect  to  the  ground  and  the  antenna. 

j  It  is  the  horizontal  component  of  the  wave  that  is  parallel  to 

! 

! 


the  antenna  which  generates  the  emf  on  the  wire. 

For  lower  frequencies  (i.e.,  <  300  Khz),  the  physics  of  a 
Beverage  antenna  is  much  more  complex.  In  lower  frequency 
applications,  the  ground  wave  is  the  principle  propagation  path. 
In  this  situation  there  is  negligible  tilt  to  the  wave  as  a 
result  of  propagation  path  geometry.  However,  as  the  vertically 
polarized  wave  travels  over  an  imperfect  conductor  the  electric 
field  closest  to  ground  begins  to  develop  a  forward  tilt  as 
pictured  in  Figure  4.3.  As  in  the  case  of  the  higher  frequency 
applications,  the  tilted  electric  field  of  the  low  frequency  wave 
has  a  horizontal  component  parallel  to  the  antenna  which  induces 
an  emf  on  the  wire. 

Originally,  use  of  the  Beverage  antenna  was  restricted  to 
low  frequencies  propagating  over  very  poor  ground.  The  Beverage 
antenna  is  now  often  used  in  an  attempt  to  reduce  noise 
interference  in  high  frequency  operations  over  excellent  ground. 
The  higher  frequency  skywave  propagation  path  provides  the 
necessary  tilt  to  receive  vertically  polarized  signals,  but  the 
near  perfect  ground  does  tilt  vertically  polarized  ground  waves 
from  nearby  noise  sources.  Thus,  the  high  frequency  Beverage 
antenna  becomes  a  highly  directive,  low  noise  device. 
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Ground 


FIGURE  4.3  E  Field  Over  an  Imperfect  Ground 

23 


Development  of  the  electric  field  pattern  of  a  Beverage 
antenna  begins  with  an  understanding  of  the  current  on  the 
antenna.  If  one  assumes  low  ohmic  losses,  matched  termination, 
and  negligible  attenuation  along  the  wire,  the  phasor  current 
amplitude  is  constant  and  the  phase  velocity  is  that  of  free 
space;  hence  [Ref  5:  p.  240]: 

J(z)=Joe-^*^  (A)  (4.4) 

In  (4.4)  the  antenna  is  assumed  to  lie  along  the  +z  axis,  (Iq)  is 
current  at  the  transmitter's  terminals,  and  (k)  is  the  free  space 
wavenumber  given  by: 

(;n'i)  (4.5) 

In  (4.5),  (1)  is  wavelength  of  the  frequency  (f)  of  interest. 
Beverage  antenna  Mathcad  applications  assume  (Iq)  is  normalized 
to  one  amp. 

With  the  current  defined  by  (4.4)  the  magnitude  of  the 
Beverage  antenna's  electric  field  is  obtained  by  [Ref  8:  pp.  315- 
316]  : 

l^l^30jcLJ„sln(e)|SinU)|  (4.e) 

r  JC 

In  (4.6),  (r)  is  the  distance  from  the  antenna  to  the  far-field 

observation  point,  L  is  the  length  of  the  antenna,  and  (X)  is 

given  by: 
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A’=-^  (1 -cos  (6) )  (radians) 


(4.7) 


The  electric  field  pattern  given  by  (4.6)  is  rotated  about 
the  +z  axis  to  form  the  three-dimensional  field  pattern  above  the 
ground  plane.  The  pattern  is  only  valid,  however,  in  the  far- 
field.  Therefore,  all  of  the  following  conditions  must  be 
satisfied  for  (4.6)  tn  apply  [Ref  3:  p.  92]: 

r^l.eX  (m)  (4.8) 

r^5L  (w)  (4.9) 


(m)  (4.10) 

Improvements  to  the  accuracy  of  (4.6)  can  be  made  for  far- 
field  radiation  patterns  if  one  accounts  for  the  effects  of  real 
ground.  Through  use  of  image  theory  and  the  fact  that  Beverage 
antennas  excite  vertically  polarized  waves,  the  electric  field 
pattern  equation  is  modified  as  follows  [Ref  5:  pp.  229-235]: 


30icLJ(jSin  (0) 


-j2khzos  ( -5  -6) 

'  ] 


(V/m)  (4.11) 


In  (4.11),  (Ty)  is  the  vertic:ai  reflection  coefficient  of 
ground.  Antenna  height  is  typically  less  than  one  wavelength. 
The  vertical  reflection  coefficient  is  given  by: 

In  (4.12),  (ej.  )  is  the  relative  complex  permittivity  of  the 
ground  under  the  antenna  and  is  calculated  as  follows: 

In  (4.13),  (ej.)  is  the  relative  permittivity  of  the  ground  and 
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e^/cos  ( -0)  -J  e^/-sin^  ( -0) 

r^= - = - ^  ^ -  {dimensionless)  (4.12) 

e^/cos  ( -j  -0)  e^z-sin^  ( -j  -0) 

e-/=e,-j— — %—  (dimensionless)  (4.13) 

'  ^  2uf€o 

(a)  is  the  conductivity  of  the  ground. 

The  direction  of  maximum  radiation  of  a  Beverage  antenna  may 
be  determined  from  (4.11).  However,  it  can  also  be  estimated 
quickly  by  the  following  empirical  formula  [Ref  5:  p.  241]: 

0max=cos-M3--^^)  (zadians)  (4.14) 

T 

with  the  magnitude  of  the  vertically  polarized  electric 
field  given  by  (4.11),  radiation  intensity  is  computed  as  follows 
[Ref  3:  pp.  28] : 

U=-^— [\e\^]  (w/solidang)  (4.15) 

2no 

In  (4.15),  (f7o)  is  the  intrinsic  impedance  of  free  space. 

The  radiated  power  and  directivity  of  a  Beverage  antenna  are 
determined  by  applying  radiation  intensity  to  standard  antenna 
formulas  as  follows  [Ref  3:  pp.  28-30]: 
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^rad=fV’'^sin(8)  ded4>  {W) 

J  0  J  0 


(4.16) 


47lt7(0  ) 

^o~ -  (dimensionless)  (4.17) 

^zad 

The  characteristic  impedance  of  a  Beverage  antenna  can 
be  estimated  by  its  dimensions  and  is  generally  resistive. 
Characteristic  impedance  of  a  Beverage  antenna  over  perfect 
ground  is  given  by  [Ref  7:  pp  19-21]: 

Zo=138log(44)  (Ci)  (4.18) 

a 

In  (4.18),  (d)  is  the  diameter  of  the  wire  in  the  same  units  as 

(h) .  Caution  must  be  exercised  when  using  this  value  of 
characteristic  impedance  since  any  sharp  transition  in  the  wire 
(i.e.,  vertical  downleads)  or  real  ground  effects  can  reduce  the 
accuracy  of  the  calculation.  Typical  values  for  a  Beverage 
antenna's  characteristic  impedance  are  200-300  ohms. 

If  the  impedance  (Zi)  of  the  terminating  load  of  a  Beverage 
antenna  is  known  or  can  be  estimated,  reflection  efficiency  (e^) 
can  be  determined  from  the  voltage  reflection  coefficient  (T)  as 
follows: 


r=  (dimensionless)  (4.19) 

Zj+Zq 

ej,=l-|rp  (dimensionless)  (4.20) 

Other  than  reflection  efficiency,  accurate  estimates  of 
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other  Beverage  antenna  losses  cannot  be  determined  by  geometry 
alone.  Nevertheless,  there  are  other  sources  of  lost  power. 
Since  a  Beverage  antenna  is  a  relatively  long  antenna  with  a 
matched  termination,  relatively  little  power  is  reflected  by  the 
load.  Instead,  most  of  the  power  supplied  by  the  transmitter 
that  is  not  radiated  is  absorbed  by  the  load  or  lost  as  heat  to 
the  ground. 

Gain  (G)  is  the  product  of  antenna's  directivity  and 
efficiencies.  The  Mathcad  Beverage  antenna  applications  express 
gain  as  [Ref  3:  p.43]: 

G=e^o  (dimensionless)  (4.21) 

Effective  isotropic  radiated  power  for  the  Beverage  antenna 
(EIRP)  is  the  product  of  the  power  radiated  by  the  antenna  and 
the  directivity.  (EIRP)  is  computed  by  [Ref  5:  p.  62]: 

EIRP=P,^^o  (4.22) 

Electromagnetic  waves  incident  upon  a  Beverage  antenna  are 

normally  assumed  to  be  vertically  polarized.  If  the  incoming 
wave  is  not  vertically  polarized,  a  polarization  mismatch 
occurs  with  the  antenna  and  losses  result.  Polarization  losses 
are  determined  from  a  polarization  loss  factor  (PLF)  given  as: 

PLF=\a^‘dl\^  (dimensionless)  (4.23) 

In  (4.23),  (aj  and  (a„)  are  the  unit  polarization  vectors  of  the 
antenna  and  wave,  respectively. 

Maximum  effective  aperture  (A^^)  is  estimated  from 
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directivity,  (PLF) ,  and  reflection  efficiency  as  follows  [Ref  3: 
pp.  51-63 ] : 

A^=[PLF]  im^)  (4.24) 

471 

The  maximum  effective  height  (heo)  of  a  Beverage  antenna  can 
be  determined  using  the  results  of  (4.24)  as  follows  [Ref  5;  p. 
42]  : 


(w) 


(4.25) 


In  (4.25),  (Rj.)  is  radiation  resistance  and  is  written  as: 


i?,= 


2P, 


zad 


(Q) 


(4.26) 


The  conductor  diameter  must  be  much  less  than  the  length  of 
the  Beverage  antenna  to  avoid  unwanted  radiation  from  vertical 
sections.  For  the  purpose  of  Beverage  antenna  Mathcad 
applications,  the  following  is  assumed  for  proper  antenna 
operations : 

Okl,  if  di.OlL  iin)  (4.27) 

Table  4.2  and  Figure  4.4  compare  measured  data  to  that 
calculated  by  the  Mathcad  applications  for  a  Beverage  antenna 
(L=110.4  meters  and  h=1.23  meters)  operating  over  dry  soil 
(a=.003  S/m  and  e^=12)  at  18  MHz.  Table  4.3  compares  measured 
and  calculated  data  for  a  Beverage  antenna  with  L=110.4  meters 
and  h=1.13  meters  operating  over  wet  soil  (a=.01  S/m  and  €,.=17) 


29 


at  5  MHz  [Ref  9:  pp.  22-26] 


TABLE  4.2  Beverage  Antenna  Data  Comparison 


ANTENNA 

MEASURED 

CALCULATED 

PARAMETER 

DATA 

DATA 

Omax 

18.5° 

19.2° 

TABLE  4.3  Beverage  Antenna  Data  Comparison 


ANTENNA 

MEASURED 

CALCULATED 

PARAMETER 

data 

DATA 

Zo 

450  n 

408  D 

P 

.93 

.91 

30 


Electric  Field  Magnitude  in  V/ m 
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V.  THE  LOOP  ANTENNA 


A  loop  antenna  is  a  coil  of  one  or  more  turns.  It  is 
commonly  used  as  a  receiving  antenna  for  operations  in  the  lower 
frequency  regions.  The  loop  antenna  is  also  used  for  direction 
finding  and  UHF  transmissions  [Ref  10:  p.  6-1].  Loop  antennas 
may  have  an  air  core  or  ferrite  core.  They  may  also  be 
electrically  large  or  small.  For  the  purpose  of  the  Mathcad 
applications,  a  loop  antenna  is  considered  electrically  small  if 
its  radius  (a)  satisfies  the  following  [Ref  3:  p.  181]: 

(w)  (5.1) 

6ji 

In  (5.1),  (A.)  is  wavelength.  The  geometry  of  both  large  and 

small  loop  antennas  is  illustrated  in  Figure  5.1. 
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The  radius  of  the  conductor  is  (b)  in  Figure  5.1.  For  all  loop 
antenna  Mathcad  applications  the  center  of  the  loop  is  the  origin 
and  the  antenna's  axis  is  aligned  parallel  to  the  +z  axis.  When 
Mathcad  applications  examine  the  performance  of  a  loop  over  a 
ground  plane,  the  coordinate  system  is  rotated  with  the  antenna 
as  necessary  to  obtain  the  desired  geometry  (i.e.,  the  axis  of  a 
vertical  loop  is  parallel  to  the  ground  and  the  +z  axis) . 

A.  THE  ELECTRICALLY  SMALL  LOOP 

Electrically  small  loops  are  normally  used  for  low  frequency 
reception  or  direction  finding.  Small  loops  are  poor 
transmitters  due  to  small  radiation  resistance  (Rj.)  and  low 
conduction-dielectric  efficiency  (Ccd)*  Transmitter  performance 
can  be  improved  with  increased  perimeter,  adding  additional 
turns,  or  insertion  of  a  ferrite  core  [Ref  3:  p.  164]. 

Two  key  assumptions  are  made  in  the  analysis  of  a  small 
loop.  First,  it  is  assumed  that  current  around  the  loop  is 
constant.  This  supposition  allows  the  loop  to  be  approximated  by 
an  infinitesimal  magnetic  dipole  centered  at  the  origin  and 
parallel  to  the  +z  axis.  Second,  it  is  presumed  that  the  various 
resistances  and  reactances  of  the  loop  can  be  computed  from 
dimensional  information  and  knowledge  of  the  antenna's  material 
properties . 

Given  the  above  assumptions  in  free  space,  the  electric 
field  (E)  for  a  small  loop  in  the  far-field  is  determined  by  [Ref 
3 :  pp.  168-169]  : 

In  (5.3),  (S)  is  the  cross-sectional  area  of  the  loop,  (f)  is  the 
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£,=£•0=0  (V/in) 


(5.2) 


^  _ic5fn„J„sin(e)e-^*^ 
= - 27 - 


(V/w) 


(5.3) 


frequency  of  interest,  (r)  is  the  distance  from  the  origin  to  the 
observation  point  in  the  far-field,  (Iq)  is  the  antenna  feed 
current,  (Mo)  is  the  permeability  of  free  space,  and  (k)  is  the 
free  space  wavenumber  given  by: 


(m-i)  (5.4) 

As  is  the  case  in  all  Mathcad  applications,  current  in  the 
loop  is  normalized  to  one  amp.  Since  electric  field  is  not  a 
function  of  (<^)  ,  the  field  pattern  is  symmetric  when  rotated 
about  the  antenna's  axis. 

It  should  be  noted  that  (5.2)  and  (5.3)  apply  to  all  small 
loops,  regardless  of  shape.  Thus,  (5.2)  and  (5.3)  can  be  used 
for  small  square  loops.  However,  it  should  also  be  noted  that 
loop  antenna  Mathcad  applications  assume  a  circular  loop  is  being 
analyzed  and  calculate  cross-sectional  area  based  on  the  radius 
provided  by  the  user.  In  order  to  use  the  applications  with 
rectangular  loops  or  loops  of  an  odd  shape,  the  user  must  compute 
an  equivalent  radius  (a)  that  will  yield  the  correct  area. 

An  observation  point  for  any  loop  antenna  is  assumed  to  be 
in  the  far-field  if  the  following  conditions  are  valid  [Ref  3: 
P-92]: 

ril.eX  im)  (5.5) 
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(w) 


(5.6) 


(/n)  (5.7) 

In  (5. 5) -(5. 7),  (D)  is  the  largest  dimension  of  the  loop.  The 

largest  dimension  of  the  loop  is  assumed  to  be  the  diameter. 

For  small  loops  in  free  space,  radiated  power  (Pr.d)  is 
estimated  as  follows: 


(ka)*\lj^  m  (5.8) 

In  (5.8),  (ry^)  is  the  intrinsic  impedance  of  free  space. 

The  directivity  (D^)  of  a  small  loop  in  free  space  is  1.5 
and,  ignoring  polarization  any  mismatches,  the  maximum  effective 
aperture  (Ag„)  of  a  lossless  loop  is  written  as  [Ref  3:  p.  175]: 

The  ohmic  resistance  of  any  loop  antenna  (Rohmic)  i  including 
multiple  turn  antennas,  is  estimated  by  the  following  [Ref  3: 
pp. 171-172] : 

(fl)  (5.10) 

In  (5.10),  (N)  is  the  number  of  turns,  (R^)  is  the  surface 
impedance  of  the  conductor,  (Rp)  is  the  ohmic  resistance  due  to 
proximity  effect,  and  (R^)  is  the  ohmic  skin  effect  resistance 
per  unit  length.  If  the  conductivity  (aj  of  the  conductor  is 
known,  the  surface  impedance  of  the  conductor  is  computed  by: 
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nfii. 


(Q) 


(5.11) 


Given  the  spacing  between  turns  (q)  has  been  measured,  the 
ratio  of  (Rp)  to  (R^)  is  estimated  using  Figure  5.2  [Ref  3:  p. 
172]  . 


FIGURE  5.2  Rp/R„ 

The  radiation  resistance  of  a  small  loop  in  free  space  is 
deteinmined  using  the  circumference  of  the  loop  (C)  as  follows 
[Ref  3:  pp.  170-171]; 


R^=20ti^{£)*N^  (Q) 


(5.12) 


When  the  radiation  and  ohmic  resistance  of  any  antenna 
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has  been  calculated,  the  conduction-dielectric  efficiency  of  the 
antenna  is  determined  by: 

^cd~~B - (dimemsionless)  (5.13) 

From  (5.10)  and  (5.12)  it  can  be  seen  that  ohmic  resistance 
is  directly  proportional  to  the  number  of  turns  while  radiation 
resistance  is  proportional  to  the  square  of  the  number  of  turns. 
Thus,  as  shown  in  (5.13)  conduction-dielectric  efficiency  can  be 
improved  by  increasing  the  number  of  turns  in  a  loop  antenna.  It 
can  also  be  seen  in  (5.10)  and  (5.12)  that  increasing  the  radius 
of  the  loop  improves  conduction-dielectric  efficiency. 

Additional  improvements  to  conduction-dielectric  efficiency 
may  be  made  by  inserting  a  ferrite  core  in  the  loop  antenna.  If 
a  core  is  added,  (5.12)  is  modified  as  follows: 

i?-  =  2  07^2  (^)«  (Jif)  2/^2  (Q)  (5.14) 


Effective  permeability  of  the  ferrite  core  in  (5.14)  is 

computed  by: 


^  ^ demag 


{H/m) 


(5.15) 


In  (5.15),  (Mf)  is  the  actual  permeability  of  the  core  material 
and  (Dde„ag)  is  an  experimentally  derived  demagnetization  factor. 
Demagnetization  factor  as  a  function  of  the  ratio  of  core  length 
to  diameter  is  shown  in  Figure  5.3  [Ref  3:  pp.  196-197]. 
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Cor#  L»ngth/Diom»t$r 

FIGURE  5.3  Demagnetization  Factor 

Any  antenna  parameter  that  requires  use  of  permeability  in 
its  formula  must  be  approximated  in  the  Mathcad  applications  by 
replacing  (Mj,)  with  (nj  [Ref  11:  pp.  86-89]. 

When,  a  loop  antenna  is  actually  employed,  it  is  not  in  free 
space  and  real  ground  must  be  considered.  Although  real  ground 
does  not  change  the  components  of  electric  field  given  by  (5.2), 
it  does  modify  (E^)  in  (5.3).  Since  the  orientation  of  the  loop 
with  respect  to  ground  determines  the  polarization  of  the  loop's 
electromagnetic  wave,  alignment  of  the  antenna  must  be  known 
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before  electric  field  can  be  correctly  computed.  In  order  to 
keep  the  loop  antenna  Mathcad  applications  moderately  simple, 
only  horizontal  and  vertical  loops  are  considered. 

Although  modeled  by  an  infinitesimal  vertical  magnetic 
dipole,  a  small  horizontal  loop  (i.e.,  +z  axis  perpendicular  to 
ground)  has  horizontally  polarized  electromagnetic  waves.  The 
horizontal  reflection  coefficient  (Th)  is  [Ref  5:  pp.  229-230]: 


cos  (0)  -^e^/-sin  (0)  ^ 
cos  (0)  +^e^i-sin{Q)  ^ 


(dimensionless) 


(5.16) 


In  (5.16),  the  relative  complex  permittivity  (e^-)  of  the  ground 
is  calculated  using  the  relative  permittivity  (e^)  and 
conductivity  (a)  of  the  ground  by: 


€./=€■- j  ^  (dimensionless)  (5.17) 

In  (5.17),  (€o)  is  the  relative  permittivity  of  free  space. 

The  total  electric  field  of  a  small  horizontal  loop  over 
ground  is  the  sum  of  the  direct  path  signal  and  ground  reflected 
signal.  One  can  use  image  theory  to  estimate  the  contribution  to 
the  far-field  pattern  by  ground  reflections.  The  total  electric 
field  expression  for  a  small  horizontal  loop  positioned  a 
distance  (h)  above  the  ground  is: 

^  ^  ^  e  [  1  -r^e  (6)  ]  (5.18) 

Note  in  (5.18)  that  the  image  antenna  is  a  vertical  infinitesimal 
magnetic  dipole.  Thus,  the  contribution  of  the  image  is 
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subtracted  from  the  contribution  of  the  actual  antenna.  The 
Mathcad  loop  antenna  applications  ignore  the  minor  contribution 
of  the  surface  wave  to  the  electric  field. 

The  total  electric  field  of  a  vertical  loop  over  real  ground 
can  be  determined  in  a  manner  similar  to  (5.18),  but  the  vertical 
reflection  coefficient  (r^)  must  be  used  for  the  image  antenna's 
contribution  to  the  far-field  pattern.  The  equation  for  vertical 
reflection  coefficient  is  [Ref  5:  pp.  231-232]: 


e.-cos  ( —  -6) 

r  _ 2 _ 

e^cos  (-1-0) 


1 


e^/-sin^(^-0) 


e.'-sin^  ( —  -0) 
'  2 


( dimens i onless) 


(5.19) 


For  a  vertically  mounted  small  loop,  the  image  antenna's 
contribution  is  added  to  that  of  the  actual  antenna.  Thus,  total 
electric  field  for  a  small  vertical  loop  is  written  as; 

^  ;c5f^^J^sin(0)  (5.2  0) 

▼v.re  2r 

It  should  be  noted  by  application  users  that  (5.20)  applies  to 
loops  located  in  a  coordinate  system  that  has  been  rotated  with 
the  axis  of  the  loop  such  that  the  +z  axis  is  parallel  to  the 
loop  axis  and  the  ground. 

If  the  electric  field  of  a  small  loop  antenna  over  real 
ground  is  known,  several  parameters  can  be  computed  using  general 
antenna  formulas.  Far-field  equations  for  radiation  intensity 
(U) ,  radiated  power,  directivity,  and  radiation  resistance  are  as 
follows; 
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(5.21) 


U=^\E\^  {W/ solid  ang) 

2no 

P^ad=||^ysin(0)  d0d4>  {W)  (5.22) 

Dq=47i-^^  (dimensionless)  (5.23) 

^zad 

=  (Q)  (5.24) 

In  (5.22),  (fl)  is  a  half  sphere  that  encloses  a  loop  over  ground 
in  the  far-field.  In  (5.23),  (U„„)  is  the  maximum  radiation 
intensity  anywhere  on  the  half  sphere  as  determined  by  applying 
(5.22)  . 

A  unique  feature  of  small  loop  antennas  is  that  most 
efficiency  terms  associated  with  it  can  be  calculated  from  the 
loop's  measurements.  One  of  the  reasons  for  this  attribute  is 
the  fact  that  input  reactance  (XJ  can  be  reasonably  estimated 
from  [Ref  5;  pp.  102-103]: 

-Y-27lfa^^[ln(8  J) -1.75]  (fl)  (5.25) 

Input  impedance  (ZJ  for  the  small  loop  is  found  using  the 
formula: 

Zi=R^^jXi  (Q)  (5.26) 

In  (5.26),  input  resistance  (RJ  is  the  sum  of  radiation  and 
ohmic  resistance  and  is  written  as: 


(5.27) 


Ri=Rr-Rotn.ic 

If  the  characteristic  impedance  of  the  antenna's  feed  line 
is  known  or  can  be  presumed,  the  voltage  reflection  coefficient 
(r)  and  reflection  efficiency  (€r)  are  determined  by: 

z  ~z 

r=— — -  (dimensionless)  (5.28) 

e^=l-irp  (dimensionless)  (5.29) 

As  previously  discussed,  the  polarization  of  a  small  loop 
matches  the  orientation  of  the  loop  (i.e.,  a  horizontal  loop  is 
horizontally  polarized) .  Thus,  if  the  polarization  of  an 
incoming  wave  is  known  or  can  be  estimated,  one  can  use  the  dot 
product  of  the  unit  polarization  vector  of  a  small  loop  antenna 
(a,)  and  the  unit  polarization  vector  of  an  incoming  wave  (a„)  to 
compute  the  polarization  loss  factor  (PLF)  as  follows  [Ref  3;  p. 
51]  : 

PLF=|d^-da|^  (dimensionless)  (5.30) 

with  the  efficiency  and  loss  terms  of  the  small  loop 
estimated  by  (5.13),  (5.29),  and  (5.30);  gain  (G) ,  maximum 
effective  aperture  (Ae„)  ,  and  effective  isotropic  radiated  power 
(EIRP)  can  be  expressed  as  [Ref  3  pp.  43-63]: 
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{dimensionless) 


(5.31) 


(^)  (/n^) 


(5.32) 


(W)  (5.33) 

Maximum  effective  height  (h^)  can  be  determined  from  the 
maximum  effective  aperture  and  radiation  resistance  of  the 
antenna  as  follows  [Ref  2:  p.  42]: 


(5.34) 


Mathcad  applications  are  developed  assuming  bandwidth  (BW) 
of  an  antenna  is  the  range  of  frequencies  over  which  all 
computations  will  be  valid.  For  the  small  loop  antenna  Mathcad 
applications  to  be  valid,  (5.1)  must  hold.  Thus,  the  radius  of 
the  loop  will  define  lowest  operating  frequency  (fnin)^ 
highest  operating  frequency  (f„«x),  and  the  bandwidth  of  the 
antenna  as  shown  below: 


4in2^0 


(5.35) 


6na 


(5.36) 


(HZ)  (5.37). 

The  electric  field  patterns  and  antenna  parameters  obtained 
using  the  small  loop  Mathcad  applications  for  a  small  loop  in 
free  space  are  identical  to  published  results  obtained  using 
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method  of  moments  techniques  [Ref  3:  pp.  169-180].  A  comparison 
of  the  results  of  the  Mathcad  applications  with  measured  data  for 
a  small  loop  (a=.25  meters,  b=.005  meters)  located  2.5  meters 
above  a  reflecting  plane  (a=6xl0^  S/m,  £r=l)  receiving  a  30  MHz 
signal  is  provided  in  Table  5.1  [Ref  6:  p.  5-14]. 


TABLE  5.1  Small  Loop  Antenna  Data  Comparison 


ANTENNA 

MEASURED 

CALCULATED 

PARAMETER 

DATA 

DATA 

DIRECTIVITY 

7.0  dB 

6.5  dB 

B.  THE  ELECTRICALLY  LARGE  LOOP 

Electrically  large  loops  are  those  loops  that  do  not  satisfy 
(5.1).  Use  of  these  antennas  is  somewhat  rare,  with  radii 
exceeding  one  wavelength  normally  not  practical.  Large  loops 
are,  in  general,  significantly  more  difficult  to  analyze  than 
small  loops  since  large  loops  cannot  be  approximated  by 
infinitesimal  magnetic  dipoles.  In  addition,  the  input  reactance 
of  a  large  loop  antenna  cannot  be  calculated  from  the  loop's 
geometry. 

The  polarization  of  a  large  horizontal  loop  over  real  ground 
is  horizontal,  but  the  polarization  of  a  large  vertical  loop  is 
not  vertical.  Due  to  the  complexities  of  computing  the 
polarization  and  ground  reflection  coefficients  of  a  large 


vertical  loop,  the  Mathcad  large  loop  antenna  applications  will 
only  examine  the  free  space  and  horizontal  cases. 


To  begin  analysis  of  the  far-field  radiation  pattern  of  a 
large  loop  in  free  space,  one  roust  assume  that  the  current  on  the 
loop  is  constant.  As  illustrated  in  Figure  5.4,  this  is  an 
approximation  that  deteriorates  with  the  size  of  the  loop  [Ref  3: 
p.  184]. 


FIGURE  5.4  Current  Magnitude  Distribution  for  a  Large  Loop 


If  one  assumes  that  loop  current  is  constant,  (5.2)  applies 
to  large  loop  antennas.  The  other  component  of  electric  field 
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(E^)  is  approximated  by  [Ref  3:  pp.  176-178]: 


2r 


ikasin  (0) ) 


(5.38) 


In  (5.38),  Ji(kasin(0))  is  a  Bessel  function  of  the  first  kind  of 
order  one.  With  the  horizontal  reflection  coefficient  of  the 
ground  computed  per  (5.16),  total  electric  field  of  the  large 
horizontal  loop  over  real  ground  is  given  by: 


gftor  =  ^ ^ ^2 r  ^ ikasinid) 


(5.39) 


Radiated  power  for  a  large  loop  in  free  space  is  computed 
using  (5.22),  where  radiation  intensity  is  given  by  (5.21)  and 
electric  field  is  computed  by  (5.38).  It  should  be  noted  that 
the  far-field  conditions  of  (5.5) -(5.7)  must  be  satisfied  if  the 
Mathcad  large  loop  antenna  electric  field  and  radiated  power 
calculations  are  to  be  valid. 

Large  loop,  free  space  approximations  for  radiation 
resistance,  directivity,  maximum  effective  aperture,  and  maximum 
radiation  intensity  are  as  follows  [Ref  3:  p.  181]  [Ref  11:  pp. 
78-79]  : 


47 


(5.40) 


R^=60n^  (Q) 

At 

D^=.6  82(-y)  {dimensionless)  (5.41) 

Ag^=.0543  (XC)  im^)  (5.42) 

U^=  (.584)^  {w/ solid  ang)  (5.43) 

For  the  horizontal  loop  over  real  ground,  radiation 
intensity  is  given  by  (5.21),  where  the  electric  field  is 
computed  in  (5.39).  With  radiation  intensity  known,  radiated 
power  is  calculated  by  (5.22)  and  effective  isotropic  radiated 
power  is  computed  using  (5.33).  Radiation  resistance  and 
directivity  for  the  horizontal  loop  are  determined  using  generic 
antenna  formulas  as  follows: 

(5-44) 

\^o\ 

Dq-  idimensionlesss)  (5.45) 

^zad 

Ohmic  resistance  for  the  large  loop  antenna  is  given  by 
(5.8).  With  radiation  resistance  given  by  (5.40)  or  (5.44),  as 
applicable,  the  conduction-dielectric  efficiency  of  the  antenna 
is  computed  using  (5.13). 

Without  knowledge  of  the  input  reactance,  it  is  impossible 
to  compute  reflection  efficiency  of  a  large  loop  antenna. 
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Therefore,  Mathcad  applications  for  large  loops  assume  a  factor 
of  unity  for  reflection  efficiency  and  calculate  gain  as  follows: 

(dimensionless)  (5.46) 

In  general,  the  polarization  loss  factor  given  in  (5.30) 
cannot  be  determined  for  an  arbitrary  large  loop.  Thus,  (PLF) 
for  large  loops  is  assumed  to  be  unity  and  maximum  effective 
aperture  is  given  by: 

('""I 

With  maximum  effective  aperture  determined  for  a  large  loop 
using  (5.47),  maximum  effective  height  is  found  by  applying 
(5.34)  . 

The  bandwidth  of  a  large  loop  is  determined  by  the  loop's 
radius  as  follows: 

-^sasA.  (m)  (5.48) 

671 

— (Hz)  (5.49) 

6a7i  a 

SV--(1-— )  (Hz)  (5.50) 

a  671 

Figure  5.5  compares  the  electric  field  pattern  obtained  from 
the  Mathcad  applications  for  a  large  loop  (a=.46  meters,  b=.05. 
meters)  in  free  space  receiving  a  330  MHz  signal  with  published 
results  obtained  from  method  of  moments  techniques  [Ref  3:  p. 
180].  Table  5.2  compares  antenna  parameters  calculated  using  the 
Mathcad  applications  with  measured  results  for  a  large  loop 
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(a=.46  meters,  b=.05  meters)  located  .72  meters  above  a 
reflecting  ground  plane  (a=6xl0^  S/m,  er=l)  receiving  a  104  MHz 
signal  [Ref  6:  p.  5-14]. 


TABLE  5.2  Large  Loop  Antenna  Data  Comparison 


ANTENNA 

MEASURED 

CALCULATED 

PARAMETER 

DATA 

DATA 

INPUT  RESISTANCE 

150  n 

207  n 

DIRECTIVITY 

7.0  dB 

6.7  dB 
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The  vertical  directivity  (Dq)  of  a  horizontal  dipole  over 
real  ground  may  be  improved  by  placing  identical  elements  in  a 
straight  line  above  the  original.  A  vertical  line  of  horizontal 
dipoles  is  commonly  referred  to  as  a  bay.  In  a  similar  manner, 
one  may  place  additional  bays  adjacent  to  the  original  and 
achieve  an  improvement  in  horizontal  directivity.  Further 
improvements  to  the  gain  (G)  of  this  array  of  bays  may  be 
realized  by  placing  a  reflector  on  one  side  to  simulate  the 
existence  of  a  conducting  plane.  The  entire  arrangement  of  bays 
and  reflector  is  referred  to  as  either  a  bedspring  or  a  curtain 
antenna.  A  bedspring  array  is  frequently  used  for  high  frequency 
(3-30  MHz)  short  wave  radio  systems  [Ref  10:  pp.  21-1  -  21-6].  A 
typical  bedspring  antenna  arrangement  is  illustrated  in  figure 
6.1. 
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To  implement  the  Mathcad  software  for  the  bedspring  antenna 
applications,  the  following  dimensions  are  needed: 

Zj  =  height  of  the  i*^^'  element  above  ground 

-  Z^-i  =  vertical  spacing  between  the  i'^^’  and  the  (i-1)^^ 
element 

Yj  =  horizontal  position  (dipole  center)  of  the 
i^^  bay  (Yq  =  0) 

Y,  -  Yi-i  =  horizontal  spacing  between  the  i*"^  and  the  (i-l)‘^‘’ 
bay 

N  =  number  of  bays 

M  =  number  of  elements  in  each  bay 

Xi  =  reflector  position 

1  =  half-length  of  each  element 

The  geometry  for  the  bedspring  antenna  is  shown  in  Figure  6.2. 
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In  Figure  6.2  ,  (i|i)  is  the  angle  between  the  +y  axis  and  the 
vector  (r)  from  the  origin  to  the  observation  point  in  the  far- 
field.  Bedspring  antenna  Mathcad  applications  assume  that  all 
antenna  elements  are  identical  horizontal  dipoles  parallel  to 
each  other  and  the  y  axis.  The  applications  also  presume  that 
the  reflector  is  a  perfect  vertical  conducting  plane,  parallel  to 
the  y-z  plane,  and  located  in  the  -x  half-space. 

The  assumption  that  the  reflector  is  a  perfect  conducting 
plane  is  fairly  reasonable,  as  experimental  data  has  shown  that 
the  reflector  typically  improves  gain  by  2.5  to  3.0  dB.  The 
reflector  screen  may  be  constructed  from  tuned  elements,  such  as 
half-wave  dipoles,  or  it  may  consist  of  a  pattern  of  closely 
spaced,  parallel  wires.  The  reflector  is  normally  located  about 
one-guarter  wavelength  behind  the  antenna  [Ref  10:  p.  21-6]. 

The  electric  field  (E)  pattern  for  the  bedspring  antenna  is 
computed  using  the  principle  of  pattern  multiplication  for 
antenna  arrays.  Using  the  identity  cos (i|i)  =sin(6) cos (0)  ,  one  may 
predict  the  electric  field  components  of  a  single  bay.  If  all 
elements  in  the  bay  are  excited  by  a  sinusoidal  current  with 
maximum  amplitude  (I„)  ,  the  electric  field  components  are  [Ref 
12:  pp.  229-231]: 


E,,  =  -j60I,±^  fe°s(J^isin„(e)5in(4.))-cos(;:i)]  (3) 

”  r  [l-sinMe)sinM(|>n 

{V/m) 


(6.1) 
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^  g  (C03  (0)  -cos  (63)  )  g -jZJcZjCOS  (0)  j 


(6.5) 


( diinensi  onl  ess) 


The  relative  amplitude  of  excitation  (CJ  of  the  i'^^  dipole  with 
respect  to  the  first  element  in  the  first  bay  is  required  for 
both  (6.4)  and  (6.5).  Equations  (6.4)  and  (6.5)  also  use  (Bg)  to 
represent  the  desired  vertical  scan  angle  of  the  antenna.  The 
vertical  scan  angle  is  approximately  equal  to  the  progressive 
phase  shift  from  one  element  to  the  next  in  the  bay.  In  (6.4), 
(R,)  is  the  vertical  reflection  coefficient  over  real  ground.  In 
(6.5),  (Rh)  is  the  horizontal  reflection  coefficient  over  real 
ground.  These  reflection  coefficients  are  determined  by  [Ref  5: 
pp.  229-235]: 
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( dimensionless) 


(6.6) 


_  e^/cos  (0)  -y'e^z-sin^  (0) 
e^/cos  (0)  -t-^e^z-sin^  (0) 


cos  (0)  -^e^i-sin(Q)  ^ 
cos  (0)  +^ej<-sin(0)  ^ 


( dimens i onl ess) 


(6.7) 


The  relative  complex  permittivity  (Cr)  the  ground  needed 
to  determine  both  reflection  coefficients  is  calculated  using  the 
relative  permittivity  (e^)  of  the  ground,  the  conductivity  of  the 
ground  (a),  and  the  permittivity  of  free  space  (Co)  as  follows: 


€,z=e,-j — (dimensionless)  (6.8) 

^  ^  2Tife^ 

The  electric  field  components  of  (6.1)  and  (6.2)  apply  only 
in  the  far-field.  Thus,  all  of  the  following  conditions  must 
hold  if  the  computed  electric  fields  are  to  be  valid  [Ref  3:  pp. 
92-93] : 

r>\.6X  (m)  (6.9) 


T>5D  (m) 


(6.10) 


r>i^  (m)  (6.11) 

A* 

In  (6.10)  and  (6.11),  (D)  is  the  largest  physical  dimension  in 

any  direction  of  the  antenna  and  is  equal  to  the  bedspring's 
diagonal  length. 

If  one  assumes  that  the  amplitude  and  phase  of  the  feed 
current  in  corresponding  elements  in  each  bay  is  the  same,  the 
array  factor  (Sy)  for  N  bays  is  written  as  [Ref  12:  pp.  229-231]: 
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( dimensi onl ess) 


(6.12) 


=e:. 


In  (6.12),  (0o)  is  the  azimuthal  scan  angle  of  the  antenna.  The 
azimuthal  scan  angle  is  approximately  equal  to  the  progressive 
phase  shift  between  bays  in  the  bedspring  antenna. 

The  array  factor  accounting  for  the  perfect  image  from  the 
reflector  (S^)  is: 


(dimensionless)  (6-13) 


The  total  electric  field  components  (Egt,  E^^)  are  computed 
by  taking  the  product  of  all  appropriate  array  factors  and  the 
pattern  for  a  single  bay  as  follows: 

(V/m)  (6.14) 

(V/lD)  (6.15) 

The  radiation  intensity  (U)  of  the  bedspring  antenna  is 

determined  from  (Ej^)  and  (E^t)  by  [Ref  3:  pp.  27-29]: 

J7=-^  [  |£'9j.|^+|£'^trPl  (W/ solid  ang)  (6.16) 

2t1o 

In  (6.16),  (r?o)  is  the  intrinsic  impedance  of  free  space. 

Using  the  radiation  intensity  calculated  in  (6.16),  one  may 
calculate  radiated  power  (Prad)  by  computing  the  following 
integral  over  the  quarter  sphere  (n)  through  which 
electromagnetic  energy  from  the  antenna  flows: 

j^Usin(d)  ddd^  (W)  (6.17) 


Mathcad  bedspring  antenna  applications  assume  that  each 
element  is  excited  by  a  feed  current  with  the  same  maximum 
amplitude.  Although  this  assumption  may  seem  restrictive,  it  has 
been  experimentally  determined  that  maximum  gain  is  obtained  from 
a  bedspring  antenna  if  all  radiator  currents  are  of  equal 
amplitude.  Therefore,  the  following  approximation  for  antenna 
feed  current  (Iq)  holds  [Ref  10:  p.  21-17]: 

(amps)  (6.18) 

It  should  be  noted  that  the  bedspring  antenna  applications 
normalize  all  element  excitation  currents  to  one  amp. 

Given  the  magnitude  of  the  feed  current  for  the  entire 
antenna,  radiation  resistance  (Rj.)  for  the  complete  assembly  is 
determined  as  follows: 

2  p 

(Q)  (6.19) 

I  ■^ol 

Directivity  is  computed  from  radiation  intensity  and 
radiated  power  by: 

£>o=  (dimensionless)  (6.20) 

^zad 

In  (6.20),  (U„„)  is  the  maximum  value  of  radiation  intensity 

anywhere  on  the  quarter  sphere  encompassing  the  antenna's 
emissions . 

Gain  is  the  product  of  total  antenna  efficiency  (€t)  and 
directivity.  Gain  is  normally  expressed  as: 
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(dimensionless) 


(6.21) 


Unfortunately,  the  total  efficiency  of  a  bedspring  antenna 
cannot  be  easily  determined  from  its  dimensions.  For  example, 
input  impedance  is  a  key  parameter  needed  to  calculate  an 
antenna's  reflection  efficiency.  However,  input  impedance  has 
several  components  that  are  extremely  difficult  to  determine. 

Self  impedance  of  each  element,  mutual  impedance  between  real 
elements  in  the  array,  mutual  impedance  between  real  elements  and 
image  elements  from  the  reflector,  and  mutual  impedance  between 
real  elements  and  imperfect  images  in  the  ground  plane  all 
contribute  to  input  impedance.  Because  of  the  complexity  of 
computations,  input  impedance  of  the  bedspring  antenna  can  only 
be  estimated  by  an  extensive  method  of  moments  algorithm. 

Other  problems  computing  antenna  efficiency  include  the 
possible  existence  of  tuning  devices  in  the  feed  lines;  probable 
lack  of  sinusoidal  current  distributions  for  elements  of 
arbitrary  length;  and  inability  to  properly  model  the  ground, 
reflector  screen,  and  characteristic  impedance  of  the  feed 
assembly.  Nevertheless,  bedspring  antenna  gain  can  be  adequately 
estimated  despite  the  extensive  calculations  associated  with 
precise  modeling.  Experience  has  shown  that  a  properly  tuned 
bedspring  antenna  operating  within  the  band  of  expected 
frequencies  exhibits  very  little  loss.  Consequently,  gain  is 
estimated  in  Mathcad  bedspring  applications  as  [Ref  10:  p.  21-3]: 

G=10log(l?„) -2  (dB)  (6.22) 
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Effective  Isotropic  Radiated  Power  (EIRP)  is  the  product  of 
radiated  power  and  directivity  and  is  computed  by  [Ref  5:  p.  62]: 

EIRP=P^^J)^  iW)  (6.23) 

The  bedspring  antenna's  unit  polarization  vector  (a,)  may  be 
found  by  converting  the  electric  field  components  of  (6.14)  and 
(6.15)  to  Cartesian  coordinates  as  follows  [Ref  4:  pp.  35-36, 
364-367] : 


E^=EqCOS  (0)  cos  ((|))  -E^sin  (<t)) 

(V/m) 

(6.24) 

Ey=EQCos  (6)  sin ((^)  +E^cos  (<t>) 

(V/m) 

(6.25) 

E^^-E^sinid)  (V/m) 

(6.26) 

a ^{x,  V.  z)  =  (dimensionless) 

(6.27) 

^\E{x,y,  z)  )2 

In  (6.27),  (a^ y j)  are  the  unit  vectors  for  the  Cartesian 
coordinate  system. 

The  polarization  loss  factor  (PLF)  of  a  bedspring  antenna  at 
a  point  in  the  far-field  for  a  given  incoming  wave  with  unit 
polarization  vector  (a„)  is  expressed  as  [Ref  3:  p.  51]; 

PLF=\a^,'dl\^  {dimensionless)  (6.28) 

Without  knowledge  of  antenna  efficiencies,  one  cannot 
exactly  predict  an  antenna's  maximum  effective  aperture  (Ae^)  or 
maximum  effective  height  (h^n,)  .  Nonetheless,  one  may  assume  a 
lossless  antenna  system  and  approximate  these  parameters  as 
follows  [Ref  2:  pp.  29-43]: 
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A^^DJPLF)  (^)  (/n2) 


(6.29) 


(6.30) 


For  a  half-wave  dipole  assejibly,  such  as  that  pictured  in 
Figure  6.1,  a  complex  current  that  achieves  maximum  gain  is  fed 
to  each  element  only  within  a  few  percent  of  the  frequency  (fx/z) 
whose  wavelength  matches  the  length  of  the  half -wavelength 
dipoles.  Therefore,  the  bandwidth  (BW)  of  a  Figure  f.l  type 
bedspring  antenna  is  [Ref  10:  p.  21-16]: 


^high~^  •  ® 2^1/2 

(Hz) 

(6.31) 

^low~  • 

(Hz) 

(6.32) 

^ bigh~ ^ low 

(Hz) 

(6.33) 

In  (6.31),  (fh,gj  is  the  upper  frequency  of  the  antenna.  In 
(6.32),  (flow)  is  the  lower  bound  on  operating  frequency. 

If  wideband  operations  are  required  for  a  bedspring  antenna, 
a  symmetrical  feed  arrangement  as  shown  Figure  6.3  may  be 
employed. 
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The  bandwidth  for  a  symmetrical 

given  by  [Ref  10:  p.  21-16]: 

feed  bedspring 

antenna  is 

{Hz) 

(6.34) 

^lou~  • 

(Hz) 

(6.35) 

{Hz) 

(6.36) 

Figure  6.4  compares  the  electric  field  pattern  computed  by 
the  bedspring  antenna  Mathcad  applications  with  measured  results 
for  a  two-bay  (Yj^=26  meters)  ,  four-stack  (Zi=13  meters)  bedspring 
antenna  with  reflector  (Xi=7  meters)  operating  at  10  MHz  over 
soil  (a=.01  S/m,  6^=10)  [Ref  12:  p.  115]. 
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Electric  Field  Magnitude  in  V/ 


Comparison  of  Bedspring  Antenna  Electric  Fields 


FIGURE  6.4  Bedspring  Antenna  Electric  Field  Patterns 
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VI 


Spiral  antennas  are  a  family  of  two-  and  three-dimensional 
structures  that  maintain  a  constant  input  impedance,  beam 
pattern,  gain,  and  polarization  as  well  as  many  other  parameters 
over  a  wide  range  of  frequencies-  Spiral  antennas  are  commonly 
referred  to  as  frequency  independent,  or  broadband,  devices.  A 
two-dimensional  spiral  is  called  a  planar  spiral,  while  a  three- 
dimensional  spiral  is  usually  termed  a  conical  spiral.  Planar 
and  conical  spiral  antennas  are  commonly  used  in  applications 
such  as  direction  finding,  missile  guidance,  and  satellite 
tracking.  [Ref  6:  pp.  14-2  -  14-3] 

Although  there  are  several  types  of  planar  and  conical 
spiral  antennas,  Mathcad  applications  will  fully  analyze  only 
those  antennas  with  reasonably  simple,  closed  form  equations: 
equiangular  planar  spirals  and  conical  log-spirals.  Mathcad 
applications  assume  that  the  base  of  all  spirals  lies  in  the  x-y 
plane  and  is  centered  at  the  origin,  that  the  axis  of  all  spirals 
is  parallel  to  the  z  axis,  and  that  the  spirals  are  in  free 
space. 

A.  THE  PLANAR  SPIRAL  ANTENNAS 

There  are  three  major  categories  of  planar  spiral  antennas: 
the  equiangular  spiral,  the  Archimedean  spiral,  and  the  log- 
periodic  spiral.  The  geometry  of  all  planar  spirals  is  pictured 
in  Figure  7.1. 
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In  Figure  7.1,  the  spiral  angle  (6)  is  the  angle  between  any 
radial  line  from  the  origin  and  a  tangent  to  any  edge  of  the 
spiral,  (r)  is  the  distance  to  any  point  on  the  spiral  from  the 
origin,  and  (r^)  is  the  distance  from  the  origin  to  the  spiral's 
feed  point.  The  Mathcad  spiral  antenna  application  user  should 
not  confuse  the  radial  distance  from  the  origin  to  any  point  on 
the  spiral  (r)  and  the  distance  from  the  origin  to  an  observation 
point  in  the  far-field  (r^f)  . 

Spiral  antennas  may  be  constructed  from  wires  or  sheets  of 
metal.  For  low  power,  receive  only  operations,  spirals  may  also 
be  built  using  printed  circuit  technology.  A  more  rugged,  all¬ 
purpose  antenna  is  constructed  by  simply  cutting  the  spiral  edges 
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from  a  sheet  of  metal  and  running  coaxial  feed  lines  along  the 
spiral  arms.  A  dummy  feed  line  may  also  be  run  on  an  opposing 
arm  for  symmetry  [Ref  6:  pp.  14-4  -  14-7]. 

The  physical  dimensions  of  the  spiral  arms  determine  the 
type  of  spiral  antenna  and  the  antenna's  parameters.  An 
equiangular  spiral  is  one  whose  edges  or  wires  satisfy  the 
following  [Ref  5:  p.  283]: 

r=r^e-^  (m)  (7.1) 

In  (7.1),  (a)  is  an  arbitrary  constant  called  the  flare  rate.  If 

the  flare  rate  is  a  negative  number,  the  spiral  is  considered 
left-handed.  If  the  flare  rate  is  positive,  the  spiral  is  right- 
handed. 

When  sheet  metal  is  used  to  construct  an  equiangular  spiral, 
(7.1)  defines  the  coordinates  of  one  edge  of  one  spiral  arm.  The 
next  edge  (Vz)  is  cut  using  the  same  spiral  curve  as  (7.1),  but 
with  an  angular  arm  width  (5)  as  follows: 

(in)  (7.2) 

Spiral  antennas  are  usually  symmetrical.  Thus,  for  a  two  arm 
spiral,  edges  (rj)  and  (rj  are  given  by: 

(m)  (7.3) 

r,  =  r^e^  (w) 

Normally,  flare  rate  is  converted  to  a  factor  called 
expansion  ratio  (e„)  which  is  written  as  [Ref  5:  p.284]: 
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^  ^  ( dimensi  onless)  (7.5) 

r(<|)) 

A  typical  value  for  the  expansion  ratio  is  4. 

The  Archimedean  spiral  has  many  of  the  same  characteristics 
as  the  equiangular  spiral,  except  that  any  point  on  the  edge  of 
an  Archimedean  spiral  is  written  as: 

r=r^(J>  im)  (7.6) 

Likewise,  the  performance  of  a  log-periodic  spiral  is 
similar  to  an  equiangular  spiral  except  that  its  edges  are 
defined  by: 

r=r^a*  (m)  (7.7) 

For  a  log-periodic  spiral  it  can  be  shown  that  the  following 
formula  is  always  satisfied  [Ref  2:  pp.  697-698]: 

<l)  =  tan (p)  In (r)  {radians)  (7.8) 

In  (7.8),  the  spiral  angle  (fi)  is  the  same  as  that  illustrated  in 
Figure  7.1  and  is  constant  at  any  point  on  the  log-periodic 
structure.  Because  of  their  broadband  characteristics  and  ease 
of  construction,  log-periodic  spirals  are  regularly  used  in  the 
lower  millimeter  wave  region  [Ref  12:  p.  17-28]. 

The  lowest  operating  frequency  (flow)  /  highest  operating 
frequency  (fhigh)/  arid  bandwidth  (BW)  of  spirals  are  functions  of 
the  antenna's  dimensions  and  the  feed  arrangement.  For 
equiangular  and  Archimedean  devices,  the  minimum  radius  (r^)  and 
wavelength  of  the  highest  operating  frequency  are  approximately 
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correlated  as  follows  [Ref  5:  pp.  284-285]: 


(m)  (7.9) 

The  wavelength  (l^ow)  corresponding  to  the  lowest  operating 
frequency  is  set  by  the  overall  radius  (R)  of  the  structure  as 
follows : 


R= 


^low 

4 


(m) 


(7.10) 


With  the  antenna's  upper  and  lower  frequency  limits 
established,  bandwidth  may  be  expressed  as; 


(fiz)  (7.11) 

The  frequency  limits  and  bandwidth  of  a  log-periodic  spiral 
are  also  determined  by  its  dimensions.  Bandwidth  of  a  log 
periodic  spiral  is  given  by  (7.11),  where  the  lower  frequency 
limit  corresponds  to  the  wavelength  computed  using  (7.10)  and  the 
upper  frequency's  wavelength  is  calculated  as  follows  [Ref  2:  p. 


700] : 


_  ^high 
20 


(/n) 


(7.12) 


A  spiral  antenna  is  frequency  independent  in  that  most 
antenna  parameters  do  not  vary  over  the  bandwidth  of  the  antenna, 
which  can  be  a  considerable  range  of  frequencies.  Input 
impedance  is  one  of  these  parameters  and  may  be  computed  using 
the  concept  of  complimentary  antennas.  The  complement  of  a 
spiral  is  formed  by  replacing  metal  with  air  and  air  with  metal. 
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The  impedances  of  the  spiral  and  its  complement  are  (Z„gt.ai)  and 
(Z.ir)  ,  respectively.  These  two  impedances  are  real,  frequency 
independent,  and  for  a  two-arm  spiral  related  as  follows  [Ref  5: 
p.  283]: 

7  7  -ilf  (7.13) 

In  (7.13),  (Ho)  is  the  intrinsic  impedance  of  free  space.  If  the 
antenna  and  its  complement  are  identical,  the  antenna  is  called 
self-complimentary  and  the  angular  arm  width  in  (7.2)  is  7r/2.  In 
the  specific  case  of  a  two-arm,  self-complimentary  spiral,  the 
impedance  of  the  antenna  is: 

Self-complimentary  spirals  are  fairly  common  because  they 
yield  desirable  radiation  patterns.  Measured  input  impedance  is 
typically  between  120-160  ohms  for  these  devices,  which  is  lower 
than  the  theoretical  results  of  (7.14).  The  presence  of  a  feed 
structure,  the  finite  size  of  the  antenna,  and  the  finite 
thickness  associated  with  the  spiral's  arms  are  the  reasons 
measured  input  impedance  is  lower  than  its  theoretical  value  [Ref 
5:  pp.  285-286] . 

Mathcad  spiral  antenna  applications  assume  that  planar 
spirals  are  equiangular  and  self -complimentary.  The  applications 
also  assume  that  input  impedance  for  a  two-arm  spiral  is  given  by 
(7.14)  . 
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An  (N)  arm,  rotationally  symmetric  spiral  has  (N-1) 
independent  radiation  modes,  where  each  mode  corresponds  to  a 
different  radiation  pattern.  Typically,  an  (N)  arm  spiral 
operating  in  mode  (M)  excites  each  arm  with  a  current  of 
identical  magnitude.  The  phase  difference  (a)  with  respect  to 
the  first  arm  for  arm  (n)  of  an  (N)  arm  spiral  excited  in  mode 
(M)  is  given  by  [Ref  6:  p.  14-4]: 

g=  iradiars)  (7.15) 

N 

Most  of  the  radiation  from  a  spiral  occurs  at  the  point 
where  the  spiral's  circumference  (C)  equals  the  product  of  the 
mode  and  wavelength  of  the  frequency  of  interest.  As  frequency 
and  wavelength  change,  the  principal  radiation  point  on  the 
antenna  changes  but  the  radiation  parameters  and  beam  patterns  do 
not.  The  shifting  radiation  point  along  the  structure  is  why 
spirals  are  broadcand  antennas.  In  addition,  as  long  as  the 
spiral  is  sufficiently  large  to  radiate  all  desired  frequencies, 
the  shape  of  the  spiral  arm  termination  does  not  effect  the 
antenna's  radiation  pattern  [Ref  6:  pp.  14-4  -  14-7]. 

All  modes  of  a  planar  spiral  whose  center  is  at  the  origin 
in  the  x-y  plane  have  a  null  along  the  z  axis,  unless  mode  one  is 
being  excited.  A  planar  spiral  radiates  in  both  positive  and 
negative  z  half  spaces  with  magnitude  patterns  that  are 
rotationally  symmetric. 

The’  electric  field  of  a  self-complimentary ,  planar  spiral  at 
an  observation  point  in  the  far-field  is  given  by  [Ref  14:  p. 
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530]  : 


E^=E^=0  iV/m) 


(7.16) 


Egk^cos  (0)  [1+jacos  (0)  ]  [tan  ( 

£  = - ± - 

sin(0)2  (7.17) 

j  (M(e..| )  -krff) 

*  ^ (  y/  jjj) 

In  (7.16)  and  (7.17),  (£<,)  is  a  source  strength  constant  and  (k) 

is  the  wavenumber  given  by: 

im-^)  (7.18) 

Rearranging  (7.17),  the  amplitude  of  the  electric  field  (A,) 
becomes : 


A  ^  —  tAn‘*(dCos(0) ) 

cos  (0)  [tan  ( )  1  * 


sin  (0)  vT+a^TcoFT0O7^ 


( V/rrt) 


(7.19) 


The  Mathcad  application  user  should  note  the  phase  variation 
as  a  function  of  azimuth  in  (7.17).  This  phase  variation  can 
result  in  azimuthal  beam  shaping  if  multiple  spirals  are  used  in 
an  array  [Ref  15:  p.  18].  However,  for  a  single  spiral  there  is 
no  change  in  electric  field  magnitude  with  azimuth  and  the 
pattern  is  considered  rotationally  symmetric. 

The  electric  field  computed  using  (7.17)  or  (7.19)  applies 
only  in  the  far-field.  Therefore,  the  following  conditions  must 
all  hold  if  the  far-field  patterns  are  to  be  accurate  [Ref  3:  p. 
92]  : 
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r^jfSl  .61 

(m) 

(7.20) 

im) 

(7.21) 

X 

im) 

(7.22) 

with  the  spiral's  electric  field  given  by  (7.16)  and 
(7.17),  the  antenna's  radiation  intensity  (U) ,  radiated  power 
(Prad)  ,  and  directivity  (D^)  are  found  using  generic  antenna 
formulas  as  follows  [Ref  3:  pp.  28-29]: 

=  (w/ solid  ang)  (7.23) 

Prad=ff^Usin{e)d6di^  (W)  (7.24) 

£>q=  {dimensionless)  (7.25) 

Plad 

In  (7.24),  (n)  is  the  entire  sphere  of  radius  (r^)  in  the  far- 
field  that  surrounds  the  spiral.  In  (7.25),  (U„a,t)  is  the  maximum 

radiation  intensity  anywhere  on  the  sphere. 

The  gain  (G)  of  a  spiral  is  the  product  of  its  directivity 
and  efficiencies.  Conduction-dielectric  efficiency  (6,.^)  of  a 
spiral  cannot  be  easily  determined  and  is,  therefore,  assumed  by 
the  Mathcad  applications  to  be  unity.  Reflection  efficiency 
(€jv),  however,  can  be  approximated.  The  impedance  of  a  two  arm, 
self -complimentary  spiral  is  given  by  (7,14).  In  the  general 
case  of  an  (N)  arm,  self -complimentary  spiral  in  free  space 
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operating  in  mode  (M) ,  the  input  impedance  is  real  and  is  given 
by  [Ref  6 :  p.  14-22 ] : 


z,=w. 


30n 


M 

sin(n^) 

N 


(Q) 


(7.26) 


If  the  characteristic  impedance  (Z^)  of  the  feed  line  is  known  or 
can  be  estimated,  the  voltage  reflection  coefficient  (T) , 
reflection  efficiency,  and  gain  are  computed  by: 


r= 


( dimens i onl ess) 


(7.27) 


e^^=l-|ri^  (dimei  jionless)  (7.28) 

(dimensionless)  (7.29) 

The  planar  spiral  radiates  in  both  directions  normal  to  its 
surface.  Improvements  to  antenna  gu^n  can  be  realized  if  the 
radiation  in  the  undesired  direction  is  reflected  or  eliminated. 

A  common  technique  used  to  achieve  an  improvement  in  gain  is  to 
place  a  cylindrical  metal  cavity  on  the  side  of  the  spiral  that 
has  the  unwanted  beam  pattern.  The  cavity  can  improve  gain  by  up 
to  4.5  dB  but  can  also  reduce  bandwidth  by  up  to  a  factor  of  5. 
The  loss  of  bandwidth  can  be  mitigated  by  filling  the  cavity  with 
electromagnetic  energy  absorbing  material.  The  absorbing 
material  will  reduce  gain  for  spiral  with  a  cavity  by  up  to  1.5 
dB  depending  on  thickness  of  the  material  and  dimensions  of  the 
cavity.  Although  Mathcad  applications  do  not  include  cavity 
parameters,  the  effect  of  cavities  on  antenna  gain  can  be 
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estimated  using  measured  results  from  an  archimedean  cavity- 
backed  spiral.  The  change  of  overall  gain  with  cavity  depth  and 
maximum  gain  with  cavity  diameter  are  pictured  in  Figures  7.2  and 
7.3,  respectively  [Ref  6:  pp.  14-17  -  14-18]. 


FIGURE  7.3  Effect  of  Cavity  Diameter  on  Gain 
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Effective  isotropic  radiated  power  (EIRP)  is  the  product  of 
gain  (dimensionless)  and  power  delivered  to  the  input  of  the 
antenna.  Effective  isotropic  radiated  power  for  a  planar  spiral 
is  given  by: 

(W^)  (7.30) 

If  the  magnitude  of  the  current  (I^)  at  the  antenna's  feed 
is  known,  its  radiation  resistance  (R.)  may  be  expressed  as: 


i?,=  . 


2P 


zad 


(Q) 


(7.31) 


The  spiral  antenna  Mathcad  application  user  should  note  that  the 
magnitude  of  the  transmitter  current  cannot  be  estimated  from  the 
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antenna's  dimensions  and  is  normalized  to  unity. 

The  polarization  of  a  planar  spiral  whose  center  is  the 
origin  is  circular  within  70*  of  the  z  axis.  The  handedness  of 
the  spiral's  polarization  is  the  same  as  the  spiral  if  measured 
in  the  +z  half-space  and  if  the  spiral  is  excited  at  the  central 
feed  point.  If  the  polarization  is  measured  in  the  -z  half-space 
or  if  the  spiral  is  fed  at  its  peripheral  termination,  the 
handedness  of  polarization  is  opposite  that  of  the  spiral.  Dual 
polarization  is  achieved  when  the  structure  is  simultaneously 
excited  at  both  the  center  and  periphery  fRef  6:  p.  14-20]. 

Mathcad  applications  assume  that  a  planar  spiral's  unit 
polarization  vector  (a„)  is  circular,  although  the  user  may 
modify  this  according  to  observed  feed  structures.  If  the 
incoming  wave  unit  polarization  vector  (a„)  can  be  determined, 
the  polarization  loss  factor  (PLF)  of  the  antenna  is  calculated 
as  [Ref  3 :  p.  51] : 

PLf’=  Id^-d'i^  (dimensionless)  (7.32) 

With  all  efficiency  and  loss  terms  computed,  the  antenna's 
maximum  effective  aperture  (A^^)  and  maximum  effective  height 
(hg^)  are  estimated  as  follows  [Ref  3:  p.  63]: 

A,„=e,^^(PLF)  (|1)  (m^)  (7.33) 


h  =2  ' 

\ 


(m) 


(7.34) 


Figure  7.4  illustrates  the  difference  between  the  electric 
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field  pattern  computed  by  the  spiral  antenna  Mathcad  applications 
and  measured  results  for  an  equiangular  spiral  antenna  (a=.3, 
6=90°,  ro=.005  meters,  R=.142  meters)  operating  at  2.8  GHz  over  a 
conducting  plane  [Ref  16:  p.l85]. 


FIGURE  7.4  Equiangular  Spiral  Electric  Field  Patterns 


B.  THE  CONICAL  SPIRAL  ANTENNA 

The  planar  spiral  antenna  offers  many  features  which  make  it 
a  very  popular  device.  However,  the  fact  that  it  radiates  in  a 
direction  normal  to  both  of  its  surfaces  is  a  major  drawback. 
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Although  cylindrical  cavities  placed  on  one  side  of  the  spiral 
can  reduce  the  effects  of  unwanted  radiation,  a  modification  to 
the  basic  planar  structure  can  accomplish  the  same  result  without 
additional  apparatus.  If  the  planar  spiral  is  altered  into  a 
conical  shape,  many  of  the  desirable  features  of  the  planar 
spiral  are  maintained,  but  radiation  occurs  primarily  in 
direction  of  the  cone's  tip. 

Any  point  on  the  edge  of  a  log-conical  spiral  antenna 
may  be  defined  by  [Ref  5:  p.  286]: 

(m)  (7.3  5) 

In  (7.35),  is  the  angular  offset  of  the  i*-*’  edge,  (B^)  is  the 

conical  half-angle,  and  (b)  is  an  arbitrary  constant  given  by: 

jb=cot(P)  (dimensionless)  (7.3C) 

The  geometry  of  a  log-spiral  is  illustrated  in  Figure  7.5. 
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FIGURE  7.5  Conical  Log-Spiral  Geometry 

In  Figure  7.5,  (B)  is  the  overall  diameter  or  twice  the  overall 

radius  (R)  and  (d)  is  twice  the  feed  radius  (rg)  . 

Analogous  to  the  planar  spiral,  the  upper  operating 
frequency  (fhigh)  of  a  conical  log-spiral  is  determined  by  the 
relationship  between  the  wavelength  (Ah^gh)  of  the  upper  frequency 
and  the  spacing  (d)  between  feed  points  as  follows  [Ref  5:  p. 


287]  ; 


d=hll^  (m)  (7.37) 

4 

The  lower  operating  frequency  (flow)  of  the  conical  spiral  is 
determined  by  the  correlation  between  the  lower  frequency's 
wavelength  (Xiow)  and  the  antenna's  overall  diameter  by  [Ref  5:  p. 
287]  : 

B=lh^  (m)  (7.38) 

O 

The  bandwidth  of  the  conical  log-spiral  is  computed  using 
(7.11)  . 

The  conical  log-spiral  radiates  a  single  lobe  in  the 
direction  of  its  apex.  The  pattern  broadens  with  increasing 
spiral  angle  (6)  and  lowering  cone  angle  (20o)  until 
irregularities  occur  and  multiple  beams  begin  to  form.  Figure 
7.6  provides  a  rapid  reference  for  the  boundary  between  usable 
and  unusable  cone  and  spiral  angles.  [Ref  13:  pp.  9-84  -  9-85] 
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Corticoi  UOQ^Soirol  D  i  On  • 


FIGURE  7.6  Conical  Log-Spiral  Useable  Dimensions 


The  far-fieid  electric  field  components  of  a  conical  log- 
spiral  of  total  arm  length  (L)  are  written  as  follows  [Ref  17: 
pp.  321-331]: 

iV/m)  (7.39) 


a  ft  f  L  JiL 

COB  (9)  cos  (9„) 

£'0  =  -jfp^— — —  r  I{^)e  °  A(i)di  {V/m) 

2t  JT  ffK  ^  Q 


^  r  L  COS  (6)  COS  (e_! 

E^=f\k— - I(^)e  °  iV/m) 

^  2r ffO  Jo 


(7.40) 


(7.41) 


In  (7.40)  and  (7.41),  the  slowness  factor  of  the  antenna  (C"  and 
the  total  arm  length  of  a  spiral  arm  (L)  are  constants  defined 


C)=Jl  +  —  {.dimensionless) 


(7.42) 


L=^  (m) 

b 


(7.43) 


In  (7.43),  (0l)  is  the  azimuth  at  the  end  of  the  spiral  arm  and 


^  1  ,  ,  sin(0„)  ,  ,  , 

^L=-^ — r7r-rln( - —)  (radians) 

^  jbsin(e^)  r. 


(7.44) 


The  coefficients  A($)  and  B{5)  of  equations  (7.40)  and  (7.41), 
respectively,  are  expressed  as: 


sin  (6.)  cos  (0)  C(5) 


-sin  (0)  cos  (0^)  ] 


(dimensionless) 


(7.45) 


»  (  -  )  C(^)  (dimensionless) 

2 


(7.46) 


In  (7.45)  and  (7.46),  (a)  and  C(0  are  given  by: 


a  =  —  ( dimensi onl ess) 

N 


(7.47) 


i?sin(0^)  bsin(Q^) 

^Q-j -*-ia)  dimensionless 


(7.48) 


The  (  +  )  sign  in  (7.48)  is  for  C($)  used  in  calculating  (E^)  , 
while  the  (-)  sign  is  for  (E^)  computations. 
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Azimuth  angle  as  a  function  of  distance  along  the  spiral  arm 
(0(0)  irj  (7.45),  (7.46),  and  (7.48)  is  calculated  using: 

4> (O  =  ■  .  ■  In ( -^ -t-l)  {radians)  (7.49) 

jbsin(0^)  ’ 

The  current  distribution  (1(0)  used  in  (7.40)  and  (7.41)  is 
that  of  the  excitation  current  along  the  spiral  arm.  This 
current  distribution  can  only  be  calculated  using  intricate 
numerical  modeling  or  a  method  of  moments  solution.  An 
examination  of  measured  current  distribution  yields  the  following 
engineering  approximation  of  (1(0)  given  a  current  at  the  feed 
point  of  (Iq)  : 

^  {Amps)  (7.50) 

Figure  7.7  provides  a  comparison  of  the  current  distribution 
computed  by  the  conical  log-spiral  Mathcad  applications  using 
(7.50)  to  that  measured  for  a  four  arm  conical  log- spiral 
(00=10°,  5=20°,  b=.46)  operating  in  mode  one.  The  application 
user  should  note  that  most  of  the  radiation  from  the  conical  log- 
spiral  occurs  at  the  location  on  the  device  where  arm  length  (L) 
in  wavelengths  approximately  equals  the  mode  number  [Ref  17:  pp. 
321-331] . 
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FIGURE  7.7  Conical  Log-Spiral  Current  Distributions 


Mathcad  spiral  antenna  application  users  should  note  the 
precise  evaluation  of  the  integrals  in  (7.40)  and  (7.41)  is  very 
time  consuming,  even  on  a  33  MHz,  386  personal  computer.  Thus, 
the  application  user  is  provided  a  trapezoidal  approximation  to 
evaluate  the  conical  log-spiral  antenna  integrals.  Also,  both 
the  exact  and  trapezoidal  applications  only  analyze  electric 
fields  and  other  parameters  in  the  +z  half-  plane. 

In  the  trapezoidal  approximation,  the  integral  is  replaced 
by  a  summation  and  (d$)  is  replaced  by: 


=  ^  (in) 


(7.51) 


In  (7.51),  (t)  is  an  operator  entered  number  of  spiral  arm 

increments.  The  application  user  should  expect  the  Mathcad 
evaluation  of  a  conical  log-spiral  using  the  trapezoidal  method 
to  take  about  25%  of  the  time  required  for  exact  method 
predictions . 

Conical  log-spiral  electric  field  calculations  are  valid 
only  if  the  conditions  of  (7.20)  through  (7.22)  are  satisfied. 
Mathcad  applications  assume  that  overall  radius  is  the  largest 
dimension  of  the  conical  log-spiral. 

Radiation  intensity  of  the  conical  log-spiral  antenna  is 
determined  by  [Ref  3:  p.  28]: 

+  (W/ solid  ang)  (7.52) 

2t1o 

The  radiated  power  and  directivity  of  the  conical  log-spiral 
are  given  by  (7.24)  and  (7.25).  Users  of  the  Mathcad  spiral 
antenna  application  should  be  aware  that  computer  solution  of 
these  formulas  may  be  very  time  consuming.  However,  one  may 
reasonably  consider  most  conical  log-spiral  antennas  in  the 
useable  region  of  Figure  7.6  to  have  a  total  electric  field  that 
is  rotationally  symmetric  with  respect  to  the  antenna's  axis. 
Given  this  fact,  an  approximation  of  directivity  may  be  made.  An 
estimate  of  half-power  beamwidth  (A0)  in  degrees  as  a  function  of 
spiral  and  cone  angles  is  provided  in  Figure  7.8  below. 
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FIGURE  7.8  Conical  Log-Spiral  Half_Power  Beamwidths 


The  half-power  beamwidth  from  Figure  7.8  is  used  to  estimate 
directivity  by  [Ref  13:  p.  9-87]: 


, 32600 
”  (Ae)2 


(dS) 


(7.53) 


It  is  impossible  to  accurately  determine  the  conduction- 
dielectric  efficiency  of  a  conical  log-spiral  antenna  based 
solely  on  measured  geometry  of  the  antenna.  Hence,  the  Mathcad 
conical  log-spiral  applications  assume  all  conduction  dielectric 
efficiency  is  unity.  The  input  impedance  (ZJ  of  a  conical  log- 
spiral  antenna  as  a  function  of  angular  arm  width  is  provided  in 
Figure  7.9. 


89 


FIGURE  7.9  Input  Impedance  of  a  Conical  Log-Spiral 


Using  the  input  impedance  from  Figure  7.9  and  an  estimate  of  the 
characteristic  impedance  of  the  feed  assembly,  the  Mathcad 
conical  log-spiral  antenna  application  user  may  compute  voltage 
reflection  coefficient,  reflection  efficiency,  gain,  and 
effective  isotropic  radiated  power  using  (7.27),  (7.28),  (7.29), 

and  (7.30),  respectively. 

If  the  polarization  unit  vector  (a„)  of  an  incoming  wave  is 
known,  precise  determination  of  polarization  loss  at  a  point  in 
the  far-field  may  be  desired.  To  accomplish  this,  Mathcad 
applications  convert  electric  field  components  of  (7.40)  and 
(7.41)  to  Cartesian  coordinates  at  a  user  defined  point  (x,y,z) 
and  compute  the  antenna's  unit  polarization  vector  (aj,)  as 


follows  [Ref  4:  pp.  35-36,  364-367]: 


E^=EqCOs  (0)  sin  (<|>)  -f^sin  (<t>) 

(V/m) 

(7.54) 

Ey^E^cos  (0)  sin  (<t»)  *E^cos  ((t>) 

(V/m) 

<7.55) 

£2  =  -£eSin  (0)  (V/m) 

(7.56) 

Of,  {x.  V.  z)  =  (dimensionless) 

y/\E(x,y,  z)  |2 


(7.57) 


The  conical  log-spiral  antenna^s  polarization  loss  factor  is 
calculated  using  (7.32).  Using  the  polarization  loss  factor, 
maximum  effective  aperture  (A,„)  and  maximum  effective  height 
(hem)  for  a  conical  log-spiral  are  given  by: 


(m^) 


(7.58) 


(7.59) 


Figure  7.10  and  Table  7.1  compare  measured  data  to  that 
calculated  by  the  Mathcad  applications  for  a  two  arm,  mode  one 
conical  log-spiral  (b=.053,  0,,=1O°,  0=73°,  d=.03  meters,  B=.30 
meters,  «S=90°)  operating  at  350  MHz  [Ref  18:  p.  332]. 
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FIGURE  7.10  Conical  Log-Spiral  Electric  Field  Pattern 


TABLE  7.1  Conical  Log-Spiral  Data  Comparison 


VIII.  THE  CONICAL  HORM  AMTENNA 

The  conical  horn  antenna  is  a  device  that  provides  a 
transition  for  an  electromagnetic  wave  from  a  circular  waveguide 
to  an  unbounded  medium  such  that  the  wavefront  phase  at  the 
horn's  aperture  is  nearly  constant  [Ref  2:  pp.  644-645].  As  a 
result  of  its  axial  symmetry,  the  conical  horn's  radiation 
pattern  is  strictly  a  function  of  the  cone's  geometry  for  a  given 
mode  of  excitation  [Ref  10:  p.  10-3].  The  conical  horn  is  widely 
employed  as  a  feed  element  for  reflector  assemblies  used  in 
satellite  tracking,  microwave  communications,  and  radar.  The 
geometry  of  a  conical  horn  is  illustrated  in  Figure  8.1. 


y 


HORN 


FIGURE  8.1  Conical  Horn  Geometry 

In  Figure  8.1,  (a)  is  the  inner  radius  of  the  circular  waveguide, 

flare  angle  (a)  is  the  included  angle  of  the  horn,  (d)  is  the 


93 


diameter  of  the  mouth  of  the  horn,  and  the  axial  height  (h)  is 
the  distance  from  the  origin  to  the  center  of  the  mouth  of  the 
horn. 

A  circular  waveguide  will  only  propagate  a  transverse, 
electric  (TE)  or  transverse  magnetic  (TM)  mode  of  an 
electromagnetic  wave  if  the  frequency  of  the  wave  is  above  a 
minimum  value  for  the  mode  called  the  cutoff  frequency  (fj.  The 
propagating  mode  with  the  lowest  cutoff  frequency  is  called  the 
dominant  mode.  For  circular  waveguides  the  (TE,i)  mode  is  the 
dominant  mode  [Ref  4:  p.570].  Mathcad  conical  horn  antenna 
applications  assume  that  the  waveguide  is  excited  only  in  the 
dominant  (TEn)  mode.  Also,  the  applications  do  not  compute 
bandwidth  for  conical  horn  antennas.  Rather,  the  software 
computes  cutoff  frequencies  for  selected  modes  such  that  the  user 
can  determine  if  the  waveguide  will  support  propagation  of  a 
specific  mode. 

Each  mode's  cutoff  frecjuency  is  a  function  of  the  circular 
waveguide's  inner  radius.  For  transverse  electric  waves,  the 
cutoff  frequency  of  an  air  filled  circular  waveguide  is  given  by 
[Ref  19:  pp.  472-473]: 

^Hz)  (8.1) 

In  (8.1),  (Mo)  is  the  permeability  of  free  space,  (e^)  is 
permittivity  of  free  space,  and  (x'mn)  is  zero  of  the 

derivative  with  respect  to  the  argument  of  the  Bessel  function  of 
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the  first  kind,  order  m.  (x'mn)  obtained  from  Table  8.1. 


TABLE  8.1  Zeroes  of  the  Bessel  Function  Derivative 


y  ' 

A  mn 

m=0 

m=l 

m=2 

m=3 

m=4 

n=l 

3.8318 

1.8412 

3.0542 

4.2012 

5.3175 

n=2 

7.0156 

5.3315 

6.7062 

8.0153 

9.2824 

n=3 

10.1735 

8.5363 

9.9695 

11.3459 

12.6819 

n=4 

13.3237 

11-7060 

13.1704 

14.5859 

15.9641 

The  cutoff  frequency  of  transverse  magnetic  waves  in  an  air 
filled  circular  waveguide  is  written  as  [Ref  19:  pp.  478-479]: 


l«z>  (8.2) 

In  (8.2),  (Xmn)  is  the  n'’^  zero  of  the  Bessel  function  of  the 
first  kind,  order  m.  (Xmn)  be  obtained  from  Table  8.2. 


TABLE  8.2  Zeros  of  the  Bessel  Function 


Xrnn 

m=0 

m=l 

m=2 

m=3 

m=4 

n=l 

2.4049 

3.8318 

5.1357 

6.3802 

7.5884 

n=2 

5.5201 

7.1056 

8.4173 

9.7610 

11.0647 

n=3 

8.6537 

10.1735 

11.6199 

13.0152 

14 . 3726 

n=4 

11.7915 

13 . 3237 

14.7960 

16.2235 

17 . 6160 
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A  conical  horn  is  said  to  be  optimum  if  the  diameter  of  the 
mouth  of  the  horn  satisfies  the  following  [Ref  10:  p.  10-9]: 

d=-|Xsin(-|)  (m)  (8.3) 

To  allow  for  comparison  between  the  antenna  being  evaluated  and 
an  optimum  horn,  the  conical  horn  Mathcad  applications  calculate 
optimum  (d)  for  a  user  provided  frequency  of  interest  and  horn 
flare  angle. 

The  electric  field  calculations  of  the  conical  horn  Mathcad 
applications  use  the  magnetic  field  integral  equation  solution 
for  aperture  antennas.  This  is  similar  to  the  combined  field 
integral  equation  solution  first  postulated  by  Schorr  and  Beck 
for  conical  horns  in  1950  [Ref  20:  p.  795].  The  application  user 
can  expect  a  large  amount  of  computer  processing  time  will  be 
necessary  to  analyze  a  conical  horn  and  that  this  time  will  grow 
with  the  square  of  the  number  of  increments  (i)  into  which  the 
far-field  is  divided. 

In  order  to  compute  the  magnetic  vector  potential  (A)  of  the 
aperture  field  of  a  conical  horn,  several  preliminary  functions 
and  related  constants  must  first  be  defined.  These  include 
associated  Legendre  functions  of  the  first  kind  of  order  (v) 
[P^tmjcos  (0)  ]  ,  spherical  Hankel  functions  of  the  second  kind  of 
order  (v)  [h^,‘^’ (kr)  ]  ,  and  the  derivative  with  respect  to  (kr)  of 

the  spherical  Hankel  functions  [h'„'^’ (kr)  ] .  The  arguments  of  the 
spherical  Hankel  functions  are  (r) ,  the  distance  from  the  origin 
to  an  observation  point  in  the  far-field,  and  the  wavenumber  (k) 
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for  a  wavelength  (A.)  corresponding  to  a  frequency  of  interest. 
Wavenumber  is  given  by: 


(8.4) 


The  order  (v)  of  the  Legendre  and  Hankel  functions  may  be 
approximated  using  the  constant  (b^)  as  follows  [Ref  21:  p.  521] 


(1--) 


log  [■ 


{dimensionless) 


(8.5) 


log  [cos  ( ] 


V  =  - . 5  +  . 


( dimens i onless) 


(8.6) 


The  associated  Legendre  function  of  the  first  kind  with  m=l 
may  be  estimated  using  the  gamma  function  (r(z))  for  angles  less 
than  (7r/3)  as  follows  [Ref  22:  pp.  336]: 

(cos  (9)  )  =  [.  4sin(e)]cos((v^.5)e^^) 

r(v+1.5)\2  4  (8.7) 

( dimens i onl ess) 

The  value  of  the  spherical  Hankel  function  and  its 
derivative  with  respect  to  its  argument  in  the  far-field  may  be 
estimated  by  [Ref  17  p.  796]: 


hr'  (kr) 

kr 


( dimens i onless) 


(8.8) 


'  (kr) 


-T^-e  ^  ^  ’  [1-J  ]  {dimensionless)  (8-9) 

kr  kr 


The  components  of  the  magnetic  vector  potential  are  given  by 
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[Ref  20:  p.798] : 


♦  [sin (56' )  P,  (cos  (6^) )  cos^  (<|)' )  cos  O' ) 

^_6co£O^^P  (cos  (6  ) )  sin^  (4>' )  sin^  O' )  1 
sinO')  ° 

♦  d0 '  d<t> '  (Wb/w) 


(8.10) 


gjJcAicos  (p) 

*  [sin(60 ' )  (cos  (0q)  )  cos  ((J**  )  sin  (<|)  ' )  cos  (0 ' ) 

•  (cos  (6^) )  sin(4>' )  cos  (0' )  sin^  (6' )  ] 


sin (0 ' ) 


*d6 '  d4>’  (Wb/m) ) 


(8.11) 


‘gjJcAcos  (p) 


*  [sin  (50 ' )  P,  (cos  (0^) )  cos  (<|)' )  sin  (0 ' )  ]  d0 '  d<J)' 

( Wb/m) 


(8.12) 


The  integrals  in  (8.10)  through  (8.12)  are  performed  over  the 
aperture  of  the  horn,  thus  the  primed  components  in  these 
integrals  indicate  source  coordinates.  The  terms  cos(B),  (6), 
and  (Bh)  in  the  magnetic  vector  potential  integrals  are  given  by: 


cos  (p)  =cos  (6)  cos  (0' )  +sin(0)  sin(0' )  cos  (4>-(t>' ) 

(dimensionless) 


(8.13) 


6  =  —  (dimensionless)  (8.14) 

a 

fl  - e-^‘^^[ihJ^Ukh)  '  (kh)]  (wb/m)  (8.15) 

j(2n}^fe^r  h 

In  (8.15),  (Bq)  is  an  amplitude  scaling  constant.  The  magnetic 
vector  potential  integrals  contained  in  (8.10)  through 
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(8.12)  include  approximations  of  the  associated  Legendre  function 
and  its  first  derivative  with  respect  to  cos (6')  as  follows  [Ref 
20:  p.798]: 

cos  (0 ' )  =sin  (60 ' )  Pv  )  {dimensionless)  (8.16) 


,  (1) 


cos  (0 '  ) 


6cos (50 ' ) 
sin  (0 ' ) 

( dimensionless) 


Pj^’cos  (-|) 


(8.17) 


Far-field  electric  field  components  are  approximated  by 
converting  the  magnetic  vector  potential  components  to  spherical 
coordinates  and  multiplying  by  (k^)  as  follows  [Ref  3:  pp.  92, 
102  j  : 


Ae=i^^cos  (0)  cos  (<t)) +AyCos  (0)  sin (({>) -A^sin (0)  iWb/m)  (8.18) 


i^^=-.A^sin (<{)) +^yCos  (<}>)  {Wb/m)  (8.19) 


£’e=Jc%  (V/m)  (8.20) 

E^=k-A^  (V/m)  (8.21) 

For  the  electric  field  components  of  (8.20)  and  (8.21)  to  be 
valid,  the  observation  point  (r)  must  be  in  the  far-field. 
Therefore,  assuming  the  diameter  of  the  mouth  of  the  horn  is  the 
largest  dimension  of  the  cone,  all  of  the  following  formulas  must 
be  satisfied  [Ref  3:  p.92]: 
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ril.eX  (m) 


(8.22) 


r^5d  im) 

(8.23) 

(m) 

(8.24) 

When  investigating  the  parameters  of  a  conical  horn  antenna 

with  the  magnetic  field  integral  equation,  it  must  be  presumed 

that  no  electric  field  exists  behind  the  mouth  of  the  horn. 

Thus,  radiation  intensity  (U)  ,  radiated  power  (P^ad)  /  and 

directivity  (D^)  of  the  conical  horn  antenna  are  calculated  by 

applying  generic  antenna  formulas  as  follows: 

[\Eq\^*\E^\^]  {w/ solid  ang) 

(8.25) 

P,ad=f^’'r^t/sin(0)  d0ci4>  (f'O 

Jo  Jo 

(8.26) 

^o- —  (dimensionless) 

Eiad 

(8.27) 

In  (8.25),  (rjo)  is  the  intrinsic  impedance  of  free  space  and  in 
(8.26),  (U„„)  is  the  maximum  radiation  intensity  from  the 

antenna.  For  a  conical  horn,  it  may  be  reasonably  assumed  that 
maximum  radiation  intensity  occurs  along  the  antenna's  axis 
(i.e. ,  0=0,  0=0) . 

In  order  to  lower  the  computer  processing  time  of  the 
conical  horn  Mathcad  applications,  the  user  may  choose  a  file 
that  uses  a  trapezoidal  approximation  of  the  integrals  in  (8.10) 
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through  (8.12)  and  (8.26).  The  user  may  also  vary  the  number  of 
trapezoidal  increments  (tl,  t2,  t3)  to  vary  the  extent  of  this 
approximation  and  to  adjust  computation  time. 

The  ohmic  losses  of  conical  horns  and  their  associated 
circular  waveguides  are  very  difficult  to  precisely  determine, 
but  are  normally  very  small.  Li)cewise,  reflection  efficiency  of 
a  conical  horn  cannot  be  easily  found  using  analytical 
techniques,  but  when  measured  for  horns  with  moderate  flare 
angles  and  high  directivity  it  is  usually  close  to  unity.  The 
conical  horn  antenna  Mathcad  applications  approximate  the  product 
of  ohmic  and  reflection  efficiencies  as  0.95  and  compute  antenna 
gain  (G)  and  effective  isotropic  radiated  power  (EIRP)  by  [Ref 
10:  p,  10-9]: 

G=.9  5£>a  (.dimensionless)  (8.28) 


(W)  (8.29) 

The  unit  polarization  vector  (a^)  in  Cartesian  coordinates 
of  a  conical  horn's  wave  may  be  computed  at  any  point  in  the  far- 
field  as  follows  [Ref  4:  pp.  35-40,  364-367]: 


E^=EqCos  (0)  cos  (<t))  -E^sin  (4>) 

( V/m) 

(8.30) 

Ey^E^cos  (0)  sin  (<t>)  +E^cos  ((J>) 

( V/m) 

(8.31) 

E^=-EQSin(e)  (V/m) 

(8.32) 

_  <3^  +a  £  +<3  E 

oAx.y.z)  (dimensionless) 

^l\E(x,y,  z) 


(8.33) 
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If  an  incoming  wave's  unit  polarization  vector  (a„)  is  known 
or  can  be  estimated,  the  polarization  loss  factor  (PLF)  of  the 
antenna  at  a  point  in  the  far-field  can  be  expressed  as  [Ref  3: 
p.  51]: 

PLF=|d^-c’|2  (dimensionless)  (8.34) 

The  maximum  effective  aperture  of  a  conical  horn  may  be 
determined  as  follows  [Ref  3:  p.63]: 

.9SDAPLF)  (~)  im'^)  (8.35) 

The  actual  value  of  current  (Iq)  applied  to  the  input  of  a 
conical  horn  cannot  be  determined  by  dimensional  information 
alone.  Thus,  exact  calculation  of  radiation  resistance  (Rj.)  is 

impossible.  However,  if  input  current  is  normalized  to  one  amp, 
normalized  radiation  resistance  and  maximum  effective  height 
(hgm)  may  be  written  as  [Ref  2:  p.  42]: 

2  P 

Rr  =  -r^  (CJ)  (8.36) 

l-fol 


(8.37) 


A  conical  horn  is  a  member  of  a  group  of  devices  known  as 
aperture  antennas.  A  term  frequently  used  to  analyze  the 
performance  of  aperture  antennas  is  aperture  efficiency  (Cgp). 
Aperture  efficiency  is  the  ratio  of  maximum  effrictive  aperture  to 
physical  area  at  the  mouth  of  the  horn  and  is  calculated  by  [Ref 
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3:  p.  475] 


e,_= - 3! —  {dimensionless)  ,0  -iqn 

^  d)2  (8.3B) 

2 

Aperture  efficiency  of  a  conical  horn  is  typically  about  50%. 

Table  8.3  and  Figure  8.2  compare  measured  data  to  that 
calculated  by  the  Mathcad  applications  for  a  conical  horn  (a=.045 
meters,  h=1.489  meters,  d=1.0  meters,  and  0=20°)  operating  at 
1.96  GHz  [Ref  22:  p.lOO]. 


TABLE  8.3  Conical  Horn  Data  Comparison 


ANTENNA 

MEASURED 

CALCULATED 

PARAMETER 

DATA 

DATA 

GAIN 

22.0  dB 

23.6  dB 

O  o  rm  p  o  «o  of  Oonicol  Mor-m  Ant©ririo  Electric  Eleips 


FIGURE  8.2  Conical  Horn  E-Plane  Electric  Field  Patterns 
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IX.  THE  PYRAMIDAL  HORM  ANTENNA 

The  most  popular  member  of  the  horn  family  of  antennas  is 
the  pyramidal  horn.  The  pyramidal  horn  provides  a  transition 
from  a  rectangular  waveguide  to  free  space  and  is  flared  both 
horizontally  and  vertically.  It  possesses  a  radiation  pattern 
that  is  essentially  the  combination  of  that  of  the  E-  and  H-plane 
sectoral  horns.  If  accurately  constructed,  the  pyramidal  horn 
offers  the  ability  to  control  beamwidth  in  both  principal  planes 
with  a  gain  (G)  that  closely  matches  theoretical  predictions  [Ref 
10:  p.  10-3].  The  geometry  of  the  pyramidal  horn  in  the  E-plane 
{<p=n/2)  and  H-plane  (0=0)  is  illustrated  in  Figures  9.1  and  9.2, 
respectively. 


FIGURE  9.1  E-Plane  Pyramidal  Horn  Geometry 
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RECTANGULAR 

WAVEGUIDE 


MOUTH 
OF  HORN 


FIGURE  9.2  H-Plane  Pyramidal  Horn  Geometry 


In  Figures  9.1  and  9.2,  (a)  and  (b)  are  the  dimensions  of 
the  rectangular  waveguide  used  to  excite  the  horn,  (aj  and  (b^) 
are  the  dimensions  of  the  mouth  of  the  horn,  (ilJe.h)  ^re  half  the 
flare  angles  in  the  indicated  planes,  (Pi,2,e,h)  the  distances 

from  the  imaginary  apex  of  horn  to  the  indicated  points  on  the 
perimeter  of  the  mouth  of  the  horn,  and  (Peh)  are  the  distances 
from  the  beginning  of  the  horn's  flare  to  the  center  of  the  mouth 
of  the  horn  along  the  indicated  axis. 

To  physically  construct  a  pyramidal  horn,  the  parameters 
(Pe)  and  (Ph)  should  be  equal.  These  dimensions  are  calculated 
as  follows  [Ref  3:  p.  568]: 
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Pe=  (b^-b) 


i3l  4 


(m) 


(9.1) 


p^= (a^-a) 


(m) 


(9.2) 


The  Mathcad  pyramidal  horn  applications  do  not  calculate  a 
bandwidth  of  operating  frequencies.  Rather,  the  applications 
compute  a  matrix  of  transverse  electric  (TE)  and  transverse 
magnetic  (TM)  cutoff  frequencies  (f^)  based  on  the  dimension  of 
the  waveguide.  The  cutoff  frequencies  are  the  lowest  frequencies 
of  a  given  mode  which  can  propagate  in  the  waveguide.  Cutoff 
frequency  for  transverse  electric  or  transverse  magnetic  wave  of 
mode  (m,n)  in  an  air  filled  waveguide  is  given  by  [Ref  4:  p. 

549]  ; 


(f. 


2V^ 


(->'*(-2)=  Wz) 
a  b 


(9.3) 


In  (9.3),  (Mo)  is  the  permeability  of  free  space,  (e^)  is  the 
permittivity  of  free  space,  and  (m)  and  (n)  are  integers.  For 
transverse  magnetic  waves,  (m)  and  (n)  must  be  non-zero.  For 
transverse  electric  waves,  either  (m)  or  (n)  may  be  zero,  but  not 
both.  The  mode  with  the  lowest  cutoff  frequency  is  called  the 
dominant  mode.  The  (TE^o)  mode  is  dominant  in  a  rectangular 
waveguide . 

Assuming  the  electric  fields  behind  the  mouth  of  the  horn 
are  zero  and  the  physical  dimensions  of  the  horn  are  negligible 
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in  the  far-field,  one  may  express  the  electric  field  components 
of  a  pyramidal  horn  in  the  +2  half-plane  as  [Ref  3:  pp.  565-578]: 


a  -jkz 

E^=jkE— - [sin(<t>)  (1+cos  (0)  )  Jj^l2]  (V/w)  (9-4) 

”  ®  4nr 

a  -jki 

E.=jkE^— [cos  (4>)  (1+cos  (0) )  (V/irt)  (9.5) 

4'itr 

In  (9.4)  and  (9.5),  (E^)  is  an  electric  field  amplitude  scale 

factor  set  to  unity  by  the  Mathcad  applications,  (r)  is  the 
distance  from  the  origin  to  an  observation  point  in  the  far-field 
and  (k)  is  the  wavenumber  corresponding  to  the  wavelength  (X)  of 
a  frequency  (f)  of  interest  and  is  given  by: 

k=^  (9.6) 

The  functions  (IJ  and  (I2)  in  (9.4)  and  (9.5)  are  computed  by: 


cos ( —  5) e  ^ 


(sln(6)cos(^)  ] 

* 


{dimensionless)  (9.7) 


- -jk( -4— -{sin(0)  sin(^) ) 

l2‘(  |.e 

2 


{dimensionless) 


(9.8) 


The  Mathcad  application  user  should  note  that  closed  form 
solutions  of  (9.4)  and  (9.5)  exist  but  are  very  complex  and 
tedious.  Thus,  the  applications  use  a  numerical  approximation  of 
these  integrals  to  compute  electric  field  components  for  the 
horn. 

The  electric  field  components  of  (9.4)  and  (9.5)  are 
applicable  only  in  the  far-field  of  the  horn.  Therefore,  the 
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distance  to  the  far-field  observation  point  must  satisfy  all  of 
the  following  [Ref  3:  p.92]: 

ril.eX  (m)  (9.9) 

ri5£>  (/n)  (9.10) 

21)2 

im)  (9.11) 

In  (9.10)  and  (9.11) ,  (D)  is  the  distance  between  opposite 
corners  of  the  mouth  of  the  horn. 

The  radiation  intensity  (U)  ,  radiated  power  (P^ad)  / 
directivity  (Dq)  ,  and  effective  isotropic  radiated  power  (EIRP) 
of  a  pyramidal  horn  are  calculated  using  generic  antenna 
equations  as  follows: 


{W/solidang)  (9.12) 

^rad^/Z^^sinO)  ded<l>  (fiO  (9.13) 

^o=~  p  (dimensionless)  (9.14) 

^rad 

EIRP=P^^^^  (W)  (9.15) 


In  (9.12),  (ho)  is  the  intrinsic  impedance  of  free  space,  in 
(9.13),  (n)  is  the  hemisphere  of  radius  (r)  in  the  +z  half-space, 
and  in  (9.14)  ,  (Un,aj()  is  the  maximum  radiation  intensity  anywhere 
on  that  hemisphere. 

The  efficiency  of  a  pyramxGai.  .»orn  is  difficult  to  predict 
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based  on  measured  geometry  alone.  However,  the  typical  horn  can 
reasonably  be  assumed  to  be  approximately  50%  efficient.  Thus, 
Mathcad  pyramidal  horn  applications  compute  gain  as  follows; 

G=.5Do  (.dimensionless)  (9.16) 

The  unit  polarization  vector  (a,)  of  a  pyramidal  horn's 
electromagnetic  wave  at  a  point  in  the  far-field  is  determined 
using  the  electric  field  components  from  (9.4)  and  (9.5).  The 
antenna's  unit  polarization  vector  in  Cartesian  coordinates  is 
computed  by  [Ref  4:  pp.  34,  364-367]: 


£^=£■0003  (0)  cos  (<j))  -E^sin  (4>) 

(V/m) 

(9.17) 

£^=£0003  (0)  sin  (4))  +£^cos  (<|>) 

(V/m) 

(9.18) 

£2=-£0sin(0)  (V/m) 

(9.19) 

a  ^(x,  V.  2)  =  (dimensionless) 

(9.20) 

v'|E(x,y,  z)  j 

If  an  incoming  wave's  unit  polarization  vector  (a„)  is  known 
or  can  be  estimated,  the  pyramidal  horn's  polarization  loss 
factor  (PLF)  can  be  detemnined  by  [Ref  3:  p.  51]: 

PLF=\d^'al\'^  (dimensionless)  (9.21) 

The  maximum  effective  aperture  (Aen,)  and  aperture  efficiency 
(Cap)  of  the  pyramidal  horn  is  written  as  [Ref  3:  p.63]: 


e^P= — ^  (dimensionless)  (9.23) 

The  amplitude  of  the  input  current  (1^)  used  to  excite  the 
rectangular  waveguide  cannot  be  calculated  based  on  dimensional 
information  alone.  Thus,  in  order  to  compute  radiation 
resistance  (Rj.)  and  maximum  effective  height  (h^)  ,  a  normalized 
value  of  1  amp  is  assumed  to  excite  the  waveguide.  These 
parameters  are  expressed  as  [Ref  2:  p.  42]: 

(Q)  (9.24) 

Uol 


(m) 


(9.25) 


The  Mathcad  pyramidal  horn  applications  may  also  be  used  to 
analyze  E-plane  and  H-plane  sectoral  horns.  The  term  (i^/2p)  in 
the  exponents  of  the  integrals  of  (9.7)  and  (9.8)  is  a  phase 
error  term  that  accounts  for  differences  in  phase  between  the 
center  and  any  point  in  the  aperture  of  the  horn.  The  sectoral 
horns  are  evaluated  by  eliminating  the  phase  error  term  in  the 
direction  that  is  not  flared.  To  accomplish  this  modification  to 
(9.7)  and  (9.8),  the  application  user  is  directed  to  set  (al) 
ecfual  to  (a)  for  E-plane  sectoral  horn  analysis  and  (bl)  equal  to 
(b)  for  H-plane  sectoral  horn  analysis.  If  the  application  user 
makes  these  selections,  (IJ  for  the  E-plane  sectoral  horn  or 
(I2)  for  the  H-plane  sectoral  horn  is  altered,  respectively,  as 
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follows  [Ref  3:  pp.  536,  552]: 


COS  (-^sin(0)  cos  (4>)  ) 

( (^sin(e)cos(4>) 
(dimensionless) 


(9.26) 


sin  sin  (0)  sin(4>)  ) 

J2=[ - ^ - ]  (dimensionless)  (9.27) 

-^sin(0)  sin(4>) 

Table  9.1,  Figure  9.3,  and  Figure  9.4  compare  measured  data 
to  that  calculated  by  the  Mathcad  applications  for  a  pyramidal 
horn  (Pi=.3398  meters  P2=.3198  meters,  ai=.1846  meters,  bi=.1455 
meters,  a=. 02286  meters,  b=. 01016  meters)  operating  at  9.3  GHz 
[Ref  5:  pp.  413-415] . 


TABLE  9.1  Pyramidal  Horn  Data  Comparison 


ANTENNA 

MEASURED 

CALCULATED 

PARAMETER 

DATA 

DATA 

DIRECTIVITY 

21.3  dB 

21.6  dB 
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Electric  Field  in  V/m  >4  Electric  Field  in  V/m 


IGURE  9.3  E-Plane  Pyramidal  Horn  Electric  Field  Pattern 


FIGURE  9.4  H-Plane  Pyramidal  Horn  Electric  Field  Pattern 
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X.  REMARKS  AND  COWCLDSIOW8 


The  intent  of  this  thesis  and  associated  software  was  to 
provide  the  Naval  Maritime  Intelligence  Center  (NAVMARINTCEN) 
with  a  relatively  simple,  user  friendly  set  of  Mathcad 
applications  that  would  analyze  various  types  of  antennas  based 
solely  on  dimensional  information  and  ground  characteristics. 
Although  the  software  format  achieves  the  goal  of  easy  use,  the 
nature  of  the  formulas  applicable  to  many  of  the  antennas 
necessarily  reduces  the  simplicity  of  the  programs.  Indeed,  many 
of  the  equations  used  in  this  research  project  are  so  complex 
that  they  are  not  found  in  any  existing  textbook  dealing  with 
antennas.  Nevertheless,  the  Mathcad  applications  developed  in 
conjunction  with  this  thesis  allow  the  user  to  analyze  several 
antennas  that  cannot  be  studied  with  current  electromagnetics 
software  packages  such  as  ELNEC,  NEC,  or  WIRE.  And,  in  every 
case,  the  Mathcad  applications  compute  antenna  parameters  and 
far-field  radiation  patterns  that  closely  compare  with  measured 
or  predicted  results. 

The  Mathcad  applications  are  compatible  with  any  personal 
computer  that  supports  Mathcad  3.1  for  Windows.  However,  many  of 
the  programs  require  a  numerical  solution  of  highly  complicated 
integral  equations  that  are  computationally  intensive.  Where 
necessary  to  reduce  processing  time,  the  applications  use 
trapezoidal  approximations  as  an  alternative  numerical  procedure 
to  evaluate  required  integrals.  Nevertheless,  several  programs 
still  require  several  days  to  complete  calculations  on  a  33  MHz, 


113 


386  personal  computer. 

As  previously  mentioned,  several  of  the  antennas  included  in 
this  project  are  not  adequately  reviewed  by  current  texts.  Thus, 
equations  from  many  professional  journals  and  doctoral 
dissertations  are  used  for  a  number  of  the  applications. 
Unfortunately,  many  errors  existed  in  these  sources,  and 
resolution  of  these  mistakes  significantly  slowed  the  progress  of 
our  research. 

The  most  disappointing  aspect  of  the  Mathcad  application 
software  aside  from  the  excessive  length  of  time  required  to 
analyze  some  types  of  antennas  is  the  graphical  presentation  of 
the  far-field  radiation  patterns.  Mathcad' s  limited  graphics 
flexibility,  particularly  in  regards  to  three-dimensional  and 
spherical  coordinate  plots,  precluded  better  presentation  of 
output  data. 

The  equations  and  data  of  this  antenna  analysis  package 
could  easily  be  transferred  to  another  mathematics  program,  such 
as  MATLAB.  Although  a  MATLAB  program  may  not  be  as  easy  for  a 
new  user  to  employ,  it  might  offer  advantages  with  respect  to 
lower  processing  time  and  improved  graphical  output. 
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THE  HELICAL  ANTENNA 
MATHCAD  SOFTWARE-HELIX. MCD 

When  built  to  the  proper  specifications,  the  helical  antenna  possesses 
many  qualities  which  make  it  suitable  for  a  wide  variety  of 
communications  applications.  If  the  following  conditions  are  satisfied 
the  he-ix  will  exhibit  a  highly  directional  axial  main  looe,  low  side 
lobe  level,  negngible  mutuil  interference  with  adjacent  antennas,  low 
voltage  standing  wave  ratio  (VSWR)  ,  and  resistive  input  i.mpedance  over  a 
wide  frequency  band: 

.8  <  C;^  <  1 .  15 

n  >  3 

12  <  a  <  14 

(Note:  X  in  a  subscript  indicates  the  dimension  is  in  wavelengths. 

Mathcad  equations  can  nor  use  symbolic  subscripts.  Therefore,  the  symbol 
X  will  immediately  follow  the  parameter  in  equations  (i.e.,  CX)  ). 

The  helical  antenna  Mathcad  application  will  compute  the  following 
parameters  (Items  with  *  indicate  parameters  that  are  calculatea  for 
boti*:  axial  and  peripheral  feed  geometries)  : 

C  =  circumference  of  helix 
X  =  wavelength 
a  =  pitch  angle 
Eq  =  directivity 

p  =  relative  phase  velocity 
E0^  =  Electric  Field  Components 

'f  =  Array  Factor  Phase  Shift 
U  =  Radiation  Intensity 
=  Radiated  Power 

R  =  Antenna  Input  Resistance* 
r  =  Voltage  Reflection  Coefficient* 

Cj-  =  Reflection  Efficiency* 

^em  ~  Maximum  Effective  Height* 

G  =  Gain* 

EIRP  =  Effective  isotropic  Radiated  Power 
^em  ~  Maximum  Effective  Aperture* 

AR  =  Axial  Ratio 

PLF  =  Polarization  Loss  Factor 

BW  =  Bandwidth 

fhigh  ~  Upper  Frequency  Limit 
f;;_Q;.,-  =  Lower  Frequency  Limit 
Acceptable  Conductor  Diameter 

Ex, y, z  =  Electric  Field  Cartesian  Components 
®p'^  ~  Unit  Polarization  Vector  Coordinate  Angles 

Og  =  Antenna  Unit  Polarization  Vector 
rmin  =  Minimum  Distance  to  the  Far-Field 
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ed  aata; 


D  =  Diameter  of  Helix  (Center  to  Center) 

S  =  Spacing  Between  Turns  (Center  to  Center) 

L  =  Length  Along  Conductor  of  One  Turn 
n  =  Number  of  Turns 
d  =  Diameter  of  Helical  Conductor 
f  =  Frequency  of  Interest 
m  =  Desired  Mode 

=  Antenna  Feed  Current 

X, y, ;  =  Antenna  Unit  Polarization  Vector  Cartesian  Coordinates 

1  =  Number  of  Increments  in  Degrees  for  Far  Field  Radiation  Pattern 
Zq  =  Characteristic  Feed  Impedance 

=  Incoming  Wave  Electric  Field  Unit  Vector 


Enter  input  data  here: 


ow-  I 


(dimensionless) 


D  )=  1074  (meters  ) 


L  :=  .766  (meters) 


n-=10  (turns) 


f  =  9,25  10*  (HZ) 


( dimensionless ) 


Zp:=150  (Q) 


S  •=  U766  (meters ) 


1  ■=  360  ( degrees  ) 


d;=  005  (meters) 


Io:=l  (A) 


(meters ) 


z:=1000  (meters) 


(meters ) 
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Calculate  helical  geometric  parameters  and  define  constants: 


c:=2  9979  10*  (meters  /  sec) 


C  =iD  (meters) 


C  =0  33741  (meters) 


:=  1 20  n  ( Q ) 


c 


1  =  0  3241 


(meters  /  cycle) 


(meters  /  cycle) 


Cl  =  1  04107 

a  —  atari 

a  =  0  22324 

,  180 
ad - a 


ad  =  1 2  7908 


( dimensionless ) 


( dimensionless ) 


( radians ) 


( radians ) 


( degrees ) 


(degrees ) 


/. 


SI  =0  23635 


A. 


LI  =  2  36349 


(dimensionless ) 


(dimensionless ) 


(dimensionless ) 


(dimensionless ) 


Calculate  helical  antenna  parameters: 


Define  angular  offset  9 
from  helical  axis: 


Minimum  Distance  to  the  Far-Field 


:  0.—  2  1  (radians  ) 


Tq  =161 


(meters ) 


a  (radians) 


r,  :=  5  n  S 


.2  (nSr 

^2 - 


(meters ) 


(meters ) 


rmin  :=  iTiax( r)  (meters) 


rnnn=3  83  (meters) 
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Relative  Phase  Velocity  p: 


P  ■= 


Li 


Si+  m  + 


(dimensionless) 


p  =  1  83736 


(dimensionless ) 


Array  Factor  Phase  Shift  % 


1*1(0)  =2n- 


|si  cos(6)- 


( radians ) 


Electric  Field  Field  Components  Eq, E^: 


(V  /  m) 


E0(9)  :=  E(0’l 


(V  /  m) 


E^(0):=j  E(0) 


(V  /  m) 


Radiation  Intensity  U (0)  : 


U(0)  :=1(|E(0)|)‘ 
% 


(W  /  solid  angle) 
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Radiated  Power  P 


Prad  :=  U(9)  sin(0)  dSd^  (W) 

■  0  0 


Prad  =  1  2394  10 


Directivity  D, 


Do  ■=  12  n  Sa.  (dimensionless)  Do2  •= 


4  n  U(0) 


(dimensionless ) 


Do  =  30  73917  (dimensionless)  Do2  =26.89463  (dimensionless) 


Axial  Ratio  AR ; 


AR  =  LX  sin( 


(dimensionless) 


AR  =0.76309 


(dimensionless) 


Effective  Isotropic  Radiated  Power  EIRP: 


EIRP  =PradDo 


EIRP2  Prad  Do2 


EIRP  =0  0381 


EIRP2  =0.03333 
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Polarization  Loss  Factor  PLF: 


9p  =  alani-! - ^1  (radians) 


♦p  =  atan  -  (radians) 


9p  =  141421*  10  '  (radians)  =07854  (radians) 


E\  :=  E9(  Op )  cos(  0p )  cos(  ♦p )  -  E4(0p)  sm(4p) 


(V/m) 


Ey  :=E9(0p)  cos(6p)  sm(^p)-*- E^Op)  cos(^p)  (V/m) 


Ez  :=  F0(  Op )  sin(  Op )  •  1 


(V/m) 


oa 


_ 1 _ 

V(|Ex|  )^^(lEyl)^+(|Ez| 


'Ex 

Ey 

iE^ 


(dimensionless ) 


cja 


05“  0.5j 
0.5  +  0  5j 

-9  99999  lO”** 


(dimensionless ) 


( dimensionless ) 


PLF  -  1  (dimensionless) 

Radiation  Resistance  R^: 

Pj  2  Prad  (Q) 

( |lo| 

Rr  =  2  47881  10“^ 
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Dual  Parameters 


Axial  reed 


Peripheral  Feed 


Incut  Resistance  R: 


Ra  :=  140 


(Q) 


Rp:=ii0 


Ra=  142846-10^  (Q) 


Rp  =  147012  10^ 


Voltage  Reflection  Coef f ecient T : 


ra.=. 


R^ 

Ra-*-Z„ 


(dimensionless  1 


rp:- 


Rp-  Zq 
Rp+Z. 


fa  =-0  02443  (dimensionless)  Fp  =-0  01006 


Reflection  Efficiency 


era 


:=  1  -  ( Ifa]  (dimensionless)  erp  :=  1  -  (|rp|)' 


era  =  0.9994  (dimensionless,  erp  =0  9999 


(Q) 


(Q) 


(dimensionless ) 


(dimensionless ) 


(dimensionless ; 


(dimensionless ) 
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Gain  G: 

Ga  •»  era  Do 

(dimensionless) 

Gp  ;=  erp  Do 

(dimensionless ) 

Gadb  •=  10  log(a'a  Do) 

(dB) 

Gpdb  :=  10  log(ap  Do) 

(dB) 

Ga  =30.72082 

(dimensionless) 

Gp  =30  73606 

(dimensionless ) 

Gadb  =  14  87433 

(dB) 

Gpdb  =  14.87648 

(dB) 

Ga2  •=  Era  Do2 

(dimensionless ) 

Gp2  '  ^op  Do2 

(dimensionless ) 

Gadb2  —  1 0  log(  era  Do2 ) 

(dB) 

Gpdb2  :  =  1 0  log(  np  Do2 )( dB ) 

Ga2  =  26  87858 

(dimensionless) 

Gp2  =26.89191 

(dimensionless ) 

Gadb2  =  14  29406 

(dB) 

Gpdb2  =  14  29622 

(dB) 

Maximum  Effective  A; 

)erture  A-,™: 

.  .  era  Do  ^ 

Aema  •= - PLF 

(m^ ) 

.  Erp  Do  ,  ^ 

Aemp  - - PLF 

(m2 ) 

4it 

4  n 

Aema  =0  25679 

(m^) 

Aemp  =0  25691 

(m2 ) 

.  -  .  £ra  Do2  r,i  r- 

Aema2  ■= - PLF 

(m2) 

Aemp2  PLF  (m2)  I 

4  71 

4n 

Aema2  =0.22467 

(m2) 

Aemp2  =  0  22478 

(m2 ) 
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Maximum  Effective  Height 


hema  =  2  Rr  . 


hemp  .2  Rr^ 


hema  =2  5988  10 


hemp  =  2  59944"  10 


hema2  =2  Rr^ 


hemp2  =2 


hema2  =2.43086  10 


hemp2  =2  43146  10  ^  (m) 


Bandwidth : 

fhigh  ^  (Hz)  flow  (Hz) 

C  C 

ftugh  =  1  02179-10®  (HZ)  flow  =7  10809"  10*  (Hz) 

BW  :  fhigh -flow  (Hz) 

BW  =  3  10979  10*  (HZ) 

Acceptable  Conductor  Diameter: 

dmin  ■=  005  /.  (m)  dmax:=.05  /.  (m) 

dmin  =  1  62049  10  '  (m)  dmax  =0  0162  (m) 
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THE  BEVERAGE  ANTENNA 
MATHCAD  SOFTWARE - BEVERAGE .MCD 

The  Beverage  antenna  is  a  single  wire  structure  parallel  to  the  ground. 

It  IS  terminated  with  a  load  matching  the  characteristic  impedance  of 
the  wire.  Because  there  is  little  or  no  reflected  energy  from  the 
antenna's  termination,  the  Beverage  antenna  does  not  develop  a 
significant  standing  wave.  Therefore,  it  is  known  as  a  traveling  wave 
antenna.  The  Beverage  antenna  is  also  known  as  a  slow  wave  antenna 
since  the  relative  phase  velocity  along  the  wire  is  usually  less  than 
one . 

The  Beverage  antenna  is  used  for  a  wide  range  of  frequencies,  depending 
in  Its  length  and  the  characteristics  of  the  ground  under  the  antenna. 

It  transmits  and  receives  vertically  polarired  electromagnetic  waves 
primarily  through  a  cone  shaped  main  beam  pointing  in  the  direction  of 
the  traveling  wave.  In  the  far-field,  the  electric  field  pattern  above 
the  ground  can  be  considered  rotationally  symmetric  with  respect  to  the 
axis  of  the  antenna. 

A  Beverage  antenna  exhibits  a  highly  directional  main  lobe  and  resistive 
input  impedance  for  frequencies  corresponding  to  the  following  lengths: 

1  (wavelengths) 

(Note:  1  in  a  subscript  indicates  the  dimension  is  in  wavelengths. 

Mathcad  equations  can  not  use  symbolic  subscripts.  Therefore,  the  symbol 
).  will  immediately  follow  the  parameter  in  equations  (i.e.,  Ll)  to 
indicate  the  dimension  is  in  wavelengths.) 

The  Beverage  antenna  Mathcad  application  will  compute  the  following 
parameters : 

=  Length  of  Antenna  in  Wavelengths 
/.  =  Wavelength 

Dq  =  Directivity 
p  =  Relative  Phase  Velocity 

E  =  Electric  Field  (No  Ground  Effects) 

Et  =  Electric  Field  (Total  Field  Including  Ground  Effects) 

u  =  Radiation  Intensity 

^rad  ~  Radiated  Power 

Zq  =  Antenna  Characteristic  Impedance 
r  =  Voltage  Reflection  Coefficient 

Py.  =  Vertical  Reflection  Coefficient  (Ground  Reflection) 

=  Reflection  Efficiency 
G  =  Gain 

EIRP  =  Effective  Isotropic  Radiated  Power 

Aem  =  Maximum  Effective  Aperture 

PLF  =  Polarization  Loss  Factor 

®max  “  Angle  of  Maximum  Radiation 

BW  =  Bandwidth 

-high  “  Upper  Frequency  Limit 
flow  “  Lower  Frequency  Limit 
Acceptable  Conductor  Diameter 


125 


■  =  Ground  Relative  Complex  Permittivity 
(Note;  The  subscript  '  will  be  annotated  as  p  in  the  application) 
dj^jax  ~  Maximum  Acceptable  Conductor  Diameter 
tmin  =  Minimum  Distance  to  Far-Field 
Rj.  =  Radiation  Resistance 
X  =  Electric  Field  Function 

The  following  data  must  be  input  based  on  known  or  estimated  data: 

h  =  Height  of  Antenna  above  ground 
L  =  Length  of  Antenna 
d  =  Diameter  of  Conductor 
f  =  Frequency  of  Interest 

1  =  Number  of  Increments  in  Degrees  for  Far  Field  Radiation  Pattern 
=  Load  Impedance 

=  Incident  Wave  Electric  Field  Unit  Vector 
63  =  Incident  Wave  Arrival  Angle 
(7  =  Ground  Conductivity 
Ej-  =  Ground  Relative  Permittivity 
Iff  =  Distance  of  Field  Calculations 
Iq  =  Input  Current  at  Antenna  Terminals 


h  ■=  1  (meters)  ^0  “  *  (amps) 

Z1  =  150  (Q)  L-=200  (meters) 

6a  ^  20  (radians)  f  :=  1  6  10^  (Hz) 

180 
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Calculate  Beverage  antenna  geometric  parameters  and  define  constants: 


c  :  =  2.9979  10*  (meters  /  sec) 


.  .  c 
f 


(meters  /  cycle) 


oa  •: 


(dimensionless ) 


1  =  187369*10  (meters  /  cycle) 


’lo 


:=  120  7i  (Q) 


LI  ■=- 


(dimensionless) 


1  .9 

£  := -  10  (Farads  /  m) 

°  36  n 


LI  =  1  06741 


(dimensionless) 


Calculate  Beverage  antenna  parameters: 


Define  angular  offset  6 
from  Beverage  antenna  axis: 


Minimum  Distance  to  Far-Field 


iminc. 


0:=O,-.,  Tt  (radians) 

i 


ro:M.6  1 


(m) 


r,  =  5L 


(m) 


♦  :  =  0,  M  ^  *  radians ) 


(m) 


rmin  =max(r)  (m) 


rmin  =  1  10"' 


(m) 


Relative  Phase  Velocity  p: 


(j.  I,  .038523»21 

- 1  (dimensionless) 

1000/ 


p  =0.87551 


( dimensionless ) 
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4 


1 


Maximum  length  Lm^y: 


Lmax :  =  • 


1 


4  l-i -  cos(0a) 

\P 


( wavelengths ) 


Lmax  =  1  23457 


(wavelengths ) 


Wavenumber  k: 


1 


k  =0  03353  (m"l) 


Electric  Field  Function  X 


X(9)  —  ( 1  -  cos(6))  (radians) 

2 


Electric  Field  Without  Ground  Effects  E (0)  ; 


30  k  L  L  sin(9) 

E(0);  = - P - 

rfif 


sin(X(9)) 

X(0) 


{W;to) 


Relative  Complex  Dielectric  Coefficient  Cr '  ; 


crp  :=Er- j 


2  n  f  £„ 


(dimensionless ) 


Erp=2-  1.125*10^j 


(dimensionless) 
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Vertical  Reflection  Coefficient  f'"; 


rv(0) 


f  n  JI 

ap  cos —  0  -  erp-  sin  -■ 


(dimensionless ) 


op  cos 


(i-oj  Jap-sin(!-0) 


Electric  Field  With  Ground  Effects  Et(0) 


Et(0)  .=  E(0)Al-  r\'(0)  e 


-j  2khco<^  — 


(V/m) 


Angle  of  Maximum  Radiation  0^^ 


0max  :=  acoa  1  - 


( radians ) 


0ma\  =0.86001 


(radians) 


1 80 

0ma\  —  =49  27486 


(degrees ) 


Radiation  Intensity  U(0) 


U(0)  ■=iH_  (|E(0)l )' 


(W  /  solid  angle) 


Umax  ■=  U(0max) 


(W  /  solid  angle) 


Umax  =  19  PI 34 
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(W  /  solid  angle) 


Radiated  Power  PrArt- 


Prad  =  U(0)  sin(0)  dBd* 

Jo  .'0 


Prad  =  32  66448 


Directivity  : 


4  71  Umax 


(dimensionless ) 


=  7  37542 


(dimensionless ) 


Polarization  Loss  Factor  PLF; 


PLF  =(|crwoa|) 


(dimensionless ) 


PLF  =  1 


(dimensionless  J 


Characteristic  Impedance; 


z„  =  (it) 


Z  =2  62626  lO" 
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Voltaqe  Reflection  CoeffecientF: 
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Maximum  Effective  Aperture  A^n,: 


■2  r. 

CTV'  1  D„ 


Aem :  ®  - 


PLF 


4  !l 


(m^) 


Aem  =  190698  10'‘ 


(m^ ) 


Radiation  Resistance  Rr^ 


Rr  ■= 


2  Prad 


(Q) 


Rr=  65  32897 


Bandwidth : 


fhigh  •-  — 
La 


fhigh  =561713  10* 


(Q) 


flow  := 


■  5  c 
U 


.8 


flow-  =  1.40428*  10 


BW  :  =  fhigh  -  flow 


BW  =4. 2 1285*1  o' 


8 


(HZ) 


(HZ) 


(Hz) 


(Hz) 


(Hz) 


(Hz) 
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Maximum  Effective  Height  h,,^: 


hem  =2 

^  % 


hem  =  1  14972  10^  (m) 


Acceptable  Ccr>ductsr  Diameter: 


dmax  :=  .0 1  Li  (m) 


dmax  =0.01067  (m) 
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BEVERAGE  ANTENNA  FAR-FIELO  RADIATION  PATTERN 


For  the  purpose  of  this  far-field  radiation  pattern,  the  Beverage 
antenna  axis  is  equivalent  to  the  Ey  =  0  grid  line  (Note;  Antenna  height 
above  ground  is  negible  in  the  far-field) .  The  pattern  is  essentially 
symetric  when  rotated  about  the  antenna's  axis  above  the  ground  plane. 


Electric  Field  With  Ground  Effects 


E\(6)  ■=El(9)  cos(6) 


Ey(0)  ■  =  Et(e)  sm(9) 


THE  SMALL  LOOP  ANTENNA 
MATHCAD  SOFTWARE-SMLOOP.MCD 

The  small  loop  antenna  is  a  coil  of  one  or  more  turns  whose  radius 
(a) satisfies  the  following: 

a  <  A./6n  (m) 

where  (1)  is  the  wavelength  of  the  frequency  of  interest.  Small  loops 
are  commonly  used  to  receive  signals  in  the  lower  frequency  regions. 

They  are  also  used  for  direction  finding  and  UHF  transmissions.  The 
efficiency  of  a  transmitting  loop  antenna  is  typically  very  low. 

However,  antenna  efficiency  can  be  improved  by  inserting  a  ferrite  core 
in  the  loop,  increasing  loop  perimeter,  or  increasing  the  number  of 
turns . 

The  small  loop  Mathcad  applications  will  analyze  three  geometries:  free 
space,  horizontal  loop,  and  vertical  loop.  Each  geometry  will  also 
examine  air  and  ferrite  cores.  Several  antenna  parameters,  particularly 
those  for  loops  in  free  space,  can  be  calculated  using  more  than  one 
formula.  Where  this  occurs,  multiple  results  will  be  computed  for 
comparison.  Computations  which  are  identical  for  all  geometries  will  not 
be  repeated.  Computations  with  ferrite  cores  will  use  effective 
permeability  (ng)in  lieu  of  free  space  permeability  ((Iq)  . 

(Note:  1  in  a  subscript  indicates  the  dimension  is  in  wavelengths. 

Mathcad  equations  can  not  use  symbolic  subscripts.  Therefore,  the  symbol 
A.  will  immediately  follow  the  parameter  in  equations  (i.e.,  Ll)  to 
indicate  the  dimension  is  in  wavelengths.) 

The  small  loop  antenna  Mathcad  application  will  compute  the  following 
parameters : 

S  =  Cross-Sectional  Area 

C  =  Circumference 

k  =  Wavenumber 
1  =  Wavelength 

Dq  =  Directivity 

E  =  Electric  Field  (No  Ground  Effects) 

E^  =  Electric  Field  (Total  Field  Including  Ground  Effects) 

U  =  Radiation  Intensity 

Uj^ax  ~  Maximum  Radiation  Intensity 

Pj-ajj  =  Radiated  Power 

r  =  Voltage  Reflection  Coefficient 

Fv  =  Vertical  Reflection  Coefficient  (Ground  Reflection) 

Fh  =  Horizontal  Reflection  Coefficient  (Ground  Reflection) 

E^v  ~  Reflection  Efficiency 
G  =  Gain 

EIRP  =  Effective  Isotropic  Radiated  Power 
^em  “  Maximum  Effective  Aperture 

PLF  =  Polarization  Loss  Factor 
Bandwidth 

£r'  =  Ground  Relative  Complex  Permittivity 

(Note:  The  subscript  '  will  be  annotated  as  p  in  the  application) 
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■■^rtun  ~  Minimum  Distance  to  Far-Field 
Rohmic  “  Ohmic  Resistance 
R^  =  Radiation  Resistance 

^s  =  Surface  Impedance  of  Conductor 
‘■cd  =  Conduction-Dielectric  Efficiency 
He  =  Effective  Ferrite  Core  Permeability 
=  Input  Reactance 
Zj_  =  Input  Impedance 
hem  “  Maximum  Effective  Height 
CR  =  Core  Length  to  Diameter  Ratio 

The  following  must  be  input  based  on  known  or  estimated  data: 

h  =  Height  of  Antenna  above  ground 

a  =  Radius  of  Antenna 

b  =  Radius  of  Conductor 

N  =  Number  of  Turns 

0(~  =  Conductivity  of  Loop 

Rp/Ro  =  RR  =  Ohmic  Resistance  from  Proximity  to  Ohmic  Skin  Effect  Ratio 
Ddemag  “  Demagnetization  Factor 
Hf  =  Permeability  of  Ferrite  Core 
f  =  Frequency  of  Interest 

1  =  Number  of  Increments  in  Degrees  for  Far-Field  Radiation  Pattern 

0^  =  Incident  Wave  Electric  Field  Unit  Vector 

a  =  Ground  Conductivity 

Ej.  =  Ground  Relative  Permittivity 

tff  =  Distance  of  Field  Calculations 

Iq  =  Input  Current  at  Antenna  Terminals 

q  =  Loop  Spacing 
cl  =  Core  Length 
cd  =  Core  Diameter 
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lA 

1 

Enter  input  data  here; 

aw  '-z 

,01 

(dimensionless ) 

0  :  =  6  10’ 

(mhos  /  m) 

io) 

b  :=  .01 

(meters ) 

a  -  1 

(dimensionless ) 

1  :=  360 

(degrees ) 

rff:=  1  10* 

(meters ) 

h;=2.5 

(meters) 

II 

(dimensionless ) 

oc  :=  5.8  lo"' 

(mhos  /  m) 

cl  ;=  1 

(meters ) 

q  :=  03 


2  b 


(meters) 


(dimensionless) 


nf  :  =  4000  (4  n  10‘^  (Henrys  /  m) 


f:=3  10’  (HZ) 


cd  :*  05  (meters) 

Zp  =50  (Q) 


Ddemag  :=  4  10"^ 


(dimensionless) 


N  :  =  6 


(dimensionless ) 


Ijj  :=  1  (amps) 


a  :=  15  (meters ) 
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Calculate  small  loop  antenna  geometric  parameters  and  define  constants; 


I 


c  =2  9979  10  (meters  /  sec) 


=  120  n 


(meters  /  cycle) 


'  irr9 

t  = - 10- 

36  fl 


[Farads  /  m) 


1=9  993  (meters  /  cycle) 


Hij:  =  4iilCr  (H  /  m) 


S  :=  n  a 


D  :=2  a 


(meters ) 


S  =0.07069 


D  =0  3 


(meters) 


—  =0  53015 
6  1 


(meters ) 


CR  — 


(dimensionless ) 


C  :=2  Jia 


(meters ) 


oa  ’.=  0 


( dimensionless ) 


C  =0.94248  (meters) 


Calculate  small  loop  antenna  parameters  in  tvmm  spac*  (air  cora) 


Define  angular  offset  9 
from  small  loop  axis; 


Distance  to  Far-Field  rmin; 


9;=0, —  ..2n  (radians) 


♦  0,  ,.  2  n  (radians) 


Wavenumber  k: 


(m  ^) 


ro:=161 


r,:=5.D 


rmin  :=  max(  r) 


rmin  =  15  9888  (m) 


(m  1) 


k  =0  62876 


139 


Electric  Field  Without  Ground  Effects  E(8)  : 


kfS^„  sm(0) 


E{0).: - - e 

2  rff 


(V/m) 


Radiation  Intensity  U(0) 


U(0)  :=iH_  (lE(0)lr 
2% 


(W  /  solid  angle) 


Umax  ■: 


‘-(I 


(W  /  solid  angle) 


Umax  =9.30838  10 


(W  /  solid  angle) 


Radiated  Power  Prart^ 


Prad  =ti^_(ka)"(|g)‘ 


Prad=  7.80909  10 


Prad2  ;=  U{0)  sin(9)  d0d^ 

•  0  0 


Prad2  =7  79817  10 


Directivity  : 


Do.=  1.5 


(dimensionless ) 


Do2  •= 


4  It  Umax 


(dimensionless ) 


Do2  =  1.5 


(dimensionless ) 
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Radiation  Resistance  (Rr)  : 


141 


142 


Aem  =  11  91992 


Ac:m2  := 


ETV  8Cd  A.'  Do 


Aem2  =0  09427 


Effective  Isotropic  Radiated  Power  (EIRP) 


EIRP  :=  Prad  Do 


EIRP=0.01171 


E1RP2  :=  Prad2  Do 


EIRP2  =0  0117 


Maximum  Effective  Height  (h. 


hem  :=  2 


hem  =  0  26667 


hem2  :  = 


,  /Rt  Aem2  ^ 


hem2  =0.02371 


Bandwidth ; 


Polarization  Loss  Factor  PLF: 


Bandwidth  =- 


PLF  ;=  (|(jw  (ja|) 


(dimensionless ) 


Bandwidth  =  1  06029  10’ 


PLF  =  1 


(dimensionless ) 


SMALL  LCX3P  ANTENNA  FAR  FIELD  PATTERN  IN  FREE  SPACE  WITH  AIR  CORE 


For  the  purpose  of  this  far-fieid  radiation  pattern,  the  small  loop 
antenna  axis  is  equivalent  to  the  Ey  =  0  grid  line.  The  pattern  is 
symetric  when  rotated  in  the  (^)  dirction  about  the  antenna's  axis. 


Electric  Field  Without  Ground  Effects 


Ex(9)  :  i*-v6)|  cos(0) 


Ey(6)  =  |E(9)|  sin(9) 
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Calculate  small  loo 


Effective  Permeability 


1  4  Ddemag  (.  • 


(Henrys/m) 


ne  =  5  046t)3'  10 


(Henrys/m) 


Electric  Field  Without  Ground  Effects  E (6)  : 


kfSgel  sin(0) 

E(0)  ;= - 2 - e  ’ 

2  rflf 


(V/m) 


Radiation  Intensity  U(0) 


U(0)  :=iH_  (|E(0)|r 


(W  /  solid  angle) 


Umax  :=  U 1-^ 


(W  /  solid  angle) 


Umax  =  1.50127  10'* 


(W  /  solid  angle) 


Radiated  Power  P, 


Prad:=  U(0)  sin(0) d0d^ 

■  0  0 


Prad  =  1  2577-10^ 


Directivity  D, 


Do  .=  l  5 


(dimensionless ) 


Do2  := 


4  It  Umax 


(dimensionless ) 


Do2  =  1.5 


(dimensionless) 


Radiation  Resistance  (Rr) : 


Rr=2a/f-l 


Rr  =  9.06811  •10° 


Rr2:=^^N^ 

(KW 


Rr2  =9.05542  10 


Surface  Impedance  of  Conductor  R,;: 


Rs  :=  ^ 


Rs  =  0  09056 


;hmic  Resistance  Rp hmi ci. 


„  .  N  a  Rs  I  ^ 

Rohmic  = - (RR4-  1 ) 

b 


Rohmic  =9  37266 


Conduction-Dielectric  Ef  f iciencylCf-H)  : 


(dimensionless ) 


Rohmic  +  Rr 


ecd  =  1 


(dimensionless ) 


Input  Reactance  (Xt ) : 


Xm  :  2  xf  a  ne  |ln|8 -^1  -  1  75 


Xin  =4  3342  lO" 
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Input  Resistance  (R^ ) 


Rin  :»  Rr+  Rohmic 


Rui  =9  06812  10° 


Input  Impedance  (Zt 1 : 


Zi  :=  Rin-h  Xin  j 


Zi  =9  06812  10^  +4  3342  lo’j 


Voltage  Reflection  Coef fecient F ; 


Zi  +  z„ 


(dimensionless ) 


r  =0.99999  +5.2587  10  j 


(dimensionless ) 


Reflection  Efficiency  Ck.,: 


erv  :=l  -  (IH  T 


(dimensionless ) 


an-  =2.20048*10 


(dimensionless ) 


Gam  Gi 


G  ;=£rv  ocd  Do 


(dimensionless) 


Gdb  =  101og(a^'  ecd  Do) 


G  =3  30071  10 


(dimensionless ) 


Gdb  =-44  81392 


Maximum  Effective  Aperture  (Aq^) : 


trv  Gcd  Do 


Aem=2  62295  10 
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Effective  Isotropic  Radiated  Power  (EIRP) : 


EIRP  sPrad Do 

EIRP  =  188655  lo’ 

Bandwidth  BW : 

c 

Bandwidth  - 

6  n  a 

Bandwidth  =  1  06029  10* 

Maximum  Effective  Height  (h^n,)  : 

hem  •  ISL^  2 

•J  % 

hem  =  5  02363 

Polarization  Loss  Factor  PLF: 


PLF  =  1 


(W) 


(W) 


(Hz) 


(Hz) 


(m) 


(m) 


(dimensionless ) 

(dimensionless ) 
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E  FIELD 
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cor«) : 

Define  angular  offset  6 
from  small  loop  axis; 


2  2  I  2 


( radians ) 


( radians ) 


Relative  Complex  Permittivity  (g^ t 


op  :  =  cr-  j 


(dimensionless) 


2  n  f  E„ 


Erp  =-3.6"  10  j 


(dimensionless) 


Horizontal  Reflection  Coefficient  (Fh) 


m  :«  0  - 
•> 


( increments ) 


cos(6)-  sin(9) 

cos(Q)-t-  Jcrp-  sin(0)^ 


(dimensionless ) 


Fhl  := 

m 


-ifh 


/mnV 

)-  smi— I 


Jap  -  sin 


(dimensionless ) 


Total  Electric  Field  (Ethor) : 


Ethor(  0) 


k  S  f  Hg  Ip  sin(0)  e"-'  ‘‘ 


rh(0)  e' 


-j  2khcos(9)] 


(V/m) 


Ethor  1 

m 


~\  2  rff  / 


-J  2khl2!Lnil 


I  -  Fhl  e 


(V/m) 
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Voltage  Reflection  Coef fecient f : 


'zi  +  z„ 


(dimensionless) 


r  =0  64261  -*-0  75809] 


(dimensionless) 


Reflection  Efficiency  e^,,: 


£n  =1  -  (|r|) 


(dimensionless ) 


en-  =0  01236 


(dimensionless ) 


Gain  G I 


G  :  =  En  ocd  Do 


(dimensionless ) 


Gdb  :=  1 0  log( crv'  Ecd  Do) 


G  =  0  04728 


(dimensionless ) 


Gdb  =-13  25293 


Maximum  Effective  Aperture  (Apn,)  : 


.  erv  Bcd  1  Do  (PLF) 
Aem  •= - 


Aem  =0.37574 


Effective  Isotropic  Radiated  Power  (EIRP) : 


EIRP  :  =  Prad  Do 


EIRP  =0  04675 


Bandwidth : 


Maximum  Effective  Height  (h 


Bandwidth  - 


.  /Rr  Aem  ,  ,  ^ 

hem  = - 2  (m) 


Bandwidth  =  1  06029*  10*  (Hz) 


hem  =0.05401 
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3MIILI.  HORIZONTAL  LOOP  ANTINHA  FAR  FlgLD  PATTERN  ABOVE  GROUND  WITH  AIR 
CORZ 


cor#) : 

Incremental  Horizontal  Reflection  Coefficient  (Thl): 
m  =0  —  (increments) 


rhl 

m 


|m  nj 

(f)* 


op-  sm 

^m  n\ 

>  i 

/m  JiV 

)-  sinl - 1 


(dimensionless) 


Total  Electric  Field  (Ethor) : 


Ethor(0)  = 


k  S  f  ne  sin(  6)  e‘ 


rh(0)  e 


,-j  2khco*<C 


(V/m) 


Ethor  1 


ts  fue  e'J 

Tii  j 


.j 


1  -  rhl  e 


(V/m) 


Radiation  Intensity  U(0): 


U(0):  =  -!Sl  (|Ethor(e)|)^ 


(W  /  solid  angle) 


Umax  :  = - ( |ma\(  Ethor  1  )| 

Umax  =5.9995  lO"* 
Radiated  Power  P^aH- 


(W  /  solid  angle) 


(W  /  solid  angle) 


2  It  - 

2 

Prad  =  lU(0)  sin(0)|  d0d4  (W) 

Jo  ^0 


Prad  =  1  63918-10^ 


Directivity  D, 


4  n  Umax 


( dimensionless ) 


Do  =4.59936 


(dimensionless) 
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Radiation  Resistance  (Rr)  : 


(Q) 

Rr=l  18021  lO’ 

(Q) 

Surface  Impedance  of  Conductor  R«: 

Rs.=  S 

(Q) 

ac 

Rs  =  0  09056 

(Q) 

Ohmic  Resistance  R^Hrm,-: 

Rohmic  (RR+  1 ) 

b 

(Q) 

Rohmic  =  9  37266 

(Q) 

Conduction-Dielectric  Ef ficiencyjc^H) : 


Dcd 


Rr 


Rohmic  4  Rr 


£Cd  = 


Input  Reactance  (X^ ) : 


Xin  •:  2  nf  a  ^e  |ln|8  — I  -  1  75 
bj 


(dimensionless ) 


(dimensionless) 


5 


Xin=4  3342  10 


Input  Resistance  (Rt ) : 


Rin  ■=Rr+  Rohiruc 


Rin  =  1  18021  10 


(Q) 


(Q) 


(Q) 


(Q) 


Voltage  Reflection  Coef fecientj: 


r 


Zi+  z„ 


r  =0  64261  +0  75809J 


(dimensionless ) 


(dimensionless) 
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Reflection  Efficiency  t 


cn.  =  1  -  ( |r| 


(dimensionless) 


ETV'  =0.01236 

Gain  G: 

G  :  =  Erv  acd  Do 

Gdb  :=  10  log(£rv  ocd  Do) 

G  =  0  05684 

Gdb  =-12  45349 


(dimensionless) 


(dimensionless ) 
(dB) 

(dimensionless) 

(dB) 


Maximum  Effective  Aperture  (Ap^)  ^ 

.  .  erv  ocd  1^  Do  ( PLF )  2 

Aem  ■= - - -  (m"^) 

4  n 

Aem  =045168  (m^) 


Effective  Isotropic  Radiated 

Power  (EIRP) : 

EIRP  Prad  Do 

(W) 

EIRP=7  5392  lo’ 

(W) 

Bandwidth  BW: 

BW  :-i— 

(Hz) 

6  n  a 

BW  =  1  06029  10* 

(Hz) 

Maximum  Effective  Height 

(  h  OTTI  )  ' 

_j _ em, 

Polarization  Loss 

Factor  ' PLF) : 

,  iRr  Aem  . 

hem  :=  |  2 

(m) 

PLF  :=U<nv  <ja|J 

(dimensionless 

-J  % 

hem  =  2  37827*  10^ 

(m) 

PLF  =  1 

(dimensionless 
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SMALL  HORIZONTAL  LOOP  ANTBHHA  FAR  FIELD  PATTERN  ABOVE  GROUND  WITH 
IXRRZTE  CORE 


For  the  purpose  of  this  far-field  radiation  pattern,  the  small  loop 
antenna  axis  is  equivalent  to  the  Ex  =  0  grid  line.  The  pattern  is 
symetric  when  rotated  in  the  (^)  dirction  about  the  antenna's  axis. 


Small  Horizontal  Loop  with  Ferrite  Core 


Electric  Over  Real  Ground 


Ex{0)  :=  lEthor(0)| 


Ey{0)  :=  |Ethor(0)| 


Define  angular  offset  9 
from  small  loop  axis: 


e:  =  0,-  n 

1 


( radians ) 


♦  :  =  0.-  !T 
1 


( radians ) 


Relative  Complex  Permittivity  (c^ ■ ) : 


op  •=er-  j 


2  f  c. 


0 


(dimensionless ) 


□p  =-3  6  lO'^'j 
Vertical  Reflection  Coefficient  (Tv) : 


( dimensionless ) 


rvte)  := 


m  :  =  0  1 


(dimensionless ) 


( increments ) 


fvl  := 


ffn  mm 

l-j 

op-  sinj 

^  m  n\  ‘ 

U'  J 

//  n  mm 

'^“UrTi 

)\ 

erp-  smj 

^  n  m  ' 
^2  1  1 

(dimensionless ) 


Total  Electric  Field  (Etvert) 


Etvert(6)  ■=  j 


fk  S  f  sin(0)  e'-’ 


2  rff 


(l  +  rv(6)  eJ 


(V/m) 


Etvcrl  := 


|k  S  f4pIp  sm| 

1  1  J 

1  e'J  ''  A 

\  2rtr  / 

•  J  2  k-h  cosl^— ] 


+  Pvl  e 

m 


(V/m) 
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Radiation  Intensity  U(9): 

U(0):  =— .(|Etven(0)|r 

Umax  :  = - ( |max(  Etverl )  | 

"’V) 

Umax  =  3  72332  lO"' 

Radiated  Power  PraH- 


Prad  :  = 


|U(0)  sin(0)|  d0d4 


Prad  =001016 


directivity  D^: 


Do:  = 


4  n  Umax 
Prad 


Do  =4  60367 


Radiation  Resistance  (Rr) 


Rr:=il2l  N= 


Rr=073176 


Surface  Impedance  of  Conductor  R^ 


oc 

Rs  =  1  42898- lO"' 


(W  /  solid  angle) 


(W  /  solid  angle) 

(W  /  solid  angle) 


(W) 


(W) 


(dimensionless ) 


(dimensionless ) 


(Q) 


(Q) 


(Q) 


(Q) 
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Gain  G: 


G  =  trv  ocd  Do 


(dimensionless ) 


Gdb  ■*  10  log(tr\  ccd  Do) 


G  =0  04733 


(dimensionless ) 


Gdb  =-13  24888 


Maximum  Effective  Aperture  (A^n,)  : 


.  ETv  ocd  1  Do  (PLF) 
Aem  ■= - 


Aem  =0  37609 


Effective  Isotropic  Radiated  Power  (EIRP) 


EIRP  sPrad  Do 


EIRP  =0.04679 


Bandwidth : 


Bandwidth  =- 


Bandwidth  =  1  06029  10* 


Maximum  Effective  Height  (hpm) = 


hem  :=  2 


hem  =  0  05404 


Polarization  Loss  Factor  (PLF) 


oa  1 


(dimensionless ) 


ow  :  =  1 


(dimensionless ) 


:=  (|(TW  oa|) 


(dimensionless ) 


PLF  =  1 


(dimensionless ) 
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Total  Electric  Field  (Etvert) : 


Etvt*rt(01  ■ 


=  fc 


kS  fue  L  sin(e)  e-' 


2  rff 


(V 


Etverl  :  = 


^k  S  fue  I.  sinfcVe'J 


2  rff 


■  j  2k  h  c«|l!!-:| 


Radiation  Intensity  U (9)  : 


U(0) 


:  = - ( |Etvert{0)| ) 


rff^ 


Umax  :  =  (  |max( Etverl  )| 


Umax  =6-00502  10 
Radiated  Power  Prari- 


Prad  •= 


|U(0)  sm(0)|  d0c4 


Prad  =  1  63918M0' 


il+fvl  e 

m 


(W  /  solid  angle) 


(W  /  solid  angle) 


(W  /  solid  angle) 


(W) 


(W) 


Directiv:.ty  0^: 


_  4  n  Umax 

Do  = - 

Prad 

Do  =4  60361 

Radiation  Resistance  (Rr)  : 

R,  n' 

(I'ol)' 

Rr=  I  I80:i  lo’ 

Surface  Impedance  of  Conductor  R«;: 

Rs  = 

Rs  =0  09056 

Ohmic  Resistance  R^hmi r : 

Rohmic  :s  ^(RR.».  i ) 
b 

Rohmic  =  9  37266 


(dimensionless) 


(dimensionless ) 


(Q) 


(Q) 


(Q) 

(Q) 


(Ql 

(Q) 


Conduction-Dielectric  Ef f iciencyjgf-rt >  = 

Rr 

Ecd  := -  (dimensionless) 

Rohmic  -f  Rr 

£cd  =  1  (dimensionless) 

Input  Reactance  (Xt ) : 

Xin  =2nfane|ln|8lj-175j 

Xin=4  3342  10’  (U) 

Input  Resistance  (R^ ) : 


Rm 

:=  Rr-t-  Rohmic 

(Q) 

Rin 

=  1  18021*10'' 

(Q) 
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Voltage  Reflection  Coef fecient T: 


Zl4Z„ 

r  =0  64261  +0  75809J 

Reflection  Efficiency  g, 

=  1  -  ( Ir! 

Er\’  =  0  01236 


(dimensionless ) 


(dimensionless ) 


(dimensionless ) 


(dimensionless ) 


Ga in  G : 

G  :  =  Em  ecd  Do 

Gdb  '.=  10  log((T\  ecd  Do) 

G  =  0  05689 

Gdb  =-12  44948 


Maximum  Effective  Aperture  (A^n,)  : 


_  ETv  ocd /,  Do  (PLF) 


Aem  =0  4521 


Effective  Isotropic  Radiated  Power 
(EIRP) : 


EIRP  :=  Prad  Do 


EIRP  =7  54613  lO’ 


Polarization  Loss  Factor  (PLF): 


(dimensionless ) 

(dB) 

(dimensionless ) 

(dB) 

Bandwidth : 

„  .j.  c  (H2> 

Bandwidth  := - 

6  n  a 

Bandwidth  =  1  06029- 10* 

Maximum  Effective  Height  (h^n^i  : 


,  [Rr  Aem  * 
hem  •=  f - ——  2 


hem  =  2  37936- 10“ 


0\ 

fo^ 

(ja  :  = 

1  (dimensionless) 

aw  := 

1 

0, 

lo, 

:=(U»|)^ 


(dimensionless ) 


(dimensionless ) 


PLF  =  1 
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(dimensionless ) 


CORE 


For  the  purpose  of  this  far-field  radiation  pattern,  the  small  loop 
antenna  axis  is  equivalent  to  the  Ey  =  0  grid  line.  The  pattern  is 
symetric  when  rotated  in  the  (^)  dirction  about  the  antenna's  axis. 

Small  Vertical  Loop  with  Ferrite  Core 
Electric  Over  Real  Ground 


Ex(0) |Etvert(0)|  cos(0)  Ey(0)  :=  |Etvert(0)|  sin(0) 


—  EFELD 
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THE  LARGE  LOOP  ANTENNA 
MATHCAD  SOFTWARE- LGLOOP. MOD 


The  large  loop  antenna  is  a  coil  of  one  or  more  turns  whose  radius  (a) 
satisfies  the  following: 


a  >  i/6n  (m) 

where  (1)  is  the  wavelength  of  the  frequency  of  interest.  Large  loops 
are  not  commonly  used  and  are  considered  impractical  if  the  radius 
exceeds  one  wavelength.  Like  the  small  loop,  the  efficiency  of  a 
transmitting  large  loop  can  be  low.  However,  antenna  efficiency  can  be 
improved  by  inserting  a  ferrite  core  in  the  loop  or  increasing  the 
number  of  turns. 

The  large  loop  Mathcad  applications  will  analyze  two  geometries:  free 
space  and  the  horizontal  loop.  Each  application  for  a  specific  geometry 
will  examine  air  and  farrita  cores.  Several  antenna  parameters, 
particularly  those  for  loops  in  free  space,  can  be  calculated  using  more 
than  one  formula.  Where  this  occurs,  multiple  results  will  be  computed 
for  comparison.  Computations  which  are  identical  for  all  geometries  will 
not  be  repeated.  Computations  with  ferrite  cores  will  use  effective 
permeability  (Hgiin  lieu  of  free  space  permeability (Hq) . 

(Note:  1  in  a  subscript  indicates  the  dimension  is  in  wavelengths. 

Mathcad  equations  cannot  use  symbolic  subscripts.  Therefore,  the  symbols 
will  immediately  follow  the  parameters  in  equations  (i.e.,  Li)  in  lieu  of 
subscripts . ) 

The  large  loop  antenna  Mathcad  application  will  compute  the  following 
parameters : 

S  =  Cross-Sectional  Area 
C  =  Circumference 
k  =  Wavenumber 
i  =  Wavelength 
Dq  =  Directivity 

E  =  Electric  Field  (No  Ground  Effects) 

E^  =  Electric  Field  (Total  Field  Including  Ground  Effects) 

U  =  Radiation  Intensity 

Ufjjax  “  Maximum  Radiation  Intensity 

°rad  "  Radiated  Power 

=  Vertical  Reflection  Coefficient  (Ground  Reflection) 

=  Horizontal  Reflection  Coefficient  (Ground  Reflection) 

G  =  Gain 

EIRP  =  Effective  Isotropic  Radiated  Power 
Aejn  =  Maximum  Effective  Aperture 
Bandwidth 

Ej-i  =  Ground  Relative  Complex  Permittivity 

(Note:  The  subscript  '  will  be  annotated  as  p  in  the  application) 

^min  ~  Minimum  Distance  to  Far-Field 
Rohmic  ~  Ohmic  Resistance 
Rj.  =  Radiation  Resistance 

^s  =  Surface  Impedance  of  Conductor 
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Ecd  =  Conduction-Dielectric  Efficiency 
Ug  =  Effective  Ferrite  Core  Permeability 
^em  “  Maximum  Effective  Height 
CR  =  Core  Length  to  Diameter  Ratio 

The  following  data  must  be  input  based  on  known  or  estimated  data: 

h  =  Height  of  Antenna  above  ground 

a  =  Radius  of  Antenna 

b  =  Radius  of  Conductor 

N  =  Number  of  Turns 

Oc  =  Conductivity  of  Loop 

Rp/R©  =  RR  =  Ohmic  Resistance  from  Proximity  to  Ohmic  Skin  Effect  Ratio 
^demag  ~  Demagnetization  Factor 
fif  =  Permeability  of  Ferrite  Core 
f  =  Frequency  of  Interest 

1  =  Number  of  Increments  in  Degrees  for  Far  Field  Radiation  Pattern 

o  =  Ground  Conductivity 

Ej.  =  Ground  Relative  Permittivity 

rjf  =  Distance  of  Field  Calculations 

I©  =  Input  Current  at  Antenna  Terminals 

q  =  Loop  Spacing 

cl  =  Core  Length 

cd  =  Core  Diameter 

Z©  =  Characteristic  Impedance 


Enter  input  data  here: 


a  :  = 

46 

(meters ) 

£r:=  1 

(dimensionless 

b  := 

05 

(meters ) 

rff:=  1  10’ 

(meters ) 

1  ■= 

360 

(degrees ) 

RR  =  15 

(dimensionless 

h:  = 

72 

(meters ) 

cl  ■=  1 

(meters ) 

oc  - 

=  58  10^ 

(mhos  /  m) 

Zo  :=  50 

(Q) 

q  .= 

03 

(meters ) 

Ddemag  '=  4  10'^ 

(dimensionless 

uf:: 

=  4000  (4  )T  10'^ 

(Henrys  /  m) 

( amps ) 

0  ■  = 

6  lO’ 

(mhos  /  m) 

f  ■=  3  26  10* 

(Hz) 

q 

2  b 

=  03 

( dimensionless ) 

cd  ■=  .05 

(meters ) 

N  :=  1 

(dimensionless 

169 


Calculate  large  loop  antenna  geometric  parameters  and  define  constants; 


c:=2  9979  10*  (meters  /  sec] 


:=  120  n 


(meters  /  cycle)  c  - - 


36  n 


(Farads  /  m) 


1=0  91%  (meters  /  cycle)  =4nl0' 


(H  /  m) 


S  :=  a* 


D  ;=2  a 


(meters ) 


S  =0  66476  (m*-) 


D  =0  92 


(meters ) 


—  =004879  (meters) 

6  Tl 


CR  — 


(dimensionless ) 


C  =  2  71  a  (meters ) 


C  =2  89027  (meters) 


Calculate  large  loop  antenna  parameters  in  free  space  (air  core) 


Define  angular  offset 
from  small  loop  axis; 


Distance  to  Far-Field  tn 


2  71 

0::O,—  ,  2  71 


( radians ) 


To  .  =  161 


r,  D 


Wavenumber  )<: 


( radians ) 


rmin  ~  max(  r) 


rmin  =  4  6 


k  =  6  83251 
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Directivity 


Dos  682 


Do  =2  1435 


Do2  := 


4  n  Umax 
Prad 


Do2  = 1  79871 


Radiation  Resistance  (Rr)  : 


Rr:s60  /£ 
A. 


Rj  =  1  86118'10-’ 


Rr2  •  = 


2  Prad 


2 


Rr2  =2  19898  lO’ 


Surface  Impedance  of  Conductor  R^.: 


Rs  =4  71058  10 


(dimensionless ) 


(dimensionless ) 


(dimensionless ) 


(dimensionless ) 


(Q) 


(Q) 


(fi) 


(Q) 


(Q) 


(Q) 
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Ohmic  Resistance 


Rohmic  : »  ^  ^  ^  ( RR  +  1 ) 
b 


Rohmjc  =  0  04984 


Conduction-Dielectric  E f f iciencyjc^H  )  ^ 


Rohmic  -f  Rr 


Ecd  =  0  99997 


G  3  i,  n  G  J 


G  ;=  Ecd  Do 


Gdb  :=  1 0  log(  ocd  Do ) 


G  =  2  14344 


Gdb  =  3  31111 


Maximuir.  Effective  Aperture  iAgj^i  : 

Aem  =  0543  C  ) 


Aem  =  0  14432 


^  ,  Ecd  Do 

Aem2  := - 

4  :i 


Aem2  =  0  1 4424 


(Q) 


(Q) 


(dimensionless ) 


(dimensionless ) 


(dimensionless ) 


(dB) 


(dimensionless ) 


(dB) 


Im^ ) 


(m^) 

(m*^ ) 


(m^ ) 
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Effective  Isotropic  Radiated  Power  ( 


EIRP  ;  =  Prad  Do 


EIRP=2  3S675  lO'’ 


Bandwidth : 


Bandwidth  :*  - 
a 


Bandwidth  =617143  10* 


Maximum  Effective  Height  (h^^) : 


he„  :=  Islase  , 
-I  "Vi 


hem  =  1  68821 


hem2 


JRr  Aem2  ^ 

% 


hem2  =  1  68775 


IRP)  ; 


(W) 


(W) 


(Hz) 


(Hz) 


(m) 


(m) 


(m) 


(m) 
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LARGE  LOQg  ANTENNA  FAR  FIELD  PATTERM  IN  FREE  3PACE  WITH  AIR  CORE 


For  Che  purpose  of  this  far-field  radiation  pattern,  the  large  loop 
antenna  axis  is  equivalent  to  the  Ey  =  0  grid  line.  The  pattern  is 
symetric  when  rotated  in  the  {^)  dirction  about  the  antenna's  axis. 

Electric  Field  Without  Ground  Effects 

Ext0)  :=  |E(e)|  cos(0)  Ey(0)  ;=|E(0)|  sin(0) 


•0  00.5  -0  002  -0  002 

-0  001 

-0  001 

0 

0  001 

0.001 

0.002 

0C02  000? 

-  E  FIELD 

Ex(  0) 

Calculace  large  loop  antenna  parameters  in  fr»m  sp*J«  (fTrif  cor*) 


Effective  Permeability  (u^): 


ne 


1  Ddemag  (nf-  1 ) 


(H/m) 


4e=5  04663  lO*' 


(H/m) 


Electric  Field  Without  Ground  Effects  E(9): 


E(0)  : 


2  rff  I 


g-(j  a  sin(0)) 


(V/m) 


ol  •  (  J  k  rff  )  r  I  I  I 

-I  ^  J1  k  a  sin 


(2  ;(  f  a  (ne  I 

- L|  e 

2  rfi' 


Radiation  Intensity  U(0): 


H) 


(V/m) 


U(0)  ;=_  (lE(0)|  )^ 


>J1„  (|ElJ)  = 


Umax  ;s  ma\(  U 1 ) 


(W  /  solid  angle) 


(W  /  solid  angle) 


(W  /  solid  angle) 


Umax  =  2  5382'  10"^ 


(W  /  solid  angle) 
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Radiated  Power  P 


Prad  :»  U(0)  sin(9)  dSd^ 

■  0  0 


Prad  =  I  77327-10' 


Directivity  D, 


Do:=  682- 


(dimensionless ) 


Do  =2  1435 


(dimensionless) 


Do2  .= 


4  1  Umax 


(dimensionless ) 


Do2  =  1  7987 1 


(dimensionless ) 


Radiation  Resistance  (Rr) 


Rr=rT‘60-ra 

^  Vo 


Rr  =  3  00174-10' 


Rr2.  =  iIl!iN- 

w 


Rr2  =3  54654  lo' 


Rs  •» 


Rs  =0  29852 


Ohmic  Resistance  Ri^hm-;  r  • 


Rohmic  ( RR4.  1 ) 

b 


Rohmic  =  3  1 5832 


Conduction-Dielectric  Ef ficiencyjc^^ )  : 


ecd  := 


Rr 

Rohmic  +  Rr 


Gcd  =  1 


Gam  G: 


G  :=  ecd  Do 


Gdb  =  10  log(ocd  Do) 


G  =  2  1435 


Gdb  =  3  31 123 


(Q) 


(Q) 


(Q) 


(G) 


(dimensionless ) 


(dimensionless ) 


( dimensionless ) 


(dB) 


(dimensionless ) 


(d3) 


1^8 


Aem  =0  14425 


) 


Effective  Isotropic  Radiated  Power 


EIRP  :=  Prad  Do 


EIRP=3  801  10‘“ 


Bandwidth : 


Bandwidth 


Bandwidth  =6  17143  10* 


Maximum  Effective  Height  (h^^) : 


hem  = 


(rt  Aem 

J  Ho 


hem  =6  77808- 10-' 


(EIRP) : 


(W) 


(W) 


(Hz) 


(Hz) 


(m) 


(m) 


lARGE  LOOP  ANTtNKA  FAR  FIELD  PATTERN  IN  FREE  SPACE  WITH  FERRITE  CORE 


For  the  purpose  of  this  far-field  radiation  pattern,  the  large  loop 
antenna  axis  is  equivalent  to  the  Ey  =  0  grid  line.  The  pattern  is 
symetric  when  rotated  in  the  (4)  dirction  about  the  antenna's  axis. 


Electric  Field  Without  Ground  Effects 

E\{0)  •=  1E(0)|  cos(6)  Ey(0)  •=  |E(0)|  sin(0) 


—  E  FIELD 


Calculate  large  horizontal  loop  antenna  parameters  ovr  ground  (*ir 
cor*)  : 


Define  angular  offset 
from  small  loop  axis: 


( radians ) 


( radians ) 


Relative  Complex  Permittivit\ 


L5r  '  11 


£xp  :=  (x-  I 


(dimensionless ) 


op  =  1  -  3  31288  10^1 


(dimensionless ) 


Horizontal  Reflection  Coefficient  (Fh)  : 


m  =0  —  (increments) 

1 


181 


Total  Electric  Field  (Ethor): 


Ethorie)  =  Jl(k  a  sin(0))  (l  -  rh(e)  e'^  2kh«»(e)) 

2rff  / 


Ethor  1 

m 


■J 


rhl  e 

m 


Radiation  Intensity  U(9): 


U(0)  :  =  —  (  |Ethor(e)| 
‘’lo 


(W/solid  angle) 


Umax  ■= - (  |max(Eihorl  )| 

2% 


(W/solid  angle) 


Umax  =7  03652  10‘ 


(W/solid  angle) 


Radiated  Power 


Prad  := 


1  Ji 


U(  9)  sin(9)|  d0d4 


(W) 


Prad  =  1  13445*  10"' 


(W) 


Directivity  : 


Do:= 


4  n  Umax 
Prad 


(dimensionless ) 


Do  =7.79439 


(dimensionless ) 
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( V  /  m  j 


( v/m; 


Radiation  Resistance  (Rr)  : 


r,  .  2  Prad  ^,2 
Rr  * - N 


(ig)- 


Rr  =  2  2689-10- 


Bandwidth : 


Bandwidth  :=  - 
a 


Bandwidth  =6  17143  10* 


Surface  Impedance  of  Conductor  R^: 


Rs 


Rs  =4  71058-10"^ 


Ohmic  Resistance  Rf^^rm  r  : 


Rohmic  :=.^LL^  (RJR+  | ) 
b 


Rohmic  =0  04984 


Conduccion-Dielect  ric 


Ef ficiencyiSf-H*  ' 


ccd 


Rr 

Rohmic  -►  Rr 


ecd  =0  99998 
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(Q) 


(Q) 


(Hz) 


(Hz) 


(Q) 


(Q) 


(Q) 


(Q) 


(dimensionless ) 


(dimensionless ) 


Maximum  Effective  Aperture  (A^^l : 


Aem  ■» 


£Cd  A.^  Do 


Aem  =  0  52452 


Gain  G : 


G  :=  ccd  Do 


(dimensionless ) 


Gdb  :=  10  log(ocd  Do) 


G  =  7  79422 


(dimensionless) 


Gdb  =  8  91773 


Maximum  Effective  Height  (hom) : 


hem  =  2 


hem  =  3  55347 


Effective  IsotroDic  Radiated  Power  (EIRP) : 


EIRP  ::  Prad  Do 


EIRP -8  84235  lO'’ 
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LAJlGg  HORIZONTAl.  LOOP  ANTENNA  FiVR  FIELD  PATTglW  ABOVE  GROUND  WITH  AIR 
CORE 


For  the  purpose  of  this  far-fieid  radiation  pattern,  the  large  loop 
antenna  axis  is  equivalent  to  the  Ex  =  0  grid  line.  The  pattern  is 
symetric  when  rotated  in  the  (^)  dirction  about  the  antenna's  axis. 

Large  Horizontal  Loop  with  Air  Core  Electric 
Over  Real  Ground 


E\(  0)  =  |Ethor(0)| 


cos|  0  +  - 

T 


Ey(0)  -  |Eihor(0)|  sin[0-f  — 

7i 


-0  005  -0  004  -0  003  -0  002  -OOOl  0  0  001  0  002  0  003  0  004  0  005 

E.\(0) 


E  FIELD 
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Calculate  large  horizontal  loop  ar.tenr.a  parameters  ovr  ground  (fTrxt« 
cor*)  : 


Total  Electric  Field  (Ethor): 


Ethor(  0) 


(^1 


e-o  ,11(1,  a  simG))  (l  -  rh{0)  e'J 


i  V  /  rr. ; 


Elhor  1 

m 


•J 

rhl  e 

m 


(V/m) 


Radiation  Intensity  'J(9): 


rlT^  , 

U(0)  :  =  —  ( |Elhor(9)|  )* 


(W  /  solid  angle! 


UiTiaN  : - ( |ma.\(  Ethor  1  )|  )' 

-  % 


(W  /  solid  angle) 


Umax  =  1  13486  lo''* 


(W  /  solid  angle; 


Radiated  Power 


Prad  := 


1  n 


U(0)  sini  0)|  (J0d0 


IW! 


Prad  =  I  82966  lo'® 


(W) 
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Directivity  D^, : 


_  4  n  Umax 

Do  •= - 

Prad 


Do  =7.79439 


Radiation  Resistance  (Rr )  : 


-  Prad  ^^2 

(I'or 


Rr  =  3.65931  10'° 


Surface  Impedance  of  Conductor  Rg: 


Rs  =  0.29852 


Ohmic  Resistance 


.  o .  D  c 

Rohmic:=‘  (RR-t-  1) 
b 


Rohmic  =3.15832 


(dimensionless ) 


(dimensionless) 


(Q) 


(Q) 


(Q) 


(Q) 


(Q) 


(Q) 
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(dimensionless ) 


Conduction-Dielectric  Ef f iciencyjCf.^)  : 


ecd 


Rr 


Rohmic  +  Rr 


ECd  =  1 


Gam  G* 


G  sEcd  Do 


(dimensionless ) 


(dimensionless) 


Gdb  :=  10  iog(ecd  Do) 


(dB) 


G  =  7  79439 


Gdb  =  8.91782 


(dimensionless ) 


(dB) 


Maximum  Effective  Aperture  (Aam)  •' 


Aem  ■= 


Ecd  1^  Do 
4-1! 


(m^) 


Aem  =0  52453 


(m^ ) 


Effective  Isotropic  Radiated  Power  (EIRP) : 


EIRP  Prad  Do 


(W) 


EIRP  =  1.42611  lO" 


(W) 


88 


Bandwidth : 


Bandwidth  :*  - 
a 


(Hz) 


Bandwidth  =6.17143- 10* 


(Hz) 


Maximum  Effective  Height  (h^m) : 


,  [Rr  Aem 
hem  =  I 

^  % 


2 


(m) 


hem  =  1.42709-1  o’ 


(m) 


THE  BEDSPRING  ANTENNA 
MATHCAD  SOFTWARE -BEDS PRIN.MCD 

A  bedspring  (or  curtain)  antenna  is  a  two  dimensional  array  of  identica 
horizontal  dipoles.  The  vertical  stacks  of  dipoles  are  referred  to  as 
bays.  A  bedspring  antenna  is  built  from  two  or  more  bays.  Bedspring 
assemblies  are  normally  designed  for  high  frequency  (3-30  MHz) 
operations.  The  beam  maximum  can  be  steered  in  either  the  azimuthal  or 
vertical  directions  by  adjusting  the  phase  of  the  element  feed  currents 
Bedspring  antenna  Mathcad  applications  assume  that  all  elements  lie.  in 
the  y-z  plane,  that  all  element  excitation  currents  have  identical 
maximum  amplitude,  that  there  are  only  2  dB  of  losses  associated  with 
the  antenna,  and  that  a  perfect  reflector  screen  is  located  in  the  -x 
half  space. 

(Note;  Mathcad  equations  cannot  use  symbolic  subscripts.  Therefore, 
symbols  like  a.  will  immediately  follow  the  parameter  in  equations  in 
lieu  of  subscripts.  ) 

The  bedsprinq  antenna  Mathcad  application  will  compute  the  following 
parameters : 

D  =  Maximum  Physical  Dimension  of  the  Array 
k  =  Wavenumber 
/.  =  Wavelength 

Dq  =  Directivity 

E(J2  =  Electric  Field  (9)  Component  for  an  Individual  Bay 

E^2  =  Electric  Field  (^)  Component  for  an  Individual  Bay 

Sjj  =  Array  Factor  for  Perfect  Image  Reflector 
Sy  =  Array  Factor  for  Multiple  Bays 
A,  B  =  Electric  Field  Coefficients 

Egt  =  Electric  Field  (9)  Component  (Including  Ground  Effects) 

E^t  =  Electric  Field  (♦)  Component  (Including  Ground  Effects) 

U  =  Radiation  Intensity 

Umax  ~  Maximum  Radiation  Intensity 

^rad  “  Radiated  Power 

Tv  =  Vertical  Reflection  Coefficient  (Ground  Reflection) 

=  Horizontal  Reflection  Coefficient  (Ground  Reflection) 

G  =  Gain 

EIRP  =  Effective  Isotropic  Radiated  Power 
^em  ~  Maximum  Effective  Aperture 
BW  =  Bandwidth 

Ej-1  =  Ground  Relative  Complex  Permittivity 

(Note:  The  subscript  '  will  be  annotated  as  p  in  the  application) 

^min  “  Minimum  Distance  to  Far-Field 

Rj-  =  Radiation  Resistance 

hem  ~  Maximum  Effective  Height 

IIqI  =  Magnitude  of  Antenna  Feed  Current 

^x,y,z  ~  Electric  Field  Spatial  Components 

Oa  =  Antenna  Unit  Polarization  Vector 

^high  ~  Upper  Operating  Frequency 

^low  ~  Lower  Operating  Frequency 
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PLF  -  Polarization  Loss  Factor 

<ja  =  Antenna  Unit  Polarization  Vector 

8p  =  Coaltitude  (Deflection  Angle  from  +z  Axis)  for  Polarization  Loss 
=  Azimuth  Angle  for  Polarization  Loss 


The  following  data  must  be  input  based  on  known  or  estimated  data: 

M  =  Number  of  Elements  per  Bay 
N  =  Number  of  Bays 

=  Height  of  ith  Element  Above  Ground 

=  Vertical  Spacing  of  the  ith  Element 
=  Horizontal  Position  (Center  of  Dipole)  of  i^^  Bay 

=  Horizontal  Spacing  of  the  i*-^  Bay 
Xi  =  Reflector  Position 
1  =  Half-Length  of  Each  Element 
f  =  Frequency  of  Interest 

1  =  Number  of  Increments  for  Far  Field  Radiation  Patterns 
(Note;  Due  to  Mathcad  restrictions  on  matrix  size,  i  <  125) 

0  =  Ground  Conductivity 

=  Ground  Relative  Permittivity 
Cff  =  Distance  of  Field  Calculations 
Ifjj  =  Input  Current  at  Element  Terminals 

9g  =  Coaltitude  (Deflection  Angle  from  +2  Axis)  for  Azimuth  Plot 
♦g  =  Azimuth  Angle  for  Elevation  Plot 
C^  =  Relative  Current  Amplitude  of  i^^  Element 
Gq  =  Vertical  Scan  Angle  (Coaltitude) 

♦q  =  Azimuthal  Scan  Angle 

(x,y,z)  =  Coordinates  for  Unit  Polarization  Vector 
dy,  =  Incoming  Wave  Unit  Polarization  Vector 
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Enter  input  data  here: 


M  :=4 


(elements ) 


f;=  10- 10° 


N  :=2 


(bays) 


ET  ;=  10 


1  =90 


( increments ) 


rff  =  1  10" 


Z  := 


Y  = 


C  :  = 


53  55 
42  42 
29  29 
16  16 

'O  26 
0  26 
0  26 
0  26 

1  1 
1  1 
1  1 
1  1 


av,  = 


_1_ 

/2 

_1_ 

/2 


(  —  ) 


(m) 


{ amps ) 


(dimensionless ) 


Im  :=  1 


0  :=  01 


XI  :=7 


1=11 


0o:=- 

2 


4o:=0 


X  :=  1000 


V  =  50 


z  •=  50 


(HZ) 


(dimensionless ) 


(meters ) 


( amp  s ) 


(mhos  /  m) 


(m) 


(m) 


( radians ) 


( radians ) 


(m) 


(m) 


(m) 
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loulate  bedsprinq  antenna  geometric  parameters  and  define  constants; 


c  =2.9979  10*  (meters/sec) 


=  120  n 


(meters/cycle) 


£  =_L  Iff'^ 

°  36  s 


( Farads/m) 


=29979  (meters/cycle)  ^^:=4s1(T"  (H/m) 


D  =60.83584 


Calculate  bedspring  antenna  parameters  : 


Define  angular  offset  6 
from  y-z  axis: 


Distance  to  Far-Field  r- 


0:sO,— ..- 

2i  2 


( radians ) 


ro:=1.6X 


r,  ;=5  D 


S  /  JT  ^  £ 

'"'li' 2^1/  2 


Wavenumber  k: 


( radians ) 


r  -i-2- 


mun  :=  max{  r) 


rmin  =3.04179  10  (m) 


(m“^) 


k  =  0.20959 


(m"l) 
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Ground  Relative  Complex 


Permittivity 


erp  ■=  ET- 


(dimensionless ) 


op  =  10-  l8j 


(dimensionless ) 


.^ertical  Reflection  Coefficient  f',: 


w  :  =  0 

T 


(increments) 


rvx0)  -  °P  ygp-  sin(9) 

ap  cos( 9) -h  Jap-  sm(0)^ 


(dimensionless ) 


rvl  .= 


op  cosf~j  -  Jap- 


(dimensionless ) 


ap  cosi — ^1  +  Jap-  sinl— — I 


Horizontal  Reflection  Coefficient  Fh: 


ni(0) ..  cos(9)-  /jgp-  sin(9)~ 
cos(0)  + Jap-  sin(0)^ 


(dimensionless ) 


nil  ■=- 


!ili 

(t1-^ 


I  /  w  :i 

(op  -  sm  — 


(dimensionless ) 


cosra.Jap-sinra 
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Electric  Field  Coefficients  A,B: 


( increments i 

q  =  I  M 


A(0) 


q 


kZ 


q-  1,0 


( co«(  6 )  -  cotl  do ) ) 


1.0 


-J 


rv(e)e 


1.0 


(dimensionless 


A1 


Vk 


1,0 


e 


-  co«(  00)1 


) 


q 


-j 


-  Fvl 


2“q-1.0“* 


H) 


(dimensionless 


8(9)  :=2c,. 


kZ 


q-  1.0 


(  C<M(  6 )  -  COl(  do )  ) 


1.0 


/  -j  2  k  Z  ,  »  co»<0)\ 

ln.rh(0)e  } 


(dimensionless ) 


B1 


1,0 


-  coi(  00  H  [ 


•  J  2  k  Z 


Thl 


0- 


q 


(dimensionless ) 


Electric  Field  Components  for  Individual  Bay  E0nil 


E0i(y.^)  ;= 


rff 


1  -  sin(0) 
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(V/m) 


e-J  k-ff 

60  - 

rff 


cos(  k  1  sin(  9 )  sin(  4 ) )  -  cos  k  1 ) 
1  -  sin(9)*  sin( 


cos(^)  B(0) 


i  V/m) 


Array  Factor  for  N  bays  S„: 


ql  :=  1  N 


( increments  i 


sirt  9)  (  sirt  «)-  sin(  *0)) 


ql 


(dimensionless ) 


Svi 


w.v 


ql 


‘‘'O.ql-l’*"' 


(dimensionless ! 


Array  Factor  for  Ideal  Reflector  Sy: 


Sx(9»:=l-eJ  2  “ -Xi  s.n(9)  «»(#) 


(dimensionless 


-j  2  k  Xl  sin 


Sx  1  '=1-6 

W.V 


(t1 


(dimensionless  / 
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Total  Electric  Field  Components  En^  ^  : 


E0t(0.*)  :=  £61(0, ♦)  Sx(0.4)  Sy(0.4) 


(V/m) 


E0tl  :»E011  Sxl  Svl 

w  .  V  ^  W  ,  V  -  Vli  ,  V 


(V/ir.) 


E»t(e,#)  a  E*U0,*)  Sx(0,#)  Sy(0.*) 


(V/m) 


E*tl  ::E*11  Sxl  Svl 

Vi  .\  Vi 


(V/m) 


Radiation  Intensity  U(0): 


U(0,*) 


:zI^.  l(|E0t(0,*)|r+(lE*t(0.*)|  )^1 


(W/solid  angle) 


U1 


(W/solid  angle) 


Umax  -  max(  U 1 ) 


(W/solid  angle) 


Umax  =  8  8082’  lO’ 


(W/solid  angle) 


Radiated  Power  PraH^ 


rt 

-> 


Prad  := 


1 

1 


U(0,4)  sin(0)  dOd^ 

0 


■> 


Prad  =9.79181  10^ 


(W) 


(W) 
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^  4  n  Umax 

Do  = - 

Prad 


(dimensionless ) 


Do  =  I  13041  10^ 


(dimensionless ) 


Magnitude 


of  Antenna  Feed  Current 


lo  ::M  N  Im 


( amps ) 


lo  =8 


( amps ) 


Radiation  Resistance  (Rr)  : 


Rr:= 


2  Prad 
(|lo|)^ 


(fi) 


Rr  =  30.5994 


G3in  Gr 


(Q) 


Gdb  :=  10  log(Do)-  2 


Gdb=  18  53234 


(dB) 


(dB) 


G  =  10 


Gdb 

10 


(dimensionless ) 


G  =71  32377 
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( dimensionless ) 


Polarization  Loss  Factor  PLF : 


0p  •=  atan 


( radians ) 


:= 


( radians ) 


ep  =  1.5209 


(radians)  4P“  0^4996 


( radians ) 


Ex  :=Eet(0p,4p)  cos(6p)  cos(4p)-  E4t(0p>p)  sm{^) 


(V/m) 


Ey  .-E9t(6p>p)  cos(0p)  sin(^p)-  E^t(0p,^p)  cos(^) 


(V/m) 


Ez  •=•  E9t(0p,^)  sin(0p) 


(V/m) 


/()Ex|)^+(|Ey|  )V(fEz|)^ 


(dimensionless ) 


-0  03302  +  0  03754J 
rta  =  ”0  66366  +  0.74634j 

-3  331  12-10~-’  -  4  40961  10"’j 


( dimensionless ) 


;=(lcjwoa|) 


(dimensionless ) 


PLF  =0.54991 


(dimensionless ) 


Maximum  Effective  Aperture  (A, 


Aem  ■  = 


Do  PLF 


Aem  =4  44582  lO'’ 
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Effective  Isotropic  Radiated  Power  (EIRP)  : 


EIRP  =PradC)o 


EIRP  =  1.10687  lo’ 


Maximum  Effective  Height  : 

■..pLiSSl.2 


hem  =  37.99241 


(m) 


Bandwidth  BW: 

Half-Wave  Assembly 


Symmetric  Feed  Assembly 


(Hz)  fhighl  1.5  112 


f/.2  =6  81341  10* 


(Hz) 


flow  I  =  98  112 


flugh  :=  1.02  fit2 


flow  ;=  98  fl2 


BW  =fhigh-  flow 


BW  =  2.72536  lo’ 


BWl  =  fhighl  -  flowl 


(Hz)  BWl  =3.54297  10* 


(Hz) 


(Hz) 


(Hz) 

(Hz) 

(Hz) 

(Hz) 
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BSDSPRIMG  AMTENNX  AZIMUTHAL  FAR  FIELD  PATTERN 


For  the  purpose  of  this  far-field  radiation  pattern,  the  bedspring 
antenna  lies  on  the  Ey  =  0  grid  line.  Three  field  patterns  are 
developed:  th*  (6)  coiB|>on*nt  of  oloctric  field,  the  (^)  conponent  of 
electric  field,  the  total  electric  field.  The  user  must  select  the 
desired  coaltitude  (0g)  for  which  the  patterns  will  be  graphed. 


Theta  Component  of  Azimuthal  Electric  Field 


E0tx(4)  :=  |E0t(0g,4)|  cos|4+-^|  Eety(4)  :=  |Eei(eg.4)|  sin|^  +  -^ 


Ee«y(»)  2  10 


—  E  FIELD 


Phi  Component  of  Azimuthal  Electric  Field 


E*tx(*)  ;=  lE*t(eg,*)|  cos|*  +  -iJ  E^ty(4)  :=  lE*t(eg,*)|  sin|f 


1.35* 


E«ty(*)  7.5  10 


4,5‘10”'* 


—  E  FIELD 


Total  Azimuthal  Electric  Field 


Et(*)  ;=4(  |E0t(0g.*)|  )^t  ( |E*t(eg,*)| 


Etx(*) 


•=  |Et(4)|  cos|*4--^ 


Ety(*) 


■=  |Et(*)|  sin|* 


Ety(#)  8  10 


10  ‘“-ISMO  ‘“•1.2’10  ■‘‘-e’lO  ‘^1.926M0  ^6  10  1,2M0  ‘“l.SMO  ‘®2.4M0  3n0 

E(x(*) 

—  E  FIELD 


BED3PRIMG  AMTgNNA  ELEVATION  FAR  FIEIJ)  PATTERM 


For  the  purpose  of  this  far-field  radiation  pattern,  the  bedspring 
antenna  lies  on  the  Ex  =  0  grid  line.  Three  field  patterns  are 
developed:  the  (0)  component  of  electric  field,  the  ii)  component  of 
electric  field,  the  total  electric  field.  The  user  must  select  the 
desired  aximuth  angle  (^g)  for  which  the  elevation  patterns  will  be 
graphed. 


Theta  Component  of  Elevation  Pattern  of 
Electric  Field 


♦g  -0 


E9tx(e)  := 


cos(6) 


E0ty(e)  := 


sm(9) 


1 

8 

6 

4 

2 

17 

.2 

0 

1 

1 

t 

1 

1 

1 

1 

1 

t 

1 

1 

.6 

1 

1 

• 

-1 

1 

—  E  FIELD 


Efiix(e) 


Phi  Component  of  Elevation  Pattern  of 
Electric  Field 


E4tx(6)  ;= 


cos(0) 


E*ty(0)  := 


sin(0) 


-0.005  -0.002  0.002  0.006  0.009  0.01?  0.016  0.02  0.023  0.027  0.03 


—  E  FIELD 


EMxl  9) 
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Elevation  Plot  of  Total  Electric  Field 


Et(0)  :»J(lE0t(e.te)l)^t(|E4t(8.4g)| 


Etx(0)  •= 


Et|— -0|  cos(0) 


Ety(8)  •= 


sm(0) 


-0.003  -0.002  0.002  0.006  0.009 


0.023  0.027  0.03 


E  FIELD 
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0.013  0016 

Etx<0) 


0.02 


THE  SPIRAL  ANTENNA  (EXACT  METHOD) 

MATHCAD  SOFTWARE -S PI RAL.MCD 

Spiral  antennas  are  a  feunily  of  two  or  three  dimensional  devices  that 
possess  frequency  independant  par^uneters  over  a  wide  bandwidth.  Spiral 
antennas  are  commonly  used  for  direction  finding,  satellite  trac)ting 
and  missile  guidance. 

The  planar  spiral  may  be  of  the  Archemedean,  log-spiral,  or  equiangular 
type.  All  three  radiate  two  main,  circularly  polarized  lobes 
perpendicular  to  the  plane  of  the  antenna.  Additional  gain  for  planar 
spirals  may  be  achieved  by  placing  a  metal  cavity  on  the  side  of  the 
antenna  with  the  unwanted  lobe.  The  cavity  may  be  ertpty  or  be  filled 
with  electromagnetic  energy  absorbing  material.  These  applications 
principly  examine  the  equiangular  planar  spiral  and  do  not  account  for 
cavity  backed  effects. 

The  three  dimensional,  or  conical,  spiral  exhibits  many  of  the  same 
features  as  the  planar  spiral  except  that  it  radiates  a  single  main  beam 
in  the  direction  of  its  tip,  thereby  eliminating  the  need  for  cavities. 
The  conical  log-spiral  is  the  only  three  dimensional  antenna  analyzed  by 
this  application. 

(Note:  Mathcad  equations  cannot  use  symbolic  subscripts.  Therefore, 

symbols  like  1  will  immediately  follow  the  parameter  in  equations  in 
lieu  of  subscripts.  ) 

The  spiral  antenna  Mathcad  applications  will  compute  the  following 
parameters  for  equiangular  planar  spirals  and  conical  log-spirals: 

k  =  Wavenumber 
X  =  Wavelength 
Dq  =  Directivity 

E0  =  Electric  Field  (0)  Component 
E^  =  Electric  Field  (^)  Coit^onent 

A, B,c  =  Conical  Log-Spiral  Electric  Field  Coefficients 

U  =  Radiation  Intensity 

^max  ~  Maximum  Radiation  Intensity 

P^ad  “  Radiated  Power 

G  =  Gain 

EIRP  =  Effective  Isotropic  Radiated  Power 
^em  “  Maximum  Effective  Aperture 
BW  =  Bandwidth 

^min  ~  Minimiam  Distance  to  Far-Field 

R^  =  Radiation  Resistance 

hem  “  Maximum  Effective  Height 

^x, y, z  ~  Conical  Log-Spiral  Electric  Field  Spatial  Components 
Og  =  Conical  Log-Spiral  Unit  Polarization  Vector 
^high  ~  Upper  Operating  Frequency 
^low  ~  Lower  Operating  Frequency 

tfi  =  Any  point  on  the  n^^  edge  of  a  spiral 
£ex  "  Equiangular  Planar  Spiral  Expansion  Ratio 
Zi  =  Planar  Spiral  Input  Impedance 
r  =  Voltage  Reflection  Coefficient 
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=  Reflection  Efficiency 

PLF  =  Polarization  Loss  Factor 
Q  =  Conical  Log-Spiral  Antenna  Slowness  Factor 
I  (5)  =  Conical  Log-Spiral  Current  Distribution 
4(^)  =  Conical  Log-Spiral  Azimuth 

=  Conical  Log-Spiral  Phase  Difference  of  n^^  arm 
/thigh  ~  Upper  Operating  Wavelength 
Iiqw  =  Lower  Operating  Wavelength 
P  =  Planar  Spiral  Angle 

A  =  Planar  Spiral  Electric  Field  Amplitude 
b  =  Conical  Log-Spiral  Constant 
L  =  Conical  Log-Spiral  Total  Arm  Length 

=  Conical  Log-Spiral  Azimuth  at  End  of  Arm 
9p  =  Desired  Conical  Log-Spiral  Polarization  Offset  Angle  from  z  Axis 
^  =  Desired  Conical  Log-Spiral  Polarization  Azimuth  Tvngle 
L  =  Conical  Log-Spiral  Total  Arm  Length 

The  following  data  must  be  input  based  on  known  or  estimated  data: 

M  =  Mode 

N  =  Number  of  Spiral  Arms 
f  =  Frequency  of  Interest 

i  =  Number  of  Increments  fo*-  Far  Field  Radiation  Patterns 
rff  =  Distance  of  Far-Field  Calculations 

Iq  =  Input  Current  at  Antenna  Terminals 
9g  =  Deflection  Angle  from  +z  Axis  for  Azimuth  Plot 
=  Azimuth  Angle  for  Elevation  Plot 

(x,y,z)  =  Coordinates  for  Conical  Log-Spiral  Unit  Polarization  Vector 
rQ  =  Spiral  Feed  Point 
a  =  Flare  Rate 

5n+i  =  Angular  Arm  Width  of  n^^  Spiral  Arm 
^ex  ~  Azimuth  to  Confute  Expansion  Ratio 

P  =  Conical  Log-Spiral  Angle 
R  =  Overall  Radius 

Eq  =  Source  Strength  Constant  for  Planar  spirals 
=  Conical  Log-Spiral  Input  Iit^edance 
Oyf  =  Wave  Unit  Polarization  Vector 
9q  =  Cone  Angle 

A9  =  Conical  Log-Spiral  Half-Power  Beamwidth 
Zq  =  Characteristic  Impedance  of  Feed  Assembly 

t  =  Number  of  Increments  Along  Conical  Log-Spiral  Arm 
Og  =  Equiangular  Planar  Spiral  Unit  Polarization  Vector 
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TBE  PLANAR  SPIRAL  ANTENNAS 


Enter  input  data  here: 

N  = 

2 

(arms) 

f  =3  10* 

(Hz) 

M  2 

:  1 

(mode) 

lo  =  1 

{amps) 

(Note:  is  N-1) 

rfif  =1 

10^ 

(meters) 

i  =30 

a  =.221 

(increments) 

(dimensionless ) 

ro  =. 

1 

(m) 

Eo  = 

10^ 

(V/m) 

R  = 

1 

(m) 

1 

/2 

j 

^■2 

0 

8  = 

■  0 

n 

2 

n 

! 

(radians) 

ow  = 

(dimensionless ) 

3-b 

2 

Zo  = 

100 

(Q) 

Calculate  planar  spiral  antenna 

geometric  parameters 

and  define 

constants: 

c  =2.9979- 10* 

(meters/sec) 

11^  =120-n 

(Q) 

1  =^ 
f 

(meters/cycle) 

* 

36n 

(Farads/m) 

1  =  99.93 

(meters/cycle) 

Hjj  =4iil0'’ 

(H/m) 

oa  = 


j_ 

/2 

0 


(dimensionless) 
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Calculate  planar  spiral  antenna  parameters  : 


Define  angular  offset  9 
from  y-z  axis: 


Distance  to  Far-Field 


n  It  It  n 

2’i  11 


ITg  =  1.6  1 

(radians) 

rr|  =  10  R 


(m) 

(m) 


=  0,  — , 


2n 


(radians) 


(m) 


rmin  =max(n-)  (m) 
nnin  =  1.59888*10^  (m) 


Wavenumber  )c: 


^  (ni“l) 

\ 


k  =  0.06288 


(m”l) 


Radial  Distance  to  Spiral  Edge  r: 


a* 

r(n,^r)  -roe 


(m) 


r(1.2  1)  =0.40092 

Expansion  Ratio  gay: 


£ex(n,4r) 


T<n,^r-^  2  ii) 
Kn.^r) 


eex(l,2  n)  =4.00917 


(ro) 


(dimensionless ) 


(dimensionless) 


£mini 
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Bandwidth  BW: 


Equiangular  Spiral 

Log-Periodic  Spiral 

Hugh  =4ro 

(m) 

Hughl  =20ro 

(m) 

fhigh  =  — - — 

fhighl  = — ^ — 

Hugh 

(Hz ) 

Hughl 

(Hz ) 

fhigh  =7.49475*10* 

(Hz) 

fhigh  =7  49475*10* 

(Hz) 

How  =4  R 

(m) 

flow  = - 

(Hz) 

How 

flow  =7.49475*10^ 

(Hz) 

BW  =  fliigh  -  flow 

(Hz) 

BWl  =fhighl-flow 

(Hz) 

BW  =6.74528*10* 

(HZ) 

BWl  =7.49475*10’ 

(HZ) 

Electric  Field  E  (9i)  and  Electric  Field  Amplitude  A(6)  : 


w  =0..i  (increments) 


EK*.0)  =■ 


Eo-k^cos(0)(l  +  j  8  005(0)) 


(V/m) 


sin(0)  rfif 


E* =5.97188*10  *  +2.4182*10 

\6  6j 


(V/m) 
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Electric  Field  Amplitude  A: 


A(0)  - 


/0\M  /— V»«*n(»'COi(9)) 
cos(0)  tan  - 

\2l 


sin(0)  Vl  a^  cos(0)^ 


(V/m) 


1  n  w\ 

-  —  }i-  — ' 

M 

/\ 

I  a  w) 

-tani 

2  i 

/  \  2  i/y 

cosl-  -  -1-  JI-  — ! 
\  2  1  / 

2  1 

•e 

I-  1 

1,  2  'It  w\ 

1  -f  a  COS'-  -  T  n — 

\  2  ii 

\  2  i/ 

(V/m) 


Radiation  Intensity  U(0): 


U(0)  =-L  (a(0)^; 

2% 


Umax  =max(Ul) 


Umax  =2.2285*10 


Radiated  Power  Prart- 


I  2 

Prad  =4  n  j  U(0)  sin(0)d0 

Jo 


(W  /  solid  angle) 


(W  /  solid  angle) 


(W  /  solid  angle) 


(W  /  solid  angle) 


(W) 


Prad  =7.94528*10 


(W) 
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(dimensionless) 


Do  = 


4  I  Umax 
Prad 


Do  =  3  52462 

Radiation  Resistance  Rr; 


(dimensionless ) 


Rr  = 


2  Prad 
( :io;)^ 


(Q) 


Rr  =0  01589 

Input  Impedance  Zi : 


Zi  = 


N  30  « 


M'. 

smin-  —  1 
N/ 


Zi  =  1  884%'10^ 

Voltage  Reflection  Coefficient  f: 


(Q) 


(«) 


(£2) 


r  = -  (dimensionless) 

Zi  ^  Zo 


r  =0.30675  (dimensionless) 


Reflecton  Efficiency 


erv  r  1  _  (  ji  (dimensionless) 


arv  =0.90591 
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(dimensionless) 


G  ^trv  Do 


(dimensionless) 


G  =  3.19298 


(dimensionless) 


GdB  =101og(G) 


GdB  =5.041% 


Effective  Isotropic  Radiated  Power  (EIRP) : 


EIRP  =PradDo 


EIRP  =0.028 


Polarization  Loss  Factor  (PLF) 


=  (Itrw-cra  ) 


(dimensionless) 


PLF  =  1 


(dimensionless ) 


Maximum  Effective  Aperture  (Ap^) • 


Aem  = 


X  Do-ffv  PLF 


Aem  =2.53733-10^ 
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hem  =  0.65407 


(m) 
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THE  EQPIANGOIAR  PLAMAR  SPIRAL  ANTPOa  yAR-FIBLD  ELXVATIOH  PATTEEN 


For  the  purpose  of  this  far-field  radiation  pattern,  the  spiral  antenna 
lies  parallel  to  the  Ey  =  0  grid  line  and  is  centered  at  the  origin. 

The  magitude  of  the  electric  field  pattern  is  rotationally  synmetric 
with  respect  to  the  Ex=0  grid  line.  The  equiangular  planar  spiral 
antenna  possesses  a  mirror  image  radiaton  pattern  in  the  -y  half  plane. 


E4x(9)  =  lA(0)|cos(0r^ 


E*y(0)  -lA(0)l  sm  0+^ 
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Enter  input  data  here: 


N  =2 

( arms) 

M  =1 

(mode) 

(Note: 

Mitiax 

i  =9  (increments) 


R  =  15 


8  = 


(m) 


(radians) 


f  =  280  10* 
lo  =1 

rff  =1  10* 


ro  =  .03 


Zi  =160 


ow  = 


i2 

j_ 

Ti 

0 


P  =73 

(degrees ) 

P  =P  — 

180 

(radians) 

P  =  1.27409 

(radians) 

X  =0 

(m) 

y  =0 

(m) 

z  =10* 

(m) 

Zo  =100 


00  =10 


00  =  00— ^ 
180 

00  =  0.17453 


HPBW  =80 


(Hz) 

(eUT^S) 

(meters) 

(m) 

(Q) 

(dimensionless ) 

(Q) 

(degrees) 

(radians) 

(radians) 

(degrees) 
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Calculate  conical  log-spiral  antenna  geometric  parameters  and  define 
constants : 


c  =2.9979- 10*  (meters/sec) 


Tig  =  120  n 


(Q) 


Q 

1  =-  (meters/cycle) 


X  =  1.07068  (meters/cycle) 


-L.io-’ 
36  n 


( Farads/m) 


|1q  =4s10''  (H/m) 


Calculate  conical  log-spiral  antenna  parameters  ; 


Define  angular  offset  9 
from  y-z  axis: 


n  n  n  n 

2’i  22  (radians) 


rO,  — ..2-Jl 
i 


( radians ) 


Distance  to  Far-Field  rmin; 


rr^  =1.61 

(m) 

n-j  =10-R 

(m) 

8-R^ 

(m) 

nnin  =max(rr) 

(m) 

nnin  =  1.71309  (m) 


Wavenumber  k: 


k  = 


2ji 


k  =  5.86841 


(m  1) 

(m“l) 


Radial  Distance  to  n^  Spiral  Edge  r: 


b  =cot(P)  (dimensionless) 


b  =0.30573  (dimensionless) 


t(n,W  (■”! 


K  1.4-11)  =0.04948  (m) 
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Bandwidth  BW: 


Conical  Log-Spiral 


jLhigh  =4ro 

(m) 

fhigh  =  — ^ — 

(Hz) 

Hugh 

fhigh  =2.49825-10® 

(Hz) 

How  =  -  R 

(m) 

3 

flow  ^ 

(Hz) 

How 

flow  =7.49475-10* 

(Hz) 

BW  =fhigh-flow  (Hz) 


Antenna  Slowness  Factor  Q: 


Q  = 


Q=  3.4203 


(dimensionless ) 


(dimensionless ) 


Azimuth  at  End  of  the  First  Spiral  Arm  ^L: 


♦L 


1 

bsm(0o) 


*1  =  63.29231 


In 


sin(0o) 
\  ro 


(radians) 


(radians) 
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ro  ^^♦Lbtin(9o)_ 


L  ^-  le' 

b 


L=  2.72729 


Spiral  Am  Current  I  (^) 


t  =25 


(increments) 


^  =0.-..L 

t 


I(^)  =Ioe 


(amps) 


Electric  Field  Coefficients  g, A(^) ,  B  (IL)  >  C  (ij) ,  D  (£)  > t (£)  : 


(dimensionless) 


0  =  3.14159 


(dimensionless ) 


1  =0..N-  1 


(increments) 


V  =0..- 
2 


(increments) 


w  =0..i 


(increments) 


^l(c)  =_i - 


bsm(0o)  \ro 


(radians) 


C0('>,1)  = 


1  ,  j  ]  gi  •(♦U?)-*+lo)^  /j  j  \ 
bsin(0o)/  \  bsin(9o)/ 


(dimensionless) 


C18(l^.w.l)  =  1  + 


J 


b  sin(9o)/ 


1  - 


J 


-j  •  ♦l(4)---2'*-fl'a.) 


b  sin(0o)/ 


(dimensionless ) 


J 


b  sin(eo) 


•el  ' 


-  1 


j _ Vc'j  <♦•<4) -♦+'■«) 


b  sm(to)/ 


sin(0o) 


(dimensionless ) 


Cli('.w.l)  = 


J 


bsin(9o) 


•a) 


-  1- 


J 


bsin(0o) 


sm(0o) 


(dimensionless) 


ACM)  =5]®'^ 
I 


-  sin(e)  cos(8o)! 


(dimensionless ) 


DK^.v.w)  = 


|sin(9o)  cos|- 7t+ --|•C10(^,w,I) 

^ - — - sin[--n+  --[•cosCto)) 

\  2  2/ 


(dimensionless) 


Al(^,v,w)  =^e'-' ^  \  ^  \  ■  ^^  DK^.v.w) 

I 


(dimensionless) 


I 


0)iin(  6b)coi(  ♦)( 


(dimensionless) 


•CK5.il) 
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1 


2l  [  i  I.CUKU. 


1) 


(dimensionless ) 


Electric  Field  Compontents  E9,  E4: 


E0(0,*)  -j  fUg 


g-j  krfr 

2  rff  Q 


co«(  0)  •«»(  fto ) 


I(5)e 


A(5.0»d5 


(V/m) 


E01 


'--J  H, 


e-i 

2  rfF  Q 


UeS-c 


KO-c 


-I 


•C<M(dD) 


Al(5.v,w)d5 


JO 


(V/m) 


-j  k'lfr  l-^'Cai(9)'Cai(0b) 

2  rffQ  Jo 


(V/m) 


rL 


E*1  „ 


-j  kiff 


2  rffQ  Jo 


jJEi., 


I(5)e 


«»—»+■ 


aMeo) 


(V/m) 
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Radiation  Intensity  U (9)  : 


U2(0.*) 


=  i^[(lE0(e»|)^  +  (|E4(0.*)|)^l 
2lo 


(W  /  solid  angle) 


U22 


rtf 


v,w  2 


% 


Umax  =max(U22) 


Umax  =0.25969 


Radiated  Power  PraH^ 


Prad  = 


■2n 


Jo  Jo 


U2(0,*)sm(0)d0<H 


Prad  =0.81332 


Directivity  Dq! 


Do  = 


4  a  Umax 
Prad 


Do  =4.01237 


Ddbo  =  101og(Do) 


Ddbo  =  6.03401 


DIdbo  = 


32600 

HPBW^ 


DIdbo  =  5.09375 


(W  /  solid  angle) 


(W  /  solid  angle) 


(W  /  solid  angle) 


(W) 


(W) 


(dimensionless) 

(dimensionless) 


(dB) 

(dB) 

(dB) 


(dB) 
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Voltage  Reflection  Coefficient  f: 


Zi-  Zo 
Zi  Zo 


r  =0.23077 


Reflecton  Efficiency 


trv  =  l-(irj)^ 


erv  =0.94675 


Gain  G: 


G  =crvE)o 


G  =  3.7987 


GdB  =101og(G) 


GdB  =  5.79635 


Effective  Isotropic  Radiated  Power  (EIRP) : 


EIRP  =PradDo 


EIRP  =3.26336 


(dimensionless ) 

(dimensionless ) 

(dimensionless ) 

(dimensionless) 

(dimensionless ) 

(dimensionless ) 

(dB) 

(dB) 

(W) 

(W) 
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Polarization  Loss  Factor  PLF: 


6p  -atan 


jxKy^l 


(radians) 


^  =atan  - 

\x, 


( radians) 


0p=  0 


( radians) 


♦P=  0 


(radians) 


Ex  =E0(0p,^)  cos(9p)  cos(^)-  E4(0p.^p)  sm(4p) 


(V/m) 


Ey  -Ee(ep,^)  cos(ep)  sin(4p)-  E^Cep.^p)  cosC#) 


(V/m) 


Ez  =-E0(9p,4p)sin(0p) 


(V/m) 


/(lExl  )^^(|Ey|)V(|Ez|)^  \Ez/ 


(dimensionless ) 


/-0.92973  -  0.1613J 
oa  =  l  0,29553  +  0.14919j 


(dimensionless) 


PLF  =  (  crw  oal) 


(dimensionless ) 


PLF  =0.408% 


(dimensionless) 


Maximum  Effective  Aperture  (A^m) : 


Aem  = 


rOoffvPLF 


Aem  =0.14172 


Maximum  Effective  Height  (hgjjj) 


,  Rr  Aem  , 

hem  = - 2 


hem  =  4.88816*10 


PHI  COMPONENT  OF  THE  ELEVATION  PATTERN 


0 


0 


0 


-*•10  ^-2.711*10 


TOTAL  ELEVATION  PATTERN 

4g  =-  (radians) 

2 


Et(e.*)  =^/E6(0,*)^+E*(6.*)^ 


Ex2(0)  =  |Et(0+ n.4g)l  •008(0+ - 

2, 


Ey2(0)  =  lEt(O  +  ii,te)|sinf0+ ^ 


—  EFELD 
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TOTAL  AZIMUTH  PATTERN 

Ex2(*)  =  |Et(eg,*)|  cos(*)  Ey2(4)  =  lEt(0g.*)|  sin(*) 


Ex2(») 


EFIELX) 
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THE  CONICAL  HORN  ANTENNA  (EXACT  METHOD) 

MATHCAD  SOFTWARE -HORN_CON.MCD 

Conical  horn  antennas  are  devices  used  to  provide  a  transistion  from  a 
circular  waveguide  to  an  unbounded  medium  such  that  the  wavefront  at  the 
aperture  of  the  horn  has  nearly  a  constant  phase  at  any  point  in  the 
mouth  of  the  horn. 

Conical  horns  are  commonly  used  as  feed  elements  for  reflectors  used  in 
satellite  traclcing,  microwave  communications,  and  radar. 

(Note:  Mathcad  equations  cannot  use  symbolic  subscripts.  Therefore, 

symbols  like  1  will  immediately  follow  the  parameter  in  equations  in 
lieu  of  subscripts.) 

The  conical  horn  antenna  Mathcad  applications  will  compute  the  following 
parameters : 

k  =  Wavenumber 
1  =  Wavelength 

Dq  =  Directivity 

E0  =  Electric  Field  (0)  Con^onent 
E^  =  Electric  Field  (4)  Con^onent 
U  =  Radiation  Intensity 
Uj^ax  ~  Maximum  Radiation  Intensity 
=  Radiated  Power 
G  =  Gain 

EIRP  =  Effective  Isotropic  Radiated  Power 
Aem  =  Maximum  Effective  ;^erture 
BW  =  Bandwidth 

^min  “  Minimum  Distance  to  Far-Field 

Rj.  =  Radiation  Resistance 

hem  “  Maximum  Effective  Height 

^x,y,2  ~  Electric  Field  Con^onents  in  Cartesian  Coordinates 
(jg  =  Unit  Polarization  Vector 

^cte  ~  Transverse  Electric  Cutoff  Frequencies 
^ctm  “  Transverse  Magnetic  Cutoff  Frequencies 
Egp  =  Aperture  Efficiency 
PLF  =  Polarization  Loss  Factor 

Py(cos(0))  =  Associated  Legendre  Function  of  the  First  Kind 
Spherical  Hankel  Function  and  its  Derivative 
V  =  Legendre  and  Hankel  Function  Order 
bQ,5  =  Legendre  and  Hankel  Function  Constants 

^optimum  ~  Optimum  Conical  Horn  Mouth  Diameter 
^x,y,z  “  Conical  Horn  Magnetic  Vector  Potentials 
h  =  Conical  Horn  Axial  Height 

Bfj  =  Magnetic  Vector  Potential  Integral  Coefficients 
P  =  Magnetic  Vector  Potential  Integral  Phase  Shift 
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The  following  data  must  be  input  based  on  known  or  estimated  data: 
m, n  =  Mode 

f  =  Frequency  of  Interest 

i  =  Number  of  Increments  for  Far  Field  Radiation  Patterns 
rff  =  Distance  of  Far-Field  Calculations 

Iq  =  Input  Current  at  Antenna  Terminals 

8g  =  Coaltitude  (Deflection  Angle  from  +2  Axis)  for  Azimuth  Plot 
4g  =  Azimuth  Angle  for  Elevation  Plot 
(x,y,z)  =  Coordinates  for  Unit  Polarization  Vector 
a  =  Flare  Angle 

a  =  Circular  Waveguide  Inner  Radius 

pj  =  n^^  Zero  of  Bessel  Function  of  the  First  Kind,  Order  m. 

X'm,  n  "  Zero  of  Bessel  Function  Derivative  of  the  First  Kind,  Order 

m. 

Bq  =  Electric  Field  Amplitude  Constant 

djneas  ~  Measured  Diameter  of  the  Conical  Horn's  Mouth 


n  =  1 


m  =  1 


Enter  input  data  here; 


(mode  number) 


(mode  number) 


f=l.%10’  (HZ) 

lo  1  (amps) 


i  -18  (increments) 

a  =  ,0445  (m) 

a  =2-  (radians) 

9 

(Note:  a  must  b«  less  than  n/3) 


Bo  =1  (V/m) 

dmeas  =  I  (m) 


rflF  =  1  ’0^  (meters) 


ow 


J_ 

r 

V2 

j_ 

0 


(dimensionless) 


X  =  10^  (m) 

y  =  10^  (m) 

z  =  10^  (m) 
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Calculate  conical  horn  antenna  geometric  parameters  and  define 
constants : 


c  2  9979  10* 

(meters/sec) 

Hq  =  120  « 

ID) 

1  ^  = 
f 

(meters/cycle) 

e  ^ -  10 

°  36  n 

(Farads/m) 

X=  0.1 5295 

(meters/cycle) 

»Iq  -4n  lo’ 

(H/m) 

1 

dmeas 

doptimum  '  ~ 

.  /a\ 

h  =  ^ 

(m) 

atanl- 

\2l 

doptimum  =  0.6708 1  ( m) 

h  =  1.48879 

(m) 

Define  Bessel  Function  Matrices: 

(Note;  m  =  column  number,  n  =  row  number,  matrix  index  starts  at 
n  =  m  =0.  The  n  =  0  row  has  no  physical  significance,  it  is 
only  a  placeholder.) 


[  0 

0 

0 

0 

0 

2.4049 

3.8318 

5.1357 

6.3802 

7.5884 

X  J 

5.5201 

7.1056 

8.4173 

9.7610 

11.0647 

8.6537 

10.1735 

11.6199 

13.0152 

14.3726 

1  11.7915 

13.3227 

14.7960 

16.2235 

17.6160 

j 

0 

0 

0 

0 

0 

3.8318 

1.8412 

3.0542 

4.2012 

5.3175 

XP  = 

7.0156 

5.3315 

6.7062 

8.0153 

9.2824 

10.1735 

8.5363 

9.9695 

11.3459 

12.6819 

13.3237 

11.7060 

13.1704 

14.5859 

15.9641 

(dimensionless) 


(dimensionless) 
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[  0  0  0  0  0  1 

1  I 

411135-10’  1.97552-10’  3.27702-10’  4.5077-lo’  5.70544-lo’ 

fcte=  7.52742-10’  5.72046-10’  7.19545-10’  8.60005-10’  9.95959-lo’ 

(Hz) 

1.09157-10“*  9.15906-10’  1.06968-10“*  1.21736-10“*  1.36071- 10“* 

[  1.42957-10*'*  1.256*10“*  1.4l312*lo“*  1.565-10“*  1.71288-10*'* 
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Transverse  Magnetic  (TM)  Modes 


fctm  = - * 

2 


1 

■  0 

0 

0 

0 

0 

2.58035*  10® 

4.11135*10® 

5.51037*10® 

6.84567*10® 

8.14201*10® 

fctm  = 

5.92282*10® 

7.62399*10® 

9.03138*10® 

1.04731*10*® 

1.18719*10*® 

9.28503*10® 

1.091 57*  10‘° 

1.24676*10*® 

1.39647*10*® 

1.54211*10*® 

1.26517*10*® 

1.42947*10*® 

1.58754*10*® 

1.7407!  *10*® 

1.89012*10*® 

Legendre  and  Hankel  Fuction  Constants 


log' 

i 

\  4  / 

log( 

\ 

1 

cos  - 

\2// 

(dimensionless ) 


\  =  26.32708 


(dimensionless) 


5 


jt 


a 


(dimensionless) 


6=4.5 


(dimensionless) 


Legendre  and  Han)cel  Function  Orders  y: 


V  -  -  .5  +  -S-  Jl  -t-  4  b 

0 


(dimensionless ) 


v  =  4.6553 


(Hz) 


1 


(Hz) 
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(dimensionless ) 


Legendre  Fuctions  P„cos  (6)  : 


pKG)  =  ^(''•^2) 


-  sin(0)  coa  (v+  .5)  0  + 
r(vr  1.5)  ■42 


(dimensionless ) 


=  -1.51294 


(dimensionless ) 


Hankel  Function  and  its  Derivative  ; 


(dimensionless) 


Hv2p  = — e 
kh 


kh- 


1-j 


v-h  1 

kh 


(dimensionless) 


Magnetic  Vector  Potential  Integral  Coefficient 


BH 


Boh^ 

j  (2  f 


''■'*^(^Hv2-t-kHv2p| 


(Wb/m) 


Magnetic  Vector  Potential  Integral  Phase  Shift  B-and  Other  Coefficients: 


P(9.^.9p,4p)  -cos(0)  cos(0p)  +  sin(0)  sin(0p)  cos(^-  ♦p) 


(dimensionless ) 


C(0p,^p)  -  (sm(8-0p)  cos(^p)^  cos(0p) -I- 8  cos(8-0p)  sin(^p)^  sin(0p)) 


(dimensionless) 


D(0p.4p)  =sin(8-0p)  cos(^)  sm(^p)  cos(0p)-  8  cos(8  0p)  sin(4p)  cos(^)  sin(0p) 
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(dimensionless ) 


Magnetic  Vector  Potential  Integrals 


Ax(0.*)  =-BHpl(- 


a 

2  r2  n 


JO  JO 


khP(9>.ap.%)) 


C(8p.*p)(%»dBp 


(Wb/m) 


fa 

i  2  rl  n 


Av(0,4)  -  BHplI-l  I  I  ei  '‘‘'^^•♦•®‘’-*P*D(0p.^)<%)d0p 
\2/  Jo  Jo 


(Wb/m) 


la\ 

“a 

2  r 

Ui’j 

^0  .( 

Az(0»  =-BH  pl[-j  I  e*  sin(5'0p)  cos(4p)  sm(0p)d4pd0p 


(Wb/m) 


Electric  Field  Components  Ea^E^; 


E0(0,^)  =k^  (Ax(0»  cos(0)  cos(^)  +  Ay(0,^)  cos(0)  sin{4)  -  Az(0,4)  sin(0)) 


(V/m) 


E4(0,^)  =- Ax(0»sin(4)  +  Ay(0.^)cos(^) 


(V/m) 


Radiation  Intensity  U(0) 


-ix2 

U(0,4)  =— •[(  1e0(0,4)!  )\{  |e4(0.*)1  )^1 
2% 


(W/solid  angle) 
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Radiated  Power  Praw: 


r2n 


*  s 

2 


Prad  = 


Jo 


U(e.4)  sm(e)dBd4 

0 


(W) 


Prad  =8.23092*10^ 


(W) 


Directivity  Dq: 


Do 

Prad 


(dimensionless) 


Do  =  2.33662*  10^ 


(dimensionless ) 


Radiation  Resistance  R^: 


Rr  = 


2  Prad 

(|lo|)" 


(D) 


Rr  =  1.64618*10^  (Q) 

Ga in  G : 

G  =  95  Do 

(dimensionless ) 


0=2.21979*10^ 

(dimensionless) 

GdB  =  101og(G)  (dB) 


GdB  =23.46311 


(dB) 
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Effective  Isotropic  Radiated  Power  EIRP: 


EIRP  -Prad  Do 


EIRP  =  1.92325*  10* 


Antenna  Unit  Polarization  Vector  : 


I  I - \ 

I  I  2  .,2 1 


,  Vx^y' 

atan  - - — 

(radians) 

^1  =atanf^ 

\  ^ 

\xj 

0.95532 

(radians) 

^1  =0.7854 

Ex  =E9(0pl.^pl)  cos(0pl)  cos(4pl)  -  E4(0pl,4pl)  sm(4pl) 

Ey  =E0(0pl>pl)  cos(0pl)  sin(^l)  -t-  E^0pl,4pl)  cos(^pl) 

Ez  =-E9(9pl,4pl)  sin(0pl) 


lEx\ 

(ja  r  ^  —  ■  I  Ev 

V(!Ex|  )^^(|Ey|)V(|Ezl)^  \Ez| 


/  0.05346  +  0.40449j 
oa  =  i  0.05361  +  0.40495j 
\-0. 10708  -  0.80945j 


Polarization  Loss  Factor  PLF: 


■/I  —I',  2 

PLF  =ii(jwoal/ 


PLF  =0.16663 


(W) 

(W) 

( radians) 

( radians ) 

(V/m) 

(V/m) 

(V/m) 

(dimensionless ) 

(dimensionless) 

(dimensionless ) 

(dimensionless) 
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Maximum  Effective  Aperture  (A^n^)  : 


X^  Do  95  PLF 

Aem  = - 

4  n 


(m^) 


Aem  =  0.06886 


(m2 ) 


Maximum  Effective  Height  (h^^) : 


hem 


IRr  Aem 


% 


(m) 


hem  =  1  09672 


Aperture  Efficiency  Sapj. 


Aem 

eap  = - 

/dmeas\ 


(dimensionless ) 


cap  =0.08768 


(dimensionless) 
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THI  COMICAL  HORN  AMCTMMA  rAR-TZZLD  tUVATION  PATTKHN8 

For  the  purpose  of  these  far-field  radiation  patterns,  the  conial  horn 
antenna  axis  is  parallel  to  the  Ey  »  0  grid  line  and  the  apex  of  the 
horn  is  located  at  the  origin.  Electric  field  components  behind  the 
horn's  aperture  are  assumed  to  be  zero. 

Theta  Component  of  Elevation  Pattern 

4q  = -  ( radians ) 

2 

Eex(e)  =  |Ee(0.*g)|  cos(0)  E0y(0)  =|E0(O.*g)!  sm(0) 


EFELD 


Phi  Component  of  Elevation  Pattern 


E4x(0)  =  |E4(0,*g)|  cos(0) 


E4y(0)  =  |EK0.*g)i  sin(0) 


E»y(0)-1.O84MO~*® 


-0.002  0  0.003  0.005  0.008  O.Ol  0.012  0.015  0 


—  E  FIELD 


Total  Elevation  Pattern 


Et(0.*) 


HE9(0.*r 


E*(0.*)^ 


Etx(0)  =  iEt(0.*g)|  cos(0)  Ety(0)  =  |Et(0,4g)|  sm(0) 


—  E  FIELD 


Phi  Component  of  Azxmuth  Pattern 


EW*)  nE4(eg.4)|  cos(*)  E*y(*)  =  lE*(eg.*)l  sin(*) 


-0.003  -0.002  -0.002  ~0.00I  -0.001 


0.001  0.001  0.002  0.002  0.003 


E  FIELD 
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0 

£*»<♦) 


Total  Azimuth  Pattern 

Eue.*)  =JmB.4)^^E4(e,4)^ 

Etx(4)  =  |Et(0g.*)|  cos(4)  Ety(*)  HEt(0g.*)|  sinC*) 
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THE  PYRAMIDAL  HORN  ANTENNA 
MATHCAD  SOFTWARE-HORN__PYR .  MOD 

Pyramidal  horn  antennas  are  devices  used  to  provide  a  transistion  from  a 
rectangular  waveguide  to  an  unbounded  medium  such  that  the  wavefront  at 
the  aperture  of  the  horn  has  nearly  a  constant  phase  at  any  point  in  the 
mouth  of  the  horn. 

Pyramidal  horns  are  the  most  popular  type  of  feed  elements  for 
reflectors  used  in  satellite  tracking,  microwave  communications,  and 
radar . 

The  pyramidal  horn  applications  may  be  used  to  analyze  E-  and  H-plane 
sectoral  horns.  To  analyze  an  E-plane  sectoral  horn  set  horn  dimension 
(al)  equal  to  waveguide  dimension  (a) .  To  analyze  an  H-plane  sectoral 
horn  set  horn  dimension  (bl)  equal  to  waveguide  dimension  (b) . 

(Note:  Mathcad  equations  cannot  use  symbolic  subscripts.  Therefore, 

symbols  like  1  will  immediately  follow  the  parameter  in  equations  in 
lieu  of  subscipts.) 

The  pyramidal  horn  antenna  Mathcad  applications  will  compute  the 
following  parameters: 

k  =  Wavenumber 
1  =  Wavelength 

Dq  =  Directivity 

E0  =  Electric  Fiela  (0)  Coitponent 
E|  =  Electric  Field  (4)  Coir^onent 
U  =  Radiation  Intensity 
^max  ~  Maximum  Radiation  Intensity 

Prad  “  Radiated  Power 
G  =  Gain 

EIRP  =  Effective  Isotropic  Radiated  Power 
^em  “  Maximum  Effective  Aperture 
BW  =  Bandwidth 

rmin  =  Minimum  Distance  to  Far-Field 

R^  =  Radiation  Resistance 

^em  “  Maximum  Effective  Height 

^x,y, z  ~  Electric  Field  Components  in  Cartesian  Coordinates 
Og  =  Unit  Polarization  Vector 

fc  =  Transverse  Electric  and  Transverse  Magnetic  Cutoff  Frequencies 

Egp  =  Aperture  Efficiency 

PLF  =  Polarization  Loss  Factor 

Dp  =  Pyramidal  Horn  Corner  to  Corner  Distance 

Pe,h  ~  Pyramidal  Horn  Perpendicular  Flare  to  Mouth  Distances 
Il,2  =  Electric  Field  Component  Equation  Coefficients 
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The  following  data  must  be  input  based  on  known  or  estimated  data: 

t  =  Number  of  Cutoff  Frequencies  Calculated 
m, n  =  Modes 

f  =  Frequency  of  Interest 

i  =  Number  of  Increments  for  Far  Field  Radiation  Patterns 
rff  =  Distance  of  Far-Field  Calculations 

Iq  =  Input  Current  at  Antenna  Terminals 

6g  =  Coaltitude  (Deflection  Angle  from  +z  Axis)  for  Azimuth  Plot 
=  Azimuth  Angle  for  Elevation  Plot 

(x,y, z)  =  Coordinates  for  Unit  Polarization  Vector 
a,b  =  Rectangular  Waveguide  Dimensions 
ai,bi  =  Pyramidal  Horn  Dimensions 

Pl>P2'Pe'Ph  ~  Pyramidal  Horn  Imaginary  Cone  Apex  to  Mouth  Distances 
Eq  =  Electric  Field  Amplitude  Constant 
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Enter  input  data  here; 


f  =9.3  10®  (HZ) 


t  =5  (modes) 

lo  =  1  ( an^s ) 


i  =36  (increments)  rff  =  1  10^  (meters) 

(Mote:  For  E-plan«  sectoral  horn  analysis 
sat  (al)  equal  to  (a) .  For  B-plane  sectoral 
horn  analysis  set  (bl)  equal  to  (a)) 


a  = .02286  (m) 

al  =1846  (m) 

b  =  01016  (m) 

bl  =  1455  (m) 

pe  =.3281  (m) 

ph  =  .3521  (m) 


(JW 


pi  =  3398 


p2  =.3198 

X  =10^ 


y  =10^ 


z 


=  10^ 


(dimensionless ) 

(m) 

(m) 

(m) 

(m) 


Eo  =1 


(V/m) 


(m) 


ramidal  horn  antenna  geometric  parameters  and  define 

constants : 

c  =2.9979  10*  (meters/sec) 

X  =-  (meters/cycle) 

f 

X  =0.03224  (meters/cycle) 


Dp  =Jal^+b]^ 

(m) 

Dp  =0.23505 

(m) 

pe  =(bl  -  b)  [—1  -  - 

(m) 

ph  =(al -  a) 

IS’-i 

(m) 

p€  =  0.29759 

(m) 

ph  =  0.29771 

(m) 

Calculate 
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Wavenumber  )c: 


(m 


k  =  1  94915*10^ 


Cutoff  Frequencies 


n  =0  .  t 


(modes) 


m  =0  .  t 


(modes) 


fc 


(HZ) 


r 

0 

1.47638-10“’ 

6.56168- 10® 

1  61563-10‘® 

1.31234-10“’ 

1.97533-10“’ 

l.%85-10“’ 

2.46063- 10‘® 

2.62467-10“’ 

3.01141-10“’ 

1  3.28084-10“’ 

3.59772-10’'’ 

2.95276-10*'’ 

4.42913-10*'’ 

3.02478-10*'’ 

4.47748-10*'’ 

3.23125-10*'’ 

4.61946-10*'’ 

3.54877-10*'’ 

4.84688-10*'’ 

3.95065-10*'’ 

5,14841-10*'’ 

4.41392-10*'’ 

5.51191-10*'’ 

5.90551-10*'’ 

7.38189-10*'’  j 

5.94185-10*'’ 

7  411-10*'’ 

6.04957-10*'’ 

7.49763-10*'’  ^ 

6.22496-10*'’ 

7,63985- 16”*^ ’’ 

6.46251-10*'’ 

7.83462-10*'’  ; 

6  75566-10*'’ 

8.07813-10*'’ 

(Note:  The  index  for  both  m, n  above  begins  with  zero.  TM  modes  cannot 
have  m  or  n  equal  zero.  TEqo  mode  does  not  physically  exist.) 
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Electric  Field  Component  Coefficients  It, I?: 

V  =0.1 

w  =0.1 


(increments) 


(increments ) 


I'al 

2 


lue.*)  - 


I 

J  2 


-j  t 

e 


_5I —  iio(8)'C<a(t)  4 
l2p2  / 


(dimensionless) 


Ill 


al 

T 


/  4 

COSi  It  — 
\  2 


al 

Y 


ji-. 

2p2 


'  V  _  *>  ‘  W 

"ir‘-2j'“‘iT 


•2-.  -4 


cosik  -  sm(0)  cos(^)j 

a 

If-  _  • 

\  2  / 

.ii(e.*) 

rt  •" 

2 

|k-^'sin(0)-cos(^)|  -  bj 

(dimensionless ) 


(dimensionless) 


III 


J 

=  if  al=a,-Jt  - 


I,  a  .  Iv  n  ,--6)  /w,  ,--6 

cx)s|k— sin|— 11 — -t- 10  I  cosj — 2  it -i- 10 

\i 


,  a  .  /V  It  ,_-6)  /w  -6 

k-  sm(-  it  —  +10  j  cosI— •2-it+ 10 


b] 

\2/ 


.Ill  I 

V,W  i 


(dimensionless) 


f— 

2 


I2(0.*)  = 

J 


2 


-  •in(6)  tin(4)-4 


(dimensionless) 
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121 


V,  W 


I  2 

j 

i  y 
J  2 


I  i  sui|k  -  sai(6)  sin(4)j  ! 

12(0, ♦)  ^il|bl»b.bi  '  ^ - _L|.I2(0,i) 

'  1  (k- sm(e)sin(4)i  ! 

;  2  ;  ! 


121 


=  if  bl«b.b 


i[k--'Sin| --Ji- 

1.10^ 

sm(—  2  Ji-t-  10'**] 

\  2  \i 

2  / 

\i  } 

/.  b  . 


I  “  I  V  \fs-6\  /w-  ,rt-6 

k  -  sm  -  n — 1-10  sin — 2•s-^  10 


■il 


.121. 


Electric  Field  Components  Eq, E^: 


■j  krff 

E0(0,*)  =j  kEo - (sin(*)(Ucos(9))Il(9,*)I2(0.*)) 

4  n  rflf 


-j  krff 

E01  =j  k-Eo - 

4)1  rff 


sin[— 


/  V 

COSI  -  Ji¬ 
ll 


A\ 


Ill  121 

v,w  v.wj 


-j  krff 

E*(0.*)  -J  kEo^ - (cos(*)  (l^cos(0))  Il(0.*)  I2(9.4)) 

4  Ji  rflf 


I 

V,  W 


=  j  k  Eo  - 


■  j  k  rff 


cos(—  2  iii  (  1  -h  cos/-  ji-  j-Ill 


121 


V,W 


(dimensionless ) 


(dimensionless ) 


(dimensionless) 


(V/m) 


(V/m) 


(V/m) 


(V/m) 
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(W/solid  angle) 


Radiation  Intensity  U (8)  : 


U(0.*)  =  —  [(  !E0(0.*)1  ( 1E*(0.*)|  )^j 

2% 


U1  zi??-!nE01 


(W/solia  angle) 


Umax  =max(Ul) 


(W/solid  angle) 


Radiated  Power  P, 


•2-n  i  - 


Prad  =  ;  U(0,*)  sin(0)<10d4 

JO  Jo 


Prad  =  3.39938- 10 


Directivity  D,-,: 


4  n  Umax 


(dimensionless ) 


Do  =  1. 47304- lO'^ 


(dimensionless ) 


Radiation  Resistance  R, 


2  Prad 


(I  loir 


Rr  =  6.79875- 10 
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Gain  G: 


G  =  5  Do 


0=73.65177 


GdB  =101og(G) 


GdB  =  18.67183 


Effective  Isotropic  Radiated  Power  EIRP: 

EIRP  =PradE)o 

EIRP  =  5.0074M0"^ 


Antenna  Unit  Polarization  Vector  : 


Op  -atanl 

\ 


(radians) 


♦p 


Op  =0.95532 


(radians) 


♦P  =0.7854 


Ex  =E0(ep,*p)  cos(0p)  cos(^)-  E4(0p,(Jp)  Sin(^) 


Ey  =E9{ep,*p)  cos(0p)  sm(lp)  +  E4(0p,4p)  cos(4p) 


Ez  =  -E9(ep>p)  sin(0p) 

1 

oa  = - 

v'dEx!  )^^(|Ey|)V(|Ez!)^  \Ez, 

/-0.18825  +  0.09603J  ^ 
oa  =  I  0.70255  -  0.35837j 
\-0.5143  +  0.26235J  ! 


(dimensionless ) 


(dimensionless ) 


(dB) 


(dB) 


(W) 


(W) 


(radians) 


(radians) 


(V/m) 


(V/m) 


(V/m) 


(dimensionless ) 


(dimensionless ) 


Polarization  Loss  Factor  PLF: 


PLF 


(  <rw  oa 


2 


(dimensionless ) 


PLF  =  0.33333 


(dimensionless ) 


Maximum  Effective  Aperture 


Aem 


i^  Do  .5  PLF 
4  n 


(m2) 


Aem  =2.0301 1*10 


(m2 ) 


Maximum  Effective  Height 


hem 


l' 

jRr  Aem 

■i~ 


(m) 


hem  =  3.82683*  10  * 


(m) 


Aperture  Efficiency 


cap 


Aem 
al  bl 


(dimensionless ) 


tap  =0.07558 

(dimensionless 
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THS  PTRMCIDAL  BORN 


rAR-rXXLO  AZZMDTH  PATTBRMS 


For  the  purpose  of  these  far-field  radiation  patterns,  the  pyramidal 
horn  central  axis  is  perpendicular  to  the  Ey  =  0  and  Ex  =  0  grid  lines. 

Theta  Con^onent  of  Azimuth  Pattern 


6g  (Note:  6_  must  be  less  than  a/2  radians) 

10  ^ 


Eex(*)  =  lEe(0g.4)|  cos(*)  E0y(4)  =  |E0(eg.*)|  sin(*) 
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